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Nuclear  DNA  PCR-RFLP  markers  were  sought  and  used  to  distinguish  African 

and  European  honey  bee  (Apis  mellifera  L.)  subspecies.  These  markers  along  with 

previously  described  markers  were  used  to  genetically  characterize  and  to  determine 

levels  of  hybridization  and  gene  introgression  between  New  World  African  and  European 

honey  bee  populations.  Markers  were  sought  in  RFLP  analyses,  by  using  cloned  probes 

in  Southern  blots,  in  RFLP  analysis  of  long  (>5kbp)  and  short  (<2kbp)  DNA  regions 

amplifiable  by  long  and  standard  PCR,  and  in  amplified  fragment  length  polymorphism 

(AFLP)  analysis.  To  facilitate  the  analyses,  markers  found  using  Southern  blots  and  long 

PCR  procedures  were  made  analyzable  using  standard  PCR  formats  by  amplifying  the 

regions  containing  the  informative  polymorphic  sites.  Maps  for  each  polymorphic  locus 

are  presented.  Characterization  of  markers  specific  and  predominant  to  African  and 

European  honey  bee  subspecies  is  described. 
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New  World  honey  bee  populations  from  Venezuela,  Costa  Rica,  Honduras, 
Mexico  and  the  United  States  (including  two  populations  sampled  at  two  different  years 
near  the  subtropical-temperate  boundary  in  Texas),  representing  different  lengths  of  time 
since  African  bee  invasion,  and  Old  World  European  and  African  populations  were 
genetically  characterized.  Levels  of  hybridization  between  New  World  African  and 
European  populations  were  determined. 

Significant  differences  in  the  genetic  structure  were  found  between  Old  World 
European  and  New  World  African  and  European  populations  (P<0.01)  but  not  between 
New  World  and  Old  World  African  populations  (P>0.05).  From  a  north  to  south 
direction,  in  neotropical  populations,  an  increase  in  average  frequencies  of  African- 
specific  alleles  was  seen,  whereas  frequencies  of  east  European-specific  alleles  decreased. 
West  European  allele  frequencies,  however,  remained  constant,  suggesting  a  high  level  of 
population  admixture  between  African  and  west  European  bees.  Low  levels  of 
hybridization  between  African  and  European  bees  were  found  in  neotropical  populations 
and  most  hybrids  were  of  the  African-west  European  type.  In  the  subtropical-temperate 
region  in  Texas,  an  increase  in  average  frequency  of  African-specific  alleles  was  seen  as 
reported  in  the  neotropics,  however,  a  high  level  of  hybridization  between  African  and 
European  bees  was  found. 
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CHAPTER  1 
INTRODUCTION 

Bioeeographv  and  Distribution  of  Honey  Bee  Subspecies  in  the  World. 

The  honey  bee,  Apis  mellifera  L.,  is  native  to  the  Old  World  and  was  originally 
distributed  in  Europe,  Africa,  and  the  Middle-East.  The  wide  geographical  distribution 
and  adaptation  to  different  environments  has  resulted  in  the  evolution  of  24  subspecies,  or 
races,  initially  characterized  and  described  using  multivariate  analysis  of  morphological 
characters  (Ruttner,  1988).  Four  taxonomic  lineages  (or  branches)  have  evolved:  near- 
east,  tropical  sub-Sahara-African,  north  African-west  Mediterranean  (and  northwest 
European)  and  central  Mediterranean-southeast  European  (Ruttner,  1978).  The  near-east 
branch  includes  the  subspecies  A.m.  anatolica,  A.m.  adami,  A.m.  cypria,  A.m.  syriaca, 
A.m.caucasica,  and  A.m.  meda.  The  tropical  sub-Sahara-African  group  includes  the 
subspecies  A.m.adansonii,  A.m.  litorea,A.m.  monticola,  A.m.  lamarckii,  A.m.  capensis, 
A.m.  unicolor,  and  A.m.  scutellata.  The  north  African  west-Mediteranean  group  includes 
the  subspecies  A.m.  sahariensis  and  A.m.  intermissa,  from  north  Africa,  and  A.m. 
mellifera  and  A.m.  iberica  from  western  and  northern  Europe.  The  central 
Meditteranean-southeast  European  group  includes  the  subspecies  A.m.  sicula,  A.m. 
ligustica,  A.m.  carnica,A.m.  macedonica,  and  Am.  cecropia.  Throughout  this 
dissertation,  "African  bees"  will  refer  to  the  subspecies  A.m.  scutellata,  "west-European 
bees"  will  refer  to  the  subspecies  A.m.  mellifera  and  A.m.  iberica,  "east-European  bees" 
will  refer  to  the  subspecies  A.m.  ligustica  and  A.m.  carnica,  and  "near-eastern  bees"  will 
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refer  to  the  subspecies  A.m.  caucasica.  Samples  from  the  other  African  and  European 
subspecies  were  not  available  for  this  study. 
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Importation  of  Honev  Bee  Subspecies  and  Establishment  in  the  New  World. 

The  first  honey  bee  colonies  were  brought  to  the  New  World  by  early  European 
settlers  during  the  first  half  of  the  seventeenth  century  and  were  primarily  of  the 
subspecies  A.m.  mellifera,  called  German  brown  or  north  European  black  bees  (Oertel, 
1976).  Colonies  of  A.m.  mellifera  were  also  imported  to  North  America  by  the  English  in 
1763  and  by  the  Russians  in  1809,  but  it  is  not  known  if  these  colonies  survived  (Nelson, 
1967).  In  1855,  an  attempt  to  introduce  Italian  bees  {A.m.  ligustica)  to  North  America 
was  made,  however,  not  until  1860  was  the  first  successful  introduction  of  this  race  in 
North  America  was  reported  (Pellet,  1938).  Other  honey  bee  subspecies  from  Hungary 
and  Tunesia  (near-east  group)  and  from  north  Africa  followed  but  with  little  success, 
probably  due  to  the  poor  adaptability  of  these  bees  to  temperate  climates  and, 
consequently,  their  low  performance  as  honey  producers  (Pellet,  1938).  Egyptian  (A.m. 
lamarckii),  Cyprian  (A.m.  cypria)  and  Caucasian  (A.m.  caucasica)  bees  were  also 
introduced  between  1 865  and  1 905.  Of  these,  only  the  Caucasian  bees  were  used  to  any 
extent  by  beekeepers.  Recent  mitochondrial  DNA  evidence  has  been  reported  indicating 
that  honey  bees  from  the  African  subspecies  group,  specifically  A.m.lamarckii,  had 
become  established  at  a  low  frequency  (1 .0%)  in  North  America  (Schiff  and  Sheppard, 
1993).  In  North  America,  the  most  popular  subspecies  among  beekeepers  has  become 
A.m.  ligustica,  Italian  bees,  followed  by  A.m.carnica,  Carniolan  bees,  which  are  gentler 


and  better  honey  producers  than  A.m.  mellifera.  In  temperate  regions  of  North  America, 
European  bees  established  a  substantial  feral  population  composed  primarily  of 
A.m.mellifera,  A.m.  lignstica,  and  A.m.carnica. 

In  the  neotropics,  the  first  honey  bees,  A.m.  mellifera  and  possibly  A.m.iberica, 
were  introduced  into  Cuba  from  Spain  in  1763  and,  later,  into  Central  America.  A.m. 
mellifera  was  also  brought  to  Brazil  in  1839  followed  by  A.m.  ligustica  in  1922  (Kent, 
1988).  A.m.  iberica  was  also  brought  by  the  Portuguese  to  Brazil  during  the  same  period 
(Michener  et  al,  1972).  Unlike  in  North  America,  A.m.  mellifera  remained  the 
predominant  race  among  beekeepers  in  the  neotropics  (Kent,  1988).  European  bees  were 
poorly  adapted  to  the  tropics  and  did  not  established  self-sustaining  feral  populations. 

African  bees,  A.m.  scutellata,  were  imported  to  Brazil  in  1956  with  the  intention 
of  interbreeding  with  and  improving  the  productivity  of  the  European  bees.  In  1957, 
swarms  from  26  African  colonies  escaped  from  apiaries  (Kerr,  1967),  which  resulted  in 
the  dramatic  spread  of  African  bees  in  the  New  World.  The  African  bee  migration  has 
been  an  impressive  biological  phenomenon  in  which  one  subspecies  with  greater 
ecological  advantages  has  taken  over  an  existent,  albeit  largely  managed,  population  of 
other  subspecies.  Resident  European  bees  were  quickly  replaced  by  Africanized  progeny 
as  a  result  of  paternal  introgression,  that  is,  African  drones  mating  with  European  queens 
(Michener,  1975),  and  by  nest  usurpation  or  invasion  of  managed  European  colonies  by 
African  swarms  (Danka  et  al,  1992;  Vergara  and  Perez  de  Leon,  1993).  As  the  African 
bees  expanded  their  range,  managed  colonies  expressed,  over  time,  a  significant 
reduction  of  morphometric,  behavioral,  and  genetic  characteristics  typical  of  European 
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bees  and  an  increase  in  African  traits.  (Quezada-Euan  and  Medina,  1998;  Quezada-Euan 
andPaxton,  1999;  Michener,  1975). 

African  bees  first  entered  the  United  States  in  October  1990  in  Hidalgo,  Texas 
(Hunter  et  al,  1993).  The  northernmost  limit  of  the  African  honey  bee  spread  in  the 
United  States  has  been  predicted  on  the  basis  of  a  10°C  isotherm  line,  corresponding  to 
the  average  high  temperature  in  mid  winter,  based  on  their  actual  distribution  in  South 
America  (Taylor  1977;  Taylor  and  Spivak  1984).  They  had  been  predicted  to  colonize 
most  of  the  southern  states  of  the  United  States  by  1995,  but,  currently,  they  are  only  in 
California,  Arizona,  Nevada,  New  Mexico,  and  Texas.  Factors  such  as  the  high  density 
of  European  colonies  in  Texas  and  the  effect  of  the  parasitic  mite  Varroa  jacobsonii  are 
believed  to  be  responsible  for  the  unexpected  stall  of  the  expanding  African  front  toward 
the  southeastern  states  (Sanford,  1994;  Eischen,  1999).  African  bee  movement,  however, 
has  recently  resumed  an  easterly  direction,  apparently  due  to  the  decimation  of  feral 
European  colonies  by  the  Varroa  mite  in  Texas  (Eischen,  1999).  African  swarms  have 
been  detected  near  Florida  ports  from  accidental  introductions  from  ships,  but  it  is  not 
known  if  such  introductions  will  result  in  the  establishment  of  a  substantial  feral  African 
bee  population  (Sanford,  1999).  In  Puerto  Rico,  however,  a  self-sustaining  African  feral 
population  exists  that  was  established  as  a  result  of  swarms  arriving  on  ships. 

Characteristics  of  African  and  European  Honey  Bees. 


African  bees  evolved  in  tropical  regions,  whereas  European  bees  evolved  in 
temperate  regions.  Unpredictable  rains,  as  well  as  pollen  and  nectar  availability,  have 


played  an  important  role  in  the  evolution  of  African  honey  bee  subspecies.  Most  of  the 
subspecies  variation,  based  on  morphological  features,  is  found  in  Africa  and  has  been 
attributed  to  the  diversity  of  ecosystems  and  the  variability  of  rainfall  (Ruttner,  1998; 
Rinderer  and  Hellmich,  1988). 

In  contrast  to  African  bees,  European  bees  are  adapted  to  more  predictable 
climatic  regimes,  evolving  so  that  their  biological  cycles  of  brood  rearing  and  food 
storage  coincide  with  temperate  seasonal  changes  (Rinderer  and  Hellmich,  1988).  In 
general,  European  bees  have  better  overwintering  capabilities,  lower  swarming  and 
absconding  rates  and  are  less  defensive  than  African  bees  (Michener  et  al,  1972).  The 
factors  responsible  for  good  overwintering  capabilities  of  the  European  bees  include  the 
large  size  of  the  individual  bees,  effective  thermoregulation  of  the  brood  nest  and 
economical  storage  and  consumption  of  honey  (Winston  et  al.,  1983). 

Swarming  (natural  division  of  a  colony)  by  African  bees  is  five  times  more 
frequent  than  by  European  bees.  Absconding  (abandoning  of  a  nest  by  a  colony)  occurs 
infrequently  by  European  subspecies  but  commonly  by  African  subspecies,  probably  due 
to  a  higher  predation  rate  of  African  bees  in  the  tropics  (Winston  et  al,  1983).  Other 
differences  between  African  and  European  honey  bees  include  the  following: 
developmental  period  (18  days  for  African  workers  compared  to  21  days  for  European 
workers;  Harbo  et  al.,  1981),  adult  lifespan  (shorter  in  African  bees;  Winston  et  al, 
1983),  selection  and  size  of  nesting  sites  (African  bees  are  usually  less  selective  and 
usually  nest  in  smaller  cavities;  Rinderer  et  al.,  1982),  defensive  behavior  (African  bees 
are  more  defensive;  Stort  1974;  Michener  et  al.,  1972  and  Michener,  1975),  size  of  the 


swarms  (African  bees  have  smaller  swarms;  Rinderer  et  al,  1982)  and  collection  of 
nectar  (African  bees  collect  nectar  of  lower  sugar  concentration;  Rinderer  et  al,  1982). 

High  swarming,  absconding,  and  defensive  behavior  make  African  bees  more 
difficult  to  manage  for  commercial  beekeeping  operations  (Michener,  1975).  In  the 
southern  United  States,  they  have  already  been  responsible  for  several  human  and 
domestic  animal  fatalities.  In  Florida,  they  are  expected  to  have  a  negative  impact  on  the 
tourism  and  beekeeping  industries.  Stinging  incidents  and  even  deaths  related  to  African 
bee  attacks  will  increase,  along  with  liable  lawsuits  that  result  from  these  attacks. 
Movement  of  bees  from  Florida  to  more  northern  states,  a  common  practice  involving 
thousand  of  hives,  will  be  restricted,  resulting  in  economic  losses  for  operations  that 
largely  depend  on  migratory  beekeeping. 

Methods  Used  to  Identify  African  and  European  Honey  Bees. 

African  and  European  honey  bees  have  been  discriminated  based  on  behavioral 
(Collins  et  al,  1982),  morphometric  (Daly  and  Balling,  1978;  Daly  et  al,  1982), 
biochemical  (cuticular  hydrocarbons;  Carlson  and  Bolten,  1983;  Lavine  et  al,  1988; 
Smith,  1990)  and  genetic  characteristics  (allozymes  and  DNA;  Hall,  1986;  Hall,  1992b; 
Hall  and  Smith,  1991;  Sheppard  and  Berlocher,  1985;  Smith  and  Brown,  1988;  Sylvester, 
1982).  Accurate  identification  is  important  for  regulatory  purposes,  certification  of 
honey  bee  stocks,  and  research  (Page  and  Erickson,  1985).  Reliable  identification  will  be 
needed  to  monitor  the  dynamics  of  the  expanding  African  population  in  the  United  States, 


for  example,  to  determine  the  extent  to  which  African  bee  characteristics  introgress  north 
of  the  expected  hybrid  zone  in  North  America. 

One  of  the  most  notorious  behavioral  characteristics  that  distinguish  African  and 
European  bees  is  their  high  defensiveness  (Michener,  1975).  Defensiveness  has  been 
evaluated  using  quantitative  measures  of  variables  such  as  the  time  at  which  the  first  sting 
occurred,  time  taken  for  the  colony  to  become  "very  defensive"  (African  bees  have  faster 
reaction  time),  number  of  stings  (African  bees  sting  in  greater  numbers),  and  distance  that 
the  bees  follow  an  intruder  (Collins  and  Kubasek,  1982;  Collins  et  al,  1982;  Stort,  1974). 
Another  characteristic  used  for  identification  is  the  comb  size,  which  is  usually  smaller  in 
African  bees  (Spivak,  1988).  Although  behavioral  characteristics  are  commonly  used, 
they  are  subject  to  variation  from  environmental  factors  and  management  techniques,  and 
they  cannot  be  reliably  used  to  distinguish  hybrids. 

Morphometric  analysis  has  also  been  used  to  discriminate  honey  bee  groups.  It 
measures  both  genetically  based  and  environmentally  induced  variation  using  one 
(univariate)  or  several  variables  (multivariate  analysis)  (Daly,  1988).  Multivariate 
analysis  uses  the  measurement  of  different  morphological  characters  to  obtain  a 
probability  value  for  a  particular  sample  being  African  or  European.  Such  analysis  has 
been  employed  in  the  Fast  Africanized  Bee  Identification  System  (FABIS)  (Daly  and 
Balling,  1978;  Rinderer  et  al.,  1987b).  Multivariate  analysis  requires  sophisticated 
equipment  and  trained  personnel.  Multivariate  analysis  assigns  a  probability  value  of 
being  African  or  European  (Daly  et  al,  1982).  Errors  can  result  from  the  effect  of 
environmental  factors  on  the  variation  of  measurements  of  morphological  characters.  As 
with  behavioral  methods,  morphometries  cannot  be  used  effectively  to  identify  hybrids 
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and,  therefore,  is  not  effective  in  studying  gene  introgression  and  hybridization. 
Nevertheless,  morphometries  has  been  successfully  used  to  follow  changes 
accompanying  the  Africanization  of  some  populations  in  the  New  World  (Boreham  and 
Roubik,  1987).  Morphometric  values  have  suggested  a  European  component  in  the 
neotropical  African  population  (Quezada-Euan  and  Hinsull,  1995). 

One  biochemical  method  is  based  on  the  composition  of  the  insect  cuticle 
(Lockey,  1979).  Hydrocarbons  are  extracted  by  immersing  the  insect  body  in  a  solvent 
such  as  hexane,  and  the  extracted  compounds  are  analyzed  using  gas  chromatography 
(GC).  Differences  in  the  hydrocarbon  composition  of  African  and  European  bees  have 
been  found  (Carlson  and  Bolten,  1983;  Smith,  1990).  The  cuticular  hydrocarbon 
composition,  however,  is  also  affected  by  environmental  conditions.  Bees  of  the  same 
subspecies  that  feed  on  a  particular  diet  can  have  a  different  hydrocarbon  composition 
from  bees  that  feed  on  a  different  diet.  Contamination  from  pollen,  propolis,  and  plant 
extracts  can  obscure  the  results.  Furthermore,  the  inheritance  pattern  of  the  compounds  is 
not  known. 

Genetic  differences  between  African  and  European  honey  bees  have  been  found 
in  the  form  of  protein  markers.  Protein  polymorphism,  expressed  as  allozymes,  has  been 
widely  used  because  it  is  detected  by  relatively  inexpensive  and  fast  procedures. 
Allozymes  are  inherited  in  a  clear  Mendelian  fashion  (Del  Lama  et  ai,  1990).  Allozyme 
identity  in  individuals  is  independent  of  environmental  influences,  being  the  direct 
product  of  structural  genes  (Daly,  1988).  However  the  existence  of  protein  variants  is 
subject  to  selection,  and  their  detection  depends  on  their  expression  (Hall,  1988).  Honey 
bees  and  Hymenopteran  insects  in  general  show  low  levels  of  allozyme  polymorphism 


(Crozier,  1970;  Metcalfe/  al,  1975;  Sheppard  and  Berlocher,  1985;  Snyder,  1974).  This 
low  variability  in  social  Hymenoptera  may  be  due  to  the  exposure  of  deleterious  alleles  in 
haploid  males  and  the  relatively  small  size  of  the  reproductive  populations  (Graur,  1985; 
Metcalf  et  al,  1975).  Some  of  the  polymorphic  proteins  currently  available  for  honey 
bees  include  the  enzymes  esterase,  hexokinase,  alcohol  dehydrogenase,  malate 
dehydrogenase,  and  a  non  enzymatic  protein  P-3  (Contel  et  al,  1977;  Mestriner,  1969; 
Mestriner  and  Contel,  1972;  Sheppard  and  Berlocher,  1985). 

Unlike  protein  polymorphism,  DNA  variation  is  not  necessarily  subject  to 
environmental  selective  pressures,  and  its  detection  does  not  depend  on  gene  expression. 
Therefore,  DNA  analysis  is  the  most  direct  method  to  establish  identities  and  relatedness 
among  organisms  (Hall,  1988).  DNA  analyses  include  both  mitochondrial  (mtDNA)  and 
nuclear  (nDNA)  DNA.  In  multicellular  animals,  mtDNA  is  maternally  inherited  (without 
recombination):  the  entire  mtDNA  is  transmitted  from  generation  to  generation.  Changes 
in  mtDNA  sequences  reflect  accumulated  mutations  (Brown,  1985).  Because  mtDNA  is 
uniparentally  inherited,  it  cannot  be  used  to  study  paternal  gene  introgression  and 
hybridization.  Another  disadvantage  of  mtDNA  is  the  relatively  small  size  of  the 
molecule,  which  limits  the  number  of  polymorphisms  that  can  be  found.  Nevertheless, 
mitochondrial  DNA  restriction  fragment  length  polymorphisms  (RFLPs)  have  been  found 
that  distinguish  African  and  European  honey  bees  (Hall  and  Smith,  1991;  Smith  and 
Brown,  1988).  These  markers  have  been  used  to  test  neotropical  feral  honey  bee  swarms, 
thereby  demonstrating  that  the  spread  of  the  African  bees  in  the  New  World  has  been 
through  unbroken  African  maternal  lineages  (Hall  and  Muralidharan,  1989;  Hall  and 
Smith,  1991;  Smith  et  al,  1989). 


10 

Compared  to  mtDNA,  nDNA  is  vastly  more  abundant  in  multicellular  organisms, 
allowing  for  the  screening  of  a  practically  unlimited  number  of  polymorphisms.  Nuclear 
DNA  markers,  being  maternally  and  paternally  inherited,  enable  hybridization  and  gene 
introgression  studies.  Nuclear  DNA  RFLP  markers  are  codominant,  thus,  allowing  the 
distinction  between  homozygous  and  heterozygous  individuals.  With  anonymous  cloned 
probes  in  conventional  Southern  blot  procedures  (Southern,  1975),  nDNA  RFLP  markers 
that  distinguish  African  and  European  bees  were  found  and  used  in  the  African  bee 
population  studies  described  above  (Hall,  1986;  Hall,  1990;  Hall,  1992b).  Southern  blot 
procedures,  however,  are  laborious,  expensive,  require  considerable  technical  expertise 
and,  thus,  limit  the  application  of  RFLP  markers  for  population  genetic  studies  in  which 
large  sample  sizes  are  required. 

Other  nuclear  DNA  markers  commonly  used  to  distinguish  honey  bee  subspecies 
include  random  amplified  polymorphic  DNA  (RAPD)  (Williams  et  al,  1990).  The 
simplicity  of  the  RAPD  analysis  makes  these  markers  popular  for  population  genetic 
studies.  However,  RAPD  markers  (Williams  et  al.,  1990)  are  generally  dominant,  that  is 
homozygous  and  heterozygous  individuals  cannot  be  distinguished.  Other  problems 
associated  with  the  reliability  and  reproducibility  of  these  markers  further  limit  their 
application  (Black,  1993;  MacPherson  et  al,  1993;  Schierwater  and  Ender,  1993; 
Weeden,  1992).  RAPD  markers  have  been  successfully  used  to  map  the  honey  bee 
genome  (Hunt  and  Page,  1992;  Hunt  and  Page,  1995).  RAPD  markers  that  distinguish 
African  and  European  bees  have  been  reported  (Suazo  et  al.,  1998). 

Microsatellite  markers  from  single  and  multiple  loci  specific  to  east  European, 
west  European,  and  African  bees  have  been  used  to  characterize  Old  World  European  and 
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African  honey  bee  populations  (Estoup  et  al,  1993;  Estoup  et  al,  1995).  Microsatellite 
analyses  are  technically  more  difficult,  but  simple  banding  patterns,  easy  to  interpret  are 
produced.  However,  some  microsatellite  DNA  has  been  found  to  be  evolutionary 
unstable  (Tautz,  1989). 

In  past  studies,  the  technical  difficulties  of  RFLP  analysis  through  Southern  blot 
procedures  limited  the  number  and  specificity  of  markers  that  could  be  found  and  limited 
the  number  of  samples  that  could  be  tested.  Most  seriously,  at  the  time  of  the  first 
studies,  no  markers  specific  to  west  European  {A.m.mellifera  or  A.m.iberica)  had  been 
found  (Hall,  1990)  The  predominant  European  race  in  Brazil  at  the  time  of  African  bee 
introduction  was  A.m.mellifera,  but  any  contribution  of  this  race  to  the  neotropical 
African  population  could  not  have  been  detected.  Subsequently,  in  Southern  blots,  RFLP 
alleles  were  found  at  a  highly  polymorphic  locus  characteristic  of  west  European  bees 
(McMichael  and  Hall,  1996).  Because  the  type  of  analysis  limited  the  number  of  samples 
that  could  be  tested  and  because  very  complex  band  patterns  characterized  the  alleles, 
accurate  quantification  of  allele  frequencies  was  precluded. 

To  make  RFLP  markers  more  productive  to  use,  the  polymerase  chain  reaction 
(PCR)  (Saiki  et  al.,  1988)  has  been  employed.  Small  DNA  regions  are  amplified  and 
then  digested  with  restriction  endonucleases  to  reveal  polymorphisms.  This  approach  is 
called  PCR-RFLP.  Compared  to  Southern  blot  procedures,  PCR-RFLP  requires  less 
DNA  and  no  radioactive  labeling.  It  is  relatively  fast  and  inexpensive.  Large  sample 
sizes  can  be  processed  in  less  time.  Thus,  PCR-RFLP  can  greatly  facilitate  population 
genetic  studies.  For  honey  bees,  PCR-RFLP  was  first  used  to  make  their  mtDNA  easier 
to  analyze  (Hall  and  Smith,  1991;  Crozier  et  al,  1991).  More  recently,  honey  bee  nuclear 
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DNA  RFLPs,  initially  found  with  cloned  probes,  have  been  made  PCR  analyzable:  one 
locus  with  African-specific  markers  (Hall,  1998)  and  three  loci  with  east  European 
markers  (Hall,  in  preparation). 

Processes  Involved  in  the  Spread  of  the  African  Bees  in  the  New  World. 

In  the  forty  three  years  since  the  African  bees  were  first  released,  they  have 
successfully  established  self-sustaining  feral  populations  in  the  tropical  and  subtropical 
regions  of  the  New  World.  The  processes  involved  in  the  spread  and  establishment  of  the 
African  bees  are  not  well  understood.  There  are  two  controversial  views  (reviewed  in 
Hall,  1992a).  One  view  proposes  that  the  primary  driving  force  responsible  for  the 
spread  has  been  through  African  drones  mating  with  European  queens,  resulting  in  a 
population  called  a  "hybrid  swarm"  (Rinderer  et  al,  1985).  These  hybrid  colonies,  which 
were  originally  European  and  have  been  called  "Africanized  bees",  further  perpetuated 
the  spread  of  African  genes.  According  to  this  view,  differential  drone  production  and 
social  parasitism  are  the  important  factors  that  have  contributed  to  the  Africanization 
process  (Rinderer  et  al,  1985;  Rinderer  et  al,  1986;  Rinderer  et  al,  1987a).  The  other 
view  is  that  the  primary  driving  force  has  been  African  maternal  expansion  through  feral 
swarms.  This  view  was  based  on  strong  mitochondrial  (mtDNA)  evidence  from 
independent  studies  (Hall  and  Muralidharan,  1989;  Hall  and  Smith,  1991;  Smith  et  al, 
1989).  Virtually  all  mtDNA  in  feral  swarms  in  Brazil,  Venezuela,  and  Mexico  were  of 
the  African  type.  This  view  also  holds  that  the  neotropical  population  is  not  a  hybrid 
swarm.  Hall  (1990)  used  both  mtDNA  and  east  European  nDNA  markers  to  test  a  group 
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of  managed  European  colonies  in  southern  Mexico  established  before,  but  sampled  after, 
African  bees  entered  the  region.  The  results  showed  that  the  colonies  retained  European 
mtDNA  but  had  low  levels  of  European  nDNA  markers,  probably  the  result  of  European 
queens  mating  with  African  drones.  Feral  African  swarms  were  also  tested.  Swarms 
near  the  expanding  front  had  low  levels  of  European  markers,  whereas  those  sampled 
farther  behind  the  front  had  virtually  no  European  markers.  These  findings  pointed  to 
asymmetric  gene  flow:  substantial  African  paternal  introgression  into  managed  European 
colonies  but  limited  paternal  introgression  from  the  European  colonies  into  the  feral 
African  bee  population  (Hall,  1990).  Similar  findings  were  obtained  by  others  with 
allozyme  markers  (Taylor  et  al.,  1991). 

Selection  against  European  bees  and  their  hybrids  in  tropical  conditions  probably 
contributed  to  the  rapid  Africanization  of  the  European  population.  However,  low  hybrid 
survivability  may  also  be  a  result  of  low  metabolic  capacity  compared  to  either  African  or 
European  bees  (Harrison  and  Hall,  1993).  Other  factors,  such  as  queen  developmental 
time  (DeGrandi-Hoffman  and  Watkins,  1998;  DeGrandi-Hoffman  et  al,  1998),  may  have 
also  contributed  to  the  Africanization  process.  A  honey  bee  queen  mates  with  an  average 
of  a  dozen  drones  (Adams  et  al,  1977).  Among  daughter  queens  from  mothers  mated  to 
both  African  and  European  drones,  those  with  African  paternity  will  develop  faster  and 
will  emerge  first,  killing  their  slowly  developing  sisters  still  in  their  cells.  Through  kin 
selection,  workers  with  African  paternity  may  preferentially,  and  more  successfully,  rear 
closely  related  sister  queens.  There  is  yet  no  evidence  to  support  this  possibility. 
Although  European  genes  would  be  expressed  in  the  worker  population,  these  queen 
selective  processes  would  eliminate  them  from  the  reproductive  population. 
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In  the  subtropical-temperate  boundary  region,  where  both  African  and  European 
bees  will  be  equally  adapted,  a  hybrid  zone  is  expected  be  formed  (Taylor  and  Spivak, 
1984).  Bees  with  European  genotypes  will  be  better  adapted  to  the  temperate  climate 
north  of  this  hybrid  zone,  whereas  those  with  African  genotypes  will  be  better  adapted  to 
the  subtropical  conditions  south  of  the  zone. 

Needs  for  Additional  Studies 

More  DNA  markers  that  can  be  easily  used  have  been  needed  to  distinguish 
African  and  European  bees  and  their  hybrids.  In  addition  to  research,  DNA  markers  are 
needed  for  regulatory  identification  of  African  bees  in  the  U.S. A  and  for  stock 
certification. 

Additional  studies  have  been  needed  to  obtain  a  more  accurate  measure  of 
subspecies  contributions,  introgression  patterns,  and  the  resulting  genetic  structure  of 
African  and  European  honey  bee  populations  in  the  New  World.  The  subspecies 
contribution  of  the  feral  population  in  the  neotropics  needs  to  be  better  documented. 
Also  the  extent  that  the  African  European  genotypic  composition  affects  the  fitness  of 
hybrid  colonies  remains  to  be  determined.  Hybrid  fitness  in  either  tropical  or  temperate 
conditions  must  be  established  before  implementing  breeding  programs  proposed  to 
produce  better  adapted  and  gentler  hybrids  for  commercial  beekeeping.  Many  hybrid 
studies  will  depend  upon  genetic  markers  that  are  specific  to  African  and  European  honey 
bee  subspecies,  that  involves  cost  effective  and  accurate  methods  to  screen  honey  bee 
populations. 


Goals  of  This  Research. 


My  doctoral  research  had  two  major  goals.  The  first  goal  was  to  employ  PCR  to 
augment  the  collection  of  markers.  The  second  goal  was  to  apply  the  markers  in  an 
advanced  population  genetic  study  involving  a  broader  sampling  of  New  World  bees.  In 
accomplishing  the  first  major  goal,  nuclear  DNA  PCR-RFLP  markers  were  obtained  with 
three  different  approaches.  First,  markers  found  using  conventional  Southern  blot 
procedures  were  made  PCR  analyzable  by  isolating  and  amplifying  the  regions  that 
contained  most  of  the  polymorphisms  (Chapter  2).  Second,  new  markers  were  sought 
and  found  in  DNA  initially  amplified  using  standard  and  long-PCR  rather  than  in  DNA 
detected  with  cloned  probes  using  Southern  blot  procedures  (Chapters  3  and  4).  The 
effectiveness  of  using  long  versus  standard  PCR  protocols  was  also  compared  (Chapter 
3).  To  reduce  more  stringent  requirements  used  in  long-PCR  protocols,  markers  found 
with  long-PCR  were  made  analyzable  by  standard  PCR  protocols  (Chapter  5).  Finally, 
markers  were  also  obtained  after  modifying  the  amplified  fragment  length  polyorphism 
(AFLP)  procedure,  a  technique  that  can  detect  a  large  amount  of  polymorphism  in  plant 
and  animals  (Chapter  6).  To  accomplish  the  second  main  goal,  the  PCR-RFLP  markers 
found  in  this  study  were  used  together  with  previously  described  PCR-RFLP  markers  to 
test  samples  from  Venezuela,  Costa  Rica,  Honduras,  Mexico,  and  the  U.S. A,  as  well  as 
Old  World  African  samples  of  A.m.scutellata  and  Old  World  European  samples  of 
A.m.mellifera,  A.m.iberica,  A.m.ligustica,  A.m.carnica,  and  A.m.caucasica.  Samples 
from  two  populations  were  included  in  the  analysis,  taken  at  two  different  periods  near 
the  subtropical-temperate  boundary  in  Texas. 


As  an  outcome  of  this  study,  much  more  accurate  estimates  of  gene  frequencies 
and  levels  of  hybridization  were  determined.  New  important  information  was  obtained 
that  either  confirmed  or  required  modifications  of  previous  conclusions. 


CHAPTER  2 

CHARACTERIZATION  OF  POLYMORPHIC  PROBE  P227 

Introduction 

Conventional  methods  used  to  search  for  RFLP  markers  make  use  of  Southern 
blot  procedures  (Southern,  1975)  which  typically  use  cloned  DNA  fragments  are  used  as 
probes.  Restriction  fragment  length  polymorphisms  (RFLPs)  in  mitochondrial  (mtDNA) 
and  nuclear  DNA  (nDNA)  have  been  found  that  distinguish  African  and  European  honey 
bee  subspecies  (Hall  1986;  Hall  1990;  Hall  1992b;  Hall  1998;  Hall  and  Smith  1991). 
These  markers  have  been  valuable  in  revealing  processes  involved  in  the  African  bee 
spread  in  the  New  World  and  in  detecting  accidental  introductions  of  African  bee  swarms 
in  Florida,  USA  (Hall  1990,  Hall  1992a;  Hall  and  Muralidharan  1989;  Sanford,  1999; 
Smith  et  ah,  1989).  Southern  blot  analysis,  however,  is  laborious,  time  consuming, 
expensive  and  requires  considerable  technical  expertise.  Large  amounts  of  DNA  are  also 
required,  thus,  the  number  of  enzymes  that  can  be  screened  to  reveal  polymorphisms  in  a 
particular  sample  is  limited  to  the  amount  of  DNA  that  can  be  extracted.  Therefore,  these 
analyses  are  difficult  to  use  especially  when  a  large  number  of  samples  must  be  tested,  as 
in  population  genetic  studies.  With  the  implementation  of  the  polymerase  chain  reaction 
(PCR)  (Saiki  et  al,  1988),  samples  can  now  be  amplified  in  a  relatively  short  period  of 
time  with  little  amounts  of  DNA  and  at  a  relatively  low  cost.  Large  number  of  samples 
can  be  tested  more  effectively.  Thus,  PCR  analyzable  RFLP  markers  are  more  suitable 
for  population  genetic  studies. 
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In  this  chapter,  a  highly  polymorphic  locus,  P227,  detected  using  conventional 
Southern  blot  procedures,  is  characterized,  and  its  informative  polymorphic  section  made 
analyzable  by  PCR.  This  enabled  the  section  to  be  screened  with  additional  enzymes 
than  those  used  in  the  initial  screening  using  Southern  blots.  Locus  P227  was  found  to 
have  alleles  specific  to  African,  east  and  west  European  honey  bee  subspecies  seen  with 
the  restriction  endonucleases  Alul  and  Hinfi.  Frequency  data  are  reported  for  both 
African  and  European  alleles  detected  with  Alul  and  Hinfl  in  African  and  European 
populations. 

Materials  and  Methods 

Honey  Bee  DNA  Samples  and  DNA  Isolation 

Samples  used  were  drones  (haploid  males)  and  workers  (diploid  females)  from 
four  different  locations  in  the  Transvaal  of  South  Africa:  White  River,  Louis  Trichardt, 
Warm  Baths  and  Pretoria  (collected  by  H.  G.  Hall,  University  of  Florida,  Gainesville, 
FL),  three  locations  from  the  United  States  (U.S.A):  Florida,  Kansas  (provided  by  O.  R. 
Taylor,  University  of  Kansas,  Lawrence,  KS)  and  from  a  closed  breeding  population 
representing  bee  stocks  from  across  the  United  States  in  Tucson,  Arizona  (Page  et  ah, 
1982;  Severson  et  al,  1986)  (provided  by  G.  Loper,  G.  Waller,  J.  Martin  and  E.  Erickson, 
Carl  Hayden  Bee  Research  Center,  USDA-ARS).  Adult  worker  samples  from  France 
(Apis  mellifera  mellifera),  Italy  (A.m.  ligustica),  Spain  (A.m.  iberica)  and  Austria  (A.m. 
carnicd)  (provided  by  D.  R.  Smith,  University  of  Kansas,  Lawrence,  KS)  and  from  a 
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collection  of  bees  collected  in  Russia  and  maintained  in  France  (A.  m.  caucasica)  were 
also  used  to  determine  the  specificity  of  the  European  markers.  One  to  three  workers 
were  used  per  colony  from  a  total  of  67  colonies  from  South  Africa,  44  colonies  from  the 
USA,  nine  colonies  from  France,  41  colonies  from  Italy,  26  colonies  from  Austria,  15 
colonies  from  Spain  and  7  colonies  from  France  (from  a  collection).  European  samples 
were  grouped  as  either  west  European  {A.m.  mellifera  and  A.m.  iberica)  or  east  European 
(A.m.  ligustica,  A.m.  carried).  A.m.  caucasica  was  grouped  as  a  fourth  near-east  lineage 
as  suggested  by  morphometry  and  more  recently,  by  mitochondrial  and  nuclear  DNA 
markers  (Arias  and  Sheppard,  1996;  Ruttner,  1988).  DNA  was  isolated  from  drone  pupae 
and  from  thoracic  muscles  of  adult  honey  bee  workers  as  described  in  Hall,  (1986)  and 
Hall,  (1990),  and  by  using  anionic  columns  as  suggested  by  the  manufacturer  (QIAGEN, 
Palo  Alto,  CA).  Drone  samples  (as  pupae)  were  primarily  used  in  Southern  blot  analysis 
because  of  their  simple  haploid  pattern  of  restriction  fragments. 

Probe  Description  and  Screening  for  Polymorphisms. 

Construction  of  a  plasmid  library  of  cloned  honey  bee  fragments,  electrophoresis 
of  digested  DNA,  blotting,  probe  labeling,  prehybridization,  hybridization  and 
autoradiography  were  as  described  in  Hall  1986  and  Hall  1995.  A  honey  bee  genomic 
DNA  library  was  generated  from  Pstl  digested  honey  bee  DNA  cloned  into  Pstl  digested 
pBR322  as  plasmid  vector.  The  bee  DNA  fragments  were  from  low-copy  number 
sequences  and  were  larger  than  4  kilobases  (kb)  (Hall,  1986).  Probe  227  (P227)  selected 


20 

from  this  library  had  a  9.2  kilobase  (kb)  insert.  To  screen  for  enzymes  that  reveal 
polymorphisms,  aliquots  of  5ug  of  DNA  from  one  African  and  from  one  European  drone 
sample  was  separately  digested  with  nine  different  restriction  endonucleases:  Mspl, 
Sau96l,  AM,  HaeUl,  Hhal,  Neil,  Mbol,  Hinjl  and  Ddel.  The  DNA  digested  with  the 
different  enzymes  was  electrophoresed  in  separate  wells  of  a  2%  agarose  gel  and 
transferred  to  a  nitrocellulose  membrane  (blots)  (Genescreen  plus,  NEN  Research 
products,  Boston,  MA).  Probe  was  labeled  by  Nick  translation  as  described  (Sambrook  et 
ah,  1989).  Blots  were  prehybridized,  hybridized  with  radioactively  labeled  P227  as 
probe  and  exposed  to  autoradiographic  film  (Kodak  X-Omat,  Rochester,  NY).  From  this 
initial  blot,  enzymes  that  revealed  distinctive  polymorphisms  between  the  African  and 
European  honey  bee  DNA  samples  were  selected  for  further  testing.  Additional  blots 
(three  to  four)  containing  DNA  from  nine  African  and  nine  European  samples  digested 
with  the  enzymes  that  revealed  polymorphisms  during  the  initial  screening  were  used  to 
determine  allele  frequencies. 

Subcloning 

In  this  chapter  and  in  Chapter  5,  the  term  "region"  will  refer  to  the  entire  locus 
under  study.  Alleles  will  be  the  different  banding  patterns  found  at  a  particular  locus  or 
"region".  The  conversion  of  markers  found  with  Southern  blot  procedures  and  long  PCR 
required  the  cloning  of  "sections"  of  each  region.  Alleles  resulting  from  RFLP  analysis 
of  these  sections  will  be  referred  to  as  "sub-alleles".  Thus,  the  terms  "regions"  and 
"alleles"  are  analogous  to  "sections"  and  "sub-alleles". 
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Restriction  sites  of  the  enzymes  contained  in  the  polycloning  site  of  the  plasmid 
vector  pGem3Z  (Promega  corporation,  Madison,  WI)  were  mapped  in  P227.  To 
subdivide  the  entire  probe  into  several  sections  of  sizes  amplifiable  by  standard  PCR 
protocols,  three  restriction  endonucleases  were  selected,  based  on  the  size  of  fragments 
they  generated  and  their  relative  position  in  P227:  Xhol,  Sphl  and  HindlR.  The  relative 
positions  of  the  selected  enzymes  were  mapped  in  P227  using  standard  double  digestion 
procedures  (Sambrook  et  al.,  1989).  These  enzymes  divided  P227  in  four  different 
sections:  section  1  (SI),  from  the  terminal  Pstl  site  (position  0)  to  the  Xhol  site  (position 
2500)(2.5kb),  section  2  (S2)  from  the Xhol  site  to  the  Sphl  site  (position  3900)(1.4kb), 
section  3  (S3)  from  the  Sphl  site  to  the  Hincffil  site  (position  7100)(3.2kb)  and  section  4 
(S4)  from  the  HincM  site  to  the  second  terminal  Pstl  site  (position  9200)(2.1kb)  (Figure 
2-1).  Each  section  was  subcloned  into  pGem3Z  vector  (Promega  corporation,  Madison, 
WI)  and  transformed  into  E.  coli  competent  cell  DH5a  (Lifetechnologies,  GIBCO  BRL, 
Gaithersburg,  MD).  To  clone  the  sections,  500ng  of  P227  and  500ng  of  pGem3Z 
(Promega  corporation,  Madison,  WI)  were  digested  with  the  enzymes  that  flanked  the 
corresponding  section.  The  digestion  was  done  in  50ul  of  total  volume  with  10  units  of 
each  enzyme  and  using  the  buffers  and  conditions  following  manufacturers 
recommendations.  Fragments  for  both  the  insert  (section  of  P227  to  subclone)  and  the 
vector  (pGem3Z)  were  separated  in  1%  agarose  gels,  excised  with  a  blade,  and  placed  in 
dialysis  bags.  Bags  were  filled  with  IX  concentrated  Tris  Borate  Buffer  (TBE).  To  elute 
the  DNA  fragments  from  the  agarose  gel  slice,  the  dialysis  bag  was  placed  in  an 
electrophoresis  unit  and  electrophoresed  at  4V/cm  for  2hrs.  The  buffer  containing  the 
eluted  DNA  fragments  was  collected,  and  DNA  fragments  were  purified  using  anionic 
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resin  columns  (DE52,  Whatman,  Maidstone,  England)  as  described  (Sambrook  et  al, 
1989).  DNA  fragments  were  precipitated  with  500ul  of  95%  ethanol,  washed  with  500ul 
of  75%  ethanol  and  resuspended  in  20ul  of  water.  Equal  amounts  of  both  purified  vector 
and  insert  (lng  of  each)  were  ligated  in  a  total  of  20ul  with  limit  of  T4  ligase 
(Boehringer  Mannheim,  Indianapolis,  IN)  and  buffer  conditions  as  suggested  by  the 
manufacturers.  Ligation  reactions  were  performed  at  1 5°C  for  1 6hrs.  Colony 
transformation  was  done  with  lOOul  of  competent  cells  DH5a  (Lifetechnologies,  Gibco 
BRL,  Gaithersburg,  MD).  To  transform  the  cells,  4ul  of  the  ligation  mixture  were  place 
on  top  of  the  competent  cells  and  mixed  gently.  Cells  were  incubated  on  ice  for  30  min, 
heat  shocked  at  42°C  for  45  sec  and  placed  on  ice  for  2  min.  Following  incubation  on 
ice,  900ul  of  SOC  medium  (Sambrook  et  al,  1989)  were  added,  and  the  mix  was 
incubated  at  37°C  for  lhr.  The  cells  were  plated  on  SOB  medium  (Sambrook  et  al, 
1989)  and  allowed  to  grow  overnight.  Cells  were  then  screened  for  white  or  blue 
colonies  for  the  presence  or  absence  of  inserts,  respectively.  Subclones  corresponding  to 
each  section  were  used  separately  as  radioactive  probes  in  the  original  blots,  to  find  the 
sections  that  contained  the  informative  polymorphic  sites.  Subclones  revealing 
informative  polymorphisms  were  sequenced  about  500  bases  from  their  terminal  ends 
(University  of  Florida,  Interdisciplinary  center  for  biotechnology  research  ICBR, 
Gainesville,  FL).  From  these  terminal  sequences,  primer  sequences  were  selected  using 
the  computer  software  package  OLIGO  5.0  (National  Bioscience,  Plymouth,  MN). 
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Amplification  and  RFLP  Analysis  of  PCR  Products 

Almost  the  entire  S3  of  P227  with  the  adjacent  end  of  S4  was  amplified  as  one 
segment,  referred  to  as  Section  3  extended  (P227-S3xt).  Primers  used  were  [1]  5'- 
AGAAGGAAAGAAGAAACGGATGAAC-3 '  within  S3  (position  4030)  near  the  Sphl 
site  and,  [2]  5  '-CGGAGGAGTGGTAATAATGGAAGC-3 '  within  S4  (position  7496) 
near  the  HinaW  site  (Figure  2-1).  Amplifications  were  performed  in  25  ul  of  50mM  Tris- 
HC1  pH  9.2, 16mM  (NH4)  2S04  and  1 .8mM  MgCl2  containing  200uM  of  each  dNTP's, 
250nM  of  each  primer  with  1.12  units  of  Taq  DNA  polymerase  and  0.12  units  of  Pwol 
DNA  polymerase  taken  from  a  stock  mixture  containing  nine  units  of  Taq  DNA 
polymerase  and  one  unit  of  Pwol  DNA  polymerase  (Boheringer  Mannheim,  Indianapolis, 
IN).  Fifty  to  250ng  of  genomic  DNA  was  used  per  reaction  (the  highest  amount  of  DNA 
was  used  when  DNA  was  extracted  from  thoracic  muscles  of  adult  honey  bees). 
Amplifications  were  performed  in  a  PTC- 100  or  PTC-200  thermocycler  (MJ  Research, 
Watertown,  MA)  and  the  following  profile  was  used:  95°C  for  2  min  followed  by  33 
cycles  each  at  94°C  for  45s,  57°C  for  45s  and  68°C  for  4min  with  a  final  extension  at 
68°C  for  lOmin.  A  second  section  near  the  Sphl  terminus  of  S3  (position  4030  to  4300) 
referred  to  as  P227  section  3  terminal  (P227-S3ter)  was  amplified  using  primer  [1]  from 
P227-S3xt  and  primer  [3]  5  '-CACTTTTTCAAAAGAGCGTGCAA-3 '  (position  4300) 
(Fig.  2-1).  The  reaction  solution  to  amplify  P227-S3ter  was  as  with  P227-S3xt,  except 
that  l.OmM  MgCl2  and  one  unit  of  Taq  DNA  polymerase  were  used  (Lifetechnologies 
GIBCO  BRL,  Gaithersburg,  MD).  The  amplification  profile  for  P227-S3ter  was  as 
follows:  95°C  for  2  min  followed  by  35  cycles  of  94°C  for  45sec,  60°C  for  45sec  and 
72°C  for  30sec  with  a  final  extension  of  5  min  at  72°C.  P227-S3xt  was  screened  for  22 
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restriction  endonucleases,  some  of  which  had  been  used  in  the  initial  screening  using 
southern  blots:  HaeUl,  Hinfl,  Mspl,  Dral,  Avail,  BamHl,  £coRI,  Ddel,  Ajlti.,  BgR,  NrUl, 
Stul,  Bell,  Nhel,  BstXl,  EcoRV,  Spel,  Hpal,  Nsil,  flgffl,  Ncol  and  Styl.  Five  ul  of 
amplified  products  were  digested  with  5ul  of  a  digestion  mixture  containing  1  unit  of 
enzyme,  2mM  DTT,  25ug  of  BSA,  and  restriction  digestion  buffer  to  a  final 
concentration  of  IX  as  suggested  by  the  manufacturers.  Digested  PCR  products  were 
separated  in  2%  agarose  gels,  stained  with  ethidium  bromide,  and  visualized  over  UV 
light.  The  initial  screening  was  done  with  nine  African  and  nine  European  samples  to 
increase  the  chances  of  finding  polymorphisms,  compared  to  one  African  and  one 
European  drone  used  in  the  initial  screening  with  Southern  blots.  Enzymes  that  revealed 
polymorphisms  in  all  of  one  group  and  not  the  other  were  tested  in  a  larger  sample  size 
consisting  of  48  African  and  48  European  drone  samples.  The  restriction  enzymes  Alul 
and  Hinfl  were  found  to  generate  informative  polymorphisms. 

PCR  products  were  used  to  map  the  restriction  sites  for  all  the  alleles  of  the 
enzymes  that  revealed  significant  polymorphisms  following  the  method  described  by  Her 
and  Weinshilboum,  1995  and  described  in  more  detail  in  Chapter  5. 

Sub-Allele  Count  and  Data  Analysis 

The  different  fragment  patterns  generated  by  Alul  and  Hinfl  were  identified  in 
both  drones  and  workers.  Each  pattern  identified  is  considered  a  sub-allele.  Sub-alleles 
were  named  using  the  enzyme  name  followed  by  a  letter,  for  example,  Alul- A.  Drones 
(haploid  parthenogenetic  progeny)  were  used  to  determine  the  queen  genotypes  from 
which  sub-allele  frequencies  were  obtained.  Two  to  five  drones  per  colony  were  used. 
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When  three  or  more  drones  where  found  to  have  the  same  sub-allele,  the  queen's 
genotype  was  determined  as  homozygous  for  that  sub-allele,  and  the  sub-allele  was 
counted  twice.  If  two  different  sub-alleles  were  found,  each  one  was  scored  once.  Each 
drone  was  assigned  a  combination  of  two  letters  corresponding  to  the  sub-alleles 
produced  by  Alul  and  Hinfi.  These  combinations  are  referred  to  as  composite  sub-alleles. 
The  coupled  or  repulsion  phase  of  sub-alleles  produced  by  two  restriction  enzymes,  thus 
the  identity  of  the  composite  sub-alleles,  could  not  be  determined  in  diploid  individuals 
i.e.  workers.  For  example,  a  worker  with  genotype  Alul-A:Alul-BI Hinfi- A:Hinfl-B  for 
both  Alul  (heterozygous  for  the  Alul)  and  Hinfi  (heterozygous  for  Hinfi)  will  have  four 
possible  composite  sub-alleles:  Alul-A/Hinfl-A,  Alul-A/Hinfi-B,  Alul-B/Hinfi-A  and^M- 
B/Hinfi-B.  In  this  case,  a  cis  configuration  (Alul-A:Hinfi-A)  cannot  be  distinguished  from 
a  trans  configuration  (Alul-A:  Hinfi-B).  Frequency  data  obtained  from  composite  sub- 
alleles  in  drone  samples  was  used  to  generate  a  contingency  table.  Sub-allele  association 
was  tested  by  testing  the  contingency  table  for  independence  using  a  Fisher's  exact  test 
(under  the  null  hypothesis  of  no  association)  based  on  the  metropolis  algorithm  using  the 
computer  software  package  RXC  (M.Miller,  University  of  Arizona).  Frequency  data 
from  workers  was  obtained  from  direct  counts  of  Alul  and  Hinfi  sub-alleles.  Observed 
and  Expected  heterozygosities  (H^s  and  HeXp)  were  calculated  as  described  in  Nei,  1987 
and  Nei,  1978  using  the  computer  software  package  GENEPOP  version  1.2  (Raymond 
andRousset,  1995). 
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Results 

Initial  screening,  using  P227  as  probe  in  Southern  blots  with  DNA  from  one 
African  and  one  European  sample  digested  with  9  enzymes,  revealed  informative 
polymorphisms  with  AM.  Further  analysis  using  blots  containing  nine  African  and  nine 
European  samples  showed  a  total  of  nine  alleles  of  which  six  were  African-specific  and 
three  were  European-specific  (Figure  2-2).  All  polymorphic  bands  resulting  from  AM 
digests  were  detected  when  subclones  of  both  section  3  (S3)  and  section  4  (S4)  of  P227 
were  used  as  probes.  S4  detected  only  one  AM  polymorphic  site  that  was  located  at  the 
terminal  end,  adjacent  to  S3  (determined  from  sequencing  data).  RFLP  analysis  of  the 
amplifiable  P227-S3xt  also  revealed  polymorphisms  with  HirtfL 

AM  Polymorphisms 

AM  digestions  of  P227-S3xt  region  revealed  nine  sub-alleles  (Figure  2-3).  Sub- 
alleles  AM-A,  C,  D,  and  J  were  found  predominantly  in  the  South  African  population  at  a 
collective  frequency  of  52.4%  (Table  2-1).  AM-J,  found  at  a  39.5%  frequency  in  the 
African  samples,  was  absent  from  both  east  and  west  European  types  but  present  at  a  very 
low  frequency  (0.8%)  in  the  U.S.A  samples.  Two  sub-alleles  (AM-I  and  E)  were 
predominantly  found  in  the  European  samples.  AM-I  was  found  in  east  and  west 
European  samples  at  frequencies  of  100.0%  and  79.1%  in  A.m.ligustica  and  A.m.carnica, 
respectively,  and  65.0%  and  30.7%  in  A.m.mellifera  and  A.m. iberica,  respectively.  AM-I 
was  absent  in  bees  of  the  near-east  lineage  group  of  A.m.caucasica.  AM-E  was 
A.m.caucasica-specific  and  was  found  at  a  frequency  of  92.1%  among  the  A.m.caucasica 
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samples.  Sub-alleles  Alul-I and  E  were  found  at  a  collective  frequency  of  83.5%  in  the 
U.S.A  samples. 

Sub-alleles  Alul-Al,  D  and  /have  unique  fragments  of  1000  (Al),  2834  (D)  and 
1 650  bp  (I),  respectively.  Sub-alleles  Alul-B  and  E,  Alul-E  and  F  and,  Alul-F  and  J  are 
characterized  by  fragments  of  2574,  340  and  2354bp,  respectively.  An  AM  site,  located 
220  bases  from  primer  [1]  within  S3  (close  to  the  Sphl  site)  is  present  in  sub-alleles  Alul- 
A,  Al,  C,  D,  F  and  J,  most  of  which  are  African-predominant  or  African-specific  (except 
Alul-Al  and  F),  and  absent  in  Alul-E,  /  (European-predominant  sub-alleles)  and  B 
(common  sub-allele).  Sub-alleles  Alul-J  (African-specific)  and  B  (common)  have  very 
similar  banding  patterns  but  differ  by  the  presence  and  absence,  respectively,  of  this 
restriction  site,  at  position  220,  which  generates  the  characteristic  fragment  of  2574bp  in 
Alul-B  and  2354bp  in  Alul-J.  The  220bp  fragment  in  Alul-J comigrates  with  a  separate 
225bp  fragment,  and  both  appear  in  agarose  gels  as  only  one  band,  making  the  distinction 
between  Alul-B  and  J  difficult.  Distinction  between  these  two  sub-alleles  can  be  made  by 
amplifying  and  digesting  a  shorter  section  (P227-S3ter).  The  presence  of  a  280bp 
fragment  (uncut)  identifies  the  Alul-B  sub-allele  (common  sub-allele)  whereas  its  absence 
identifies  the  Alul- J  sub-allele  (African-specific)  (Figure  2-3). 

Genotypes  that  cannot  be  distinguished  include:  A/A  from  AIC,  A1IA1  from  A 1IC 
and  All  from  CI  I.  The  inability  to  distinguish  these  genotypes  can  create  a  bias  in  the 
frequency  of  both  Alul- A  and  Alul-Al,  however,  Alul-C  was  found  at  very  low 
frequencies  in  our  samples  (only  in  one  drone,  corresponding  to  a  frequency  of  0.3%). 
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Hinfi  Polymorphisms 

Five  sub-alleles  were  revealed  with  Hinfi,  of  which  three  were  European-specific 
(Hinfi-F,  G  and  H).  Hinfi-F  was  only  found  in  east  European  samples  (56.1%  in  Am. 
lignstica  and  66.7%  in  A.m.  carnica)  and  in  35.9%  of  the  A.m.  caucasica  samples.  Hinfi- 
G  was  fixed  in  the  west  European  samples  (100.0%  in  both  A.m.mellifera  and  A.m.iberica) 
and  was  found  at  an  18.8%  frequency  in  the  A.m.  carnica  and  at  a  10.3%  frequency  in  the 
A.m  caucasica  samples.  Hinfi-F,  G  and  //were  found  at  a  collective  frequency  of  74.2% 
in  the  U.S.A  samples  (Table  2-1).  Hinfi-H  was  found  only  in  the  U.S.A  samples  at  an 
1 1.7%  frequency.  One  African-specific  sub-allele  was  found  (Hinfi-C)  at  a  frequency  of 
less  than  1 .0%  in  the  South  African  samples. 

Sub-alleles  Hinfi-C,  F,  G  and  H  were  characterized  by  fragments  of  400  (Q,  380 
(F),  250  and  225  (G),  and  470bp  (H)  respectively  (Figure  2-4).  Undistinguishable 
genotypes  include  CIC  from  CID  and  F/F  from  FID. 

A  strong  association  between  Alul  and  Hinfi  sub-alleles  was  found  (P<0.0001  ± 
0.000).  These  sub-allele  associations  reinforce  the  specificity  and  hence,  the  value  of 
these  markers  for  diagnostic  purposes.  For  example,  Alul-I,  a  predominantly  east 
European  sub-allele  was  generally  found  associated  with  Hinfi-F,  an  east  European- 
specific  sub-allele  (Table  2-2).  African-predominant  or  specific  sub-alleles  were  found 
associated  with  African  or  common  sub-alleles.  Hinfi-H  may  not  be  a  European  sub- 
allele.  It  was  not  found  in  samples  from  Europe,  and  it  was  associated  with  an  African- 
specific  sub-allele  (Alul- J). 
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Observed  heterozygosity  values  for  sub-alleles  found  with  AM  were  higher  in 
South  African  and  west  European  bees  compared  to  east  European  and  near-east  bees. 
Compared  to  the  AM  sub-alleles,  heterozygosity  values  found  with  HinfL  sub-alleles  was 
higher  in  east  European  bees  compared  to  South  African  and  west  European  (Table  2-1). 

Discussion 

We  have  employed  the  polymerase  chain  reaction  to  amplify  a  polymorphic 
region,  P227,  initially  detected  by  Southern  blot  analysis,  with  sub-alleles  specific  to  or 
predominant  in  both  African  and  European  honey  bee  subspecies.  More  importantly, 
P227  effectively  distinguishes  bees  in  the  African  {A.m.  scutellata),  west  European  {A.m. 
mellifera  and  A.m.  iberica)  and  east  European  {A.m.  ligustica,  A.m.  carried),  and  near- 
east  {A.m.  caucasica)  group  of  subspecies.  Sub-allele  frequency  data  obtained  from  AM 
and  Hinfl.  digests  of  P227-S3xt  reflects  the  genetic  composition  of  bees  from  the  United 
States  and  the  contributions  of  different  subspecies  introduced  by  European  settlers  to 
North  America  (Nelson,  1967;  Oertel,  1976;  Pellet,  1938).  The  combined  European  sub- 
allele  frequency  in  the  North  American  population  is  consistent  with  that  found  for  alleles 
at  other  polymorphic  loci  (Hall,  1990;  McMichael  and  Hall,  1996).  Most  of  these  sub- 
alleles  correspond  to  east  European-specific  sub-alleles,  which  reflects  the  preference  of 
the  Italian  honey  bee,  A.m.lignstica,  among  US  beekeepers.  West  European  bees  were 
introduced  to  North  America  prior  to  the  introduction  of  east  European  bees,  but,  because 
of  their  defensive  behavior,  they  were  not  as  popular  among  beekeepers  (Oertel,  1976; 
Pellet,  1938).  This  is  shown  by  the  lower  frequency  of  the  west  European-specific  sub- 
alleles  found  with  Hinfl.  {Hinfl-G).  Low  frequencies  of  the  A.m.caucasica  sub-allele 
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Alul-E  in  the  U.S.A  samples  also  show  a  genetic  contribution  of  this  subspecies  in  the 
U.S.A. 

Honey  bees  from  Egypt  {A.m.  lamarckii),  and  possibly  from  other  African 
subspecies,  were  introduced  in  the  U.S.A  in  the  early  18th  century  (Oertel,  1976;  Pellet, 
1938).  Low  adaptability  to  temperate  conditions  probably  resulted  in  failure  to  establish 
a  substantial  managed  or  feral  population,  as  did  the  west  and  east  European  subspecies. 
However,  it  is  possible  that  alleles  from  these  African  subspecies  could  have  introgressed 
and  persisted  at  low  frequencies  in  feral  and  managed  colonies  in  North  America.  Sub- 
alleles  found  at  very  low  frequencies  in  South  Africa  and  the  U.S.A  (example:  Alul-F) 
could  possibly  reflect  the  contribution  of  these  ncm-A.m.scutellata  African  subspecies  in 
the  New  World.  Introduction  of  non-A.m.scutellata  African  subspecies  to  America  has 
been  suggested  based  on  mitochondrial  DNA  evidence  in  South  American  populations 
(Arias  and  Sheppard,  1996).  Migratory  beekeeping  practices  could  also  contribute  to  the 
presence  of  low  frequency  sub-alleles  in  other  populations  as,  for  example,  the  regular 
movement  of  A.m.  ligustica  between  Italy  and  France  where  A.m.  mellifera  is 
predominantly  found  (Estoup  et  al,  1995).  African  subspecies  such  as  A.m.  intermissa 
brought  from  Morroco  (Northwest  Africa)  into  the  Iberian  peninsula  where  A.m.  iberica 
is  predominantly  found  (Garnery  et  al,  1978a  and  b;  Franck  et  al,  1998).  European- 
predominant  sub-alleles  are  also  found  at  low  frequencies  (less  than  4%)  in  the  South 
African  population.  European  bees  were  imported  to  South  Africa  but  did  not  establish  a 
feral  or  managed  population  because  of  their  poor  adaptability  to  tropical  conditions 
(Fletcher,  1973;  Fletcher,  1978).  From  allozyme  (Smith  and  Glenn,  1995),  mtDNA 
(Smith  et  al,  1991),  and  microsatellite  (Franck  et  al.  1998)  studies,  African-specific 
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alleles  have  been  found  in  west  European  populations,  particularly  A.m.  iberica.  Recent 
studies  have  proposed  that  the  Iberian  peninsula  could  be  an  intergradation  zone  in  which 
west  European-specific  and  African  alleles  are  found,  the  latter  having  introgressed  from 
North  Africa  (Smith  et  al.  1991  and  Franck  et  al.  1998).  Alul-Al  is  found  at  a  high 
frequency  in  A.m.  iberica,  at  intermediate  frequencies  in  A.m.scutellat,a  and  at  very  low 
frequencies  in  the  U.S  A  samples.  Alul-Al  may  be  a  sub-allele  that  originated  from  a 
north-western  African  subspecies  and  introgressed  north  into  the  Iberian  peninsula  and 
south  into  the  A.m.  scutellata  populations.  This  hypothesis  cannot  be  substantiated  with 
our  current  data. 

African  and  west  European-specific  sub-alleles  have  been  more  difficult  to  find 
than  east  European-specific  sub-alleles.  High  swarming,  high  genetic  diversity  and 
migratory  beekeeping  practices  may  have  allowed  for  a  greater  blending  of  more  diverse 
populations,  and  this  may  have  contributed  to  the  difficulties  of  finding  African-specific 
markers.  Markers  found  with  P227  include  not  only  African-specific  but  also  west  and 
east  European-specific  alleles,  which  increases  the  value  of  this  locus  thereby  facilitating 
the  use  of  these  markers  in  studying  gene  flow  or  genetic  variation  between  African  and 
European  populations  in  the  New  World. 

This  polymorphic  region,  P227,  was  successfully  detected  using  Southern  blot 
procedures,  and  the  sections  containing  the  polymorphic  sites  were  made  analyzable  by 
PCR.  PCR  facilitates  RFLP  analysis  of  the  polymorphic  sections  by  making  the 
procedure  faster,  easier  to  use,  less  expensive  and  by  reducing  the  quality  and  quantity  of 
DNA  required.  PCR  analyzable  RFLPs  of  P227  are  codominant  markers,  (homozygous 
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and  heterozygous  individuals  can  be  distinguished),  thus,  allowing  these  markers  to  be 
used  for  parental  analysis  and  hybridization  studies. 
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Table  2-2.  Contingency  table  of  composite  sub-allele  frequencies  for  Alul  (columns)  and 
Hinfi  (rows)  sub-alleles.  Sub-allele  association  frequencies  were  determined  in  drone 
samples. 


A*        C*        £>¥        /        E*         fi         F        Al  B 

C*  2 

F*  28 
G§§  4 

H**  2  2 

D  5  1  1  9  5  3  2         15  5 

¥  :  African-specific  or  predominant  sub-alleles. 

§:  East  European-specific  or  predominant  sub-allele. 

§§:West  European-specific  or  predominant  sub-allele. 

*:  A.m.caucasica  -specific  sub-allele. 

**:  Found  only  in  the  U.S.A. 

PO.0001  ±  0.0000  (Fisher's  exact  test). 


Figure  2-1.  Restriction  enzyme  cleavage  site  map  for  locus  P227  and  the  amplifiable 
polymorphic  sections.  I)  General  map  of  locus  227  with  the  enzymes  used  to  divide  and 
subclone  the  different  sections.  II)  Position  of  AM  and  Hinjl  sites  in  the  polymorphic 
section  3  extended  (P227-S3xt).  HI)  Amplified  sections  of  P227-S3ter  used  to 
distinguish  sub-alleles  Alul-B  and  Alul-J.  Locus  P227  was  divided  in  four  sections,  SI  to 
S4  defined  by  the  restriction  sites  Xhol,  Sphl,  and  HincM.  Squares  and  circles 
correspond  to  AM  and  Hinfi  sites,  respectively.  Open  and  closed  symbols  represent 
monomorphic  and  polymorphic  sites,  respectively.  Letters  in  parentheses  are  the  sub- 
alleles  that  have  the  indicated  AM  and  Hinfi  sites.  Position  of  restriction  sites,  in  base 
pairs,  are  indicated  to  the  right  (Alul)  or  left  (Hinfi)  side  of  the  map.  Primers  used  to 
amplify  each  region  are  indicated  in  brackets  above  and  below  each  diagram.  A  single 
individual  star  indicates  the  polymorphic  Hinfi  site  that  is  west  European-predominant 
{A.m.mellifera  and  A.m.iberica).  The  amplified  sections  of  P227-S3ter  (far  right) 
includes  the  polymorphic  site  that  distinguishes  sub-alleles  AM-B  (common)  andAM-J 
(African-specific). 
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Figure  2-2.  Autoradiograph  of  alleles  revealed  by  Alul  digestions  of  probe  P227  in 
African  and  European  drone  samples  as  revealed  by  Southern  blot  analysis.  Numbers  on 
top  of  each  lane  are  the  individual  alleles.  Fragment  sizes  are  indicated  in  kilobase  pairs. 
Alleles  1,  2,  3,  4,  5,  and  6  are  African-specific  and  alleles  7,  8  and  9  are  European- 
specific  (numbers  were  used  in  this  case  to  distinguish  alleles  found  in  Southern  blots  to 
those  found  with  PCR-RFLP  analyses.  These  alleles  do  not  necessarily  correspond  to  the 
PCR-RFLP  sub-alleles). 


Figure  2-3.  Restriction  digest  profiles  of  P227-S3xt  (I),  and  P227-S3ter  (II)  with  Alul. 
AF=African  sub-alleles,  COcommon  sub-alleles  and  EU=European  sub-alleles.  M50, 
Ml 00,  and  M250  correspond  to  molecular  size  standard,  50,  100,  and  250bp  ladders, 
respectively  (Lifetechnologies,  GIBCO  BRL,  Gaithersburg,  MD  ).  Fragment  sizes,  in 
base  pairs,  are  shown  to  the  right.  DNA  samples  amplified  were  all  from  drones. 
African-predominant  sub-alleles  are  characterized  by  fragments  of  1430bp  (Alul-A,  and 
Q  and  2834bp  (Alul-D).  Alul-J,  an  African-specific  sub-allele,  does  not  have  a  unique 
fragment  and  is  distinguished  from  the  others  by  its  unique  banding  pattern.  The 
European-predominant  sub-allele  Alul-1  is  distinguished  from  the  others  by  a  1650bp 
fragment  and  Alul-E  by  its  unique  banding  pattern.  Alul-E  is  an  A.m.  caucasica-specific 
sub-allele.  Sub-alleles  Alul-B  and  Alul-J  are  difficult  to  distinguish  in  RFLP  patterns  of 
P227-S3xt  (I)  but  are  clearly  distinguishable  in  Alul  digests  of  P227-S3ter  (H). 
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Figure  2-4.  Restriction  digestion  profile  for  P227-S3xt  with  HinJL  AF,  CO,  EU  and 
Ml 00  are  as  in  Figure  2-3.  European-specific  sub-alleles  are  characterized  by  fragments 
of  380  (Hinfi-F),  250  and  225  (Hinfi-G),  and  470  (Hinfi-H)bp.  Allele  Hinfi-F  is  east 
European-specific,  Hinfl-G  is  west  European-specific.  Hinfi-H  is  found  only  in  the 
U.S. A  samples. 


CHAPTER  3 

DETECTION  OF  NUCLEAR  RFLPs  IN  AFRICAN  AND  EUROPEAN  HONEY  BEE 
DNA  AMPLIFIED  WITH  STANDARD  AND  LONG  PCR. 

Introduction 

In  Chapter  2,  a  locus  was  described  with  DNA  restriction  fragment  length 
polymorphisms  (RFLPs)  that  distinguish  African  and  European  honeybees.  The  RFLPs 
were  sought  and  found  using  cloned  probes  in  Southern  blot  analyses.  Then,  to  facilitate 
DNA  tests,  the  region  with  the  RFLP  alleles  was  made  analyzable  using  the  polymerase 
chain  reaction  (PCR).  Other  loci  with  RFLPs  distinguishing  African  and  European 
honeybees  and  found  in  Southern  blots  have  also  been  made  PCR  analyzable  (Hall, 
1998).  This  conversion  is  time  consuming:  the  position(s)  of  the  informative 
polymorphic  site(s)  in  the  clone  must  be  ascertained,  so  that  the  polymorphic  region  can 
be  subcloned  and  subsequently  amplified.  In  another  approach,  regions  of  DNA  are  first 
amplified,  using  primer  sequences  obtained  from  anonymous  clones,  and  the  regions  are 
then  screened  for  RFLPs  using  a  large  collection  of  enzymes  (Karl  et  al,  1992;  Karl  and 
Avise,  1993).  The  short  length  of  DNA  amplified  by  standard  PCR  limits  the  amount  of 
variation  observed,  and,  thus,  a  large  number  of  clones  must  be  screened  to  find  useful 
polymorphisms.  With  modifications  of  the  PCR,  regions  of  DNA  of  up  to  35  kilobase 
pairs  (Kbp)  can  now  be  amplified  (Barnes,  1994;  Cheng  et  al,  1994a;  Cheng  et  al, 
1994b).  In  this  chapter,  we  describe  markers  that  distinguish  African  and  European  bees 
found  by  both  standard  and  long  PCR  approaches,  thereby  demonstrating  their  relative 
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effectiveness.  As  expected,  screening  for  RFLPs  in  long  amplified  regions  was  more 
productive  than  in  short  regions.  Associations  between  alleles  generated  by  pair  of 
enzymes  at  a  locus  were  also  determined  which  reinforce  the  African  and  European 
specificity  of  the  markers  found. 

Materials  and  Methods 

Honey  Bee  DNA  Samples 

Samples  were  drone  larvae  and  pupae  from  South  Africa  and  the  United  States 
(European)  as  described  previously  in  Chapter  2.  Drone  samples  were  primarily  used 
because  of  their  simple  haploid  pattern  of  restriction  fragments.  Other  samples  came 
from  Italy,  ten  workers  (diploid  females)  from  five  different  colonies  of  the  east 
European  subpecies  A.  m.  ligustica,  and  from  France,  ten  workers  from  five  colonies  of 
the  west  European  subspecies  A.  m.  mellifera  (provided  by  B.  Vaissiere,  and  J-M. 
Cornuet,  INRA,  France).  The  samples  of  European  workers  were  larvae  or  pupae, 
whereas  in  the  work  reported  in  Chapter  2,  the  worker  samples  were  adults.  Better 
quality  and  quantities  of  DNA  can  be  extracted  from  larvae  or  pupae  (suitable  for  long 
PCR  analysis).  The  worker  larvae  were  used  to  determine  the  origin  of  European- 
specific  alleles  but  not  to  determine  allele  frequencies. 

DNA  Isolation 

DNA  isolation  was  as  described  (Hall,  1990)  with  modifications  to  increase 
quality  and  decrease  shearing  of  the  DNA,  enabling  amplification  of  long  fragments. 
Two  phenol-chloroform  extractions  were  used  instead  of  one,  and,  in  each  extraction,  the 
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samples  were  gently  mixed  for  30min.  Wide  bore  pipette  tips  were  used  to  avoid 
mechanical  damage  of  the  DNA.  DNA  isolations  were  also  performed  using  anionic 
columns  (QIAGEN,  Valencia,  CA)  according  to  the  manufacturer's  recommendations. 
DNA  concentration  was  determined  by  UV  absorbance  at  260nm,  and  samples  were 
diluted  in  TE  (lOmM  Tris-HCl  (pH  8.0),  ImM  EDTA)  to  a  final  concentration  of  lOOng/ 

Ml. 

Clone  Selection. 

Random  pieces  of  honeybee  DNA  digested  with  Xbal  and  Pstl  were  cloned  in  the 
multicloning  site  of  pGEM3Z  (Promega,  Madison,  WI).  Plasmids  were  transformed  into 
DH5a  competent  cells  (Lifetechnologies,  Gibco  BRL,  Grand  Island,  NY)  and  screened 
on  bacterial  plates  with  X-gal  and  ampicillin  for  white  colonies  with  inserts.  Plasmid 
minipreps  were  obtained  by  the  alkaline  lysis  method  as  described  in  Sambrook  et  al, 
1 989.  Insert  sizes  were  verified  by  digesting  the  DNA  minipreps  with  Xbal  and  Pstl  and 
separating  the  fragments  by  agarose  gel  electrophoresis.  Two  libraries  were  established: 
clones  with  inserts  in  the  range  of  1-2  Kbp,  for  amplification  with  standard  PCR;  and 
clones  with  inserts  of  more  than  5  Kbp,  for  amplification  with  long  PCR.  Clones  were 
designated  as  XPS  or  XPL  (Xbal  and  Pstl,  short  or  long  inserts)  followed  by  a  number. 
Five  clones  with  short  inserts  and  five  clones  with  long  inserts  were  randomly  selected. 
About  400  bps  from  both  ends  of  each  insert  were  sequenced  (Interdisciplinary  Center  for 
Biotechnology  Research,  ICBR,  University  of  Florida).  From  these  terminal  sequences, 
primer  sequences  were  selected  using  the  computer  software  package  OLIGO  5.0 
(National  Bioscience,  Plymouth,  MN).  Primers  were  synthesized  by  Lifetechnologies 
and  by  the  ICBR,  University  of  Florida.  The  primer  sequences  for  each  amplified  region 
are  presented  in  Table  3-1. 


PCR  Amplification  and  Restriction  Enzyme  Digests 


Amplification  conditions  were  optimized  for  the  primer  pairs  selected  from  each 
clone  (Table  3-1).  If  no  amplification  was  achieved  using  the  first  pair  of  selected 
primers,  another  pair  was  designed.  Amplifications  were  performed  in  50mM  Tris-HCl 
pH  9.2, 2.0  mM  MgCl2,  (1.5mM  MgCl2  was  used  for  regions  XPS-3  and  XPL-2),  16mM 
(NH4)2S04,  with  2.0  units  of  Taq  DNA  polymerase  and  0.25  units  of  Pwol  DNA 

polymerase.  Both  polymerases  were  taken  from  a  single  stock  mixture  (Boehringer 
Mannheim,  Indianapolis,  IN).  Region  XPS-3  was  amplified  with  one  unit  of  Taq  DNA 
polymerase  without  Pwol  DNA  polymerase.  Amplification  reactions,  for  both  standard 
and  long  PCR,  were  in  50ul  total  volume,  containing  200uM  of  each  dNTP  and  250nM 
of  each  primer.  Totals  of  lOOng  and  500ng  of  template  DNA  were  used  for  the  short  and 
long  regions,  respectively.  To  increase  the  specificity  of  the  amplification  and  the  yield 
of  long  fragments,  wax  mediated  Hot  Start  PCR  was  used  (Chou  et  al,  1992).  Briefly, 
for  each  long  PCR  reaction,  1/3  of  the  total  amount  of  water,  2/3  of  the  total  volume  of 
buffer,  the  nucleotides  and  the  primers  were  aliquoted  into  one  0.5ml  thin  walled 
microcentrifuge  tube  and  overlaid  with  35ul  of  liquid  wax  (MJ  Research,  Watertown 
MA).  The  mixture  was  heated  at  75°C  for  5min  and  cooled  to  5°C  for  another  5min. 
Then  2/3  of  the  total  amount  of  water,  1/3  of  the  total  amount  of  buffer,  the  template 
DNA  and  the  enzyme  mix  were  added  on  top  of  the  solidified  wax  layer,  and  the 
optimized  PCR  program  for  the  selected  pair  of  primers  was  run.  This  procedure 
substantially  reduced  the  formation  of  primer  dimers  and  increased  the  yield  of  PCR 
products.  Hot  Start  PCR  was  needed  only  for  the  long  amplifications.  The  short  region 
and  two  long  regions  (XPL-2  and  XPL-5)  were  amplified  with  standard  PCR  profiles 
(Table  3-1).  An  autosegmented  profile  (progressive  increase  of  the  extension  time  per 
cycle)  was  used  for  one  long  region  XPL-1  (Table  3-1). 


Screening  for  RFLPs 

To  screen  for  polymorphic  sites,  5(j.l  of  PCR  products  from  one  African  and  one 
European  sample  were  digested  by  the  addition  of  5ul  of  a  restriction  digest  mixture 
containing  2  units  of  the  respective  restriction  enzyme  and  restriction  digest  buffer  to  a 
final  concentration  of  IX,  as  suggested  by  manufacturers  (for  a  final  volume  of  lOul). 
The  restriction  enzymes  used  were  Alul,  Aval,  Avail,  BamHl,  Bell,  Bglll,  BsfXL,  Ddel, 
Ecom,  EcoRV,  HaeWL,  Hhal,  HincH,  Hinfl,  Hpal,  Kpnl,  Mspl,  Nci\  Ndel,  Nhel,  Nsil, 
Nspl,  Pvull,  Sail,  Sau3A,  Sau96,  Seal,  Smal,  Spel,  Sphl,  Sstl,  Xhol,  mdXmal. 
Restriction  fragments  were  separated  by  electrophoresis  in  2%  agarose  gels,  stained  with 
ethidium  bromide  and  visualized  over  UV  light.  After  this  first  screening,  the  locus 
enzyme  combinations  (LECs)  that  generated  different  banding  patterns  between  the 
African  and  the  European  sample  were  tested  on  a  larger  sample  size  consisting  of  48 
African  and  70  European  drones.  Alleles  were  named  with  the  clone  followed  by  the 
enzyme  and  a  letter  designating  the  allele  (e.g.  XPL-l/Alul-A).  Composite  alleles  were 
determined  for  loci  XPL-2  and  XPL-5,  whereby  each  individual  drone  was  assigned  a 
combination  of  letters  corresponding  to  the  individual  alleles  produced  by  each  enzyme  at 
a  particular  locus. 

Estimation  of  Allele  Frequencies  and  Associations. 

Drones  are  haploid  parthenogenetic  progeny  of  diploid  queens.  Therefore, 
frequencies  of  both  individual  and  composite  alleles  were  estimated  in  African  and 
European  populations  on  the  basis  of  queen  genotypes  obtained  from  analyzing  their 
drone  progeny.  Two  to  six  drones  per  colony  were  analyzed.  If  two  different  alleles 
were  found  in  drones  from  the  same  colony  (heterozygous  queen),  each  allele  was  scored 
once.  If  only  one  allele  was  found,  it  was  scored  once,  or  twice  if  found  in  at  least  six 
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drones  (98%  probability  the  queen  was  homozygous)  (Hall,  1990).  Contingency  tables 
with  intralocus  association  frequencies  of  individual  alleles  produced  by  pair  of  enzymes 
was  constructed,  for  example,  alleles  produced  by  Aval  and  HaeJR  at  locus  XPL-2. 
Association  frequencies  in  all  contingency  tables  were  tested  for  independence  using 
Fisher's  exact  test  with  the  function  STRUC  in  GENEPOP  vl  .2  (Raymond  and  Rousset, 
1995). 

Results 

Comparisons  of  RFLPs  Obtained  in  DNA  Amplified  by  Standard  PCR  and  Long  PCR 

Of  the  five  clones  selected  for  standard  short  PCR,  three  were  successfully 
amplified,  one  could  not  be  amplified,  and  one  generated  PCR  products  of  the  expected 
size  but  failed  to  amplify  in  30%  of  the  samples.  Of  the  five  regions  selected  for  long 
PCR,  three  generated  PCR  products  of  the  expected  size  and  two  failed  to  amplify,  even 
when  new  sets  of  primers  were  designed.  PCR  product  sizes  ranged  from  1 .6  to  2.0  Kbp 
for  the  short  regions  and  from  5.5  to  9.0  Kbp  for  the  long  regions.  The  total  number  of 
restriction  enzymes  tested  ranged  from  26  to  3 1  for  both  short  and  long  regions.  Between 
the  initial  pair  of  African  and  European  samples,  5  of  78  LECs  (6.4%)  tested  from  the 
short  regions  generated  different  band  patterns  compared  to  25  of  85  LECs  (29.4%)  from 
the  long  regions.  After  a  large  number  of  samples  were  screened  for  those  enzymes  that 
revealed  polymorphisms  during  the  initial  screening,  only  1  of  78  of  the  initial  number  of 
LECs  (1.3%)  tested  from  the  short  regions  revealed  specific  nuclear  DNA  RFLPs 
between  African  and  European  populations,  compared  to  7  of  85  LECs  (8.2%)  from  the 
long  regions. 
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Informative  RFLPs  were  found  in  one  short  region,  XPS-3,  with  AM.  More  were 
found  within  the  long  regions:  XPL-1  with  AM,  XPL-2  with  Aval,  and  HaeTR,  and  XPL- 
5  with  HaeTR,  Ddel  and  Spel.  Fifty  individual  alleles  were  generated,  of  which  22  (44%) 
were  African-specific,  13  (26%)  were  European-specific,  two  (4%)  were  predominant  in 
European  samples,  and  13  (26%)  were  common  to  both  African  and  European 
populations.  Of  the  50  individual  alleles,  only  seven  came  from  the  short  region  XPS-3. 
For  loci  XPL-2  and  XPL-5,  composite  alleles  were  determined  from  the  combination  of 
individual  alleles.  Forty-three  composite  alleles  were  found  of  which  22  (51%)  were 
African-specific  and  19  (46%)  were  European-specific  (Table  3-2). 

Description  of  DNA  Polymorphisms 

Collective  frequencies  of  African  and  European-specific  individual  and  composite 
alleles  for  each  locus  in  the  United  States  and  South  Africa  are  presented  in  Table  3-3. 
XPS-3  is  the  only  clone  less  than  two  Kbp  in  length  analyzed  in  this  study  that  generated 
distinguishing  RFLPs  between  African  and  European  samples.  Detailed  characterization 
of  this  polymorphic  locus  will  be  reported  in  Chapter  4.  At  locus  XPL-1,  AM  generated 
eight  alleles  of  which  one  was  European-specific  (H),  one  was  European-predominant 
(F),  and  five  were  African-specific  (B,  C,  D,  E  and  G)  (Table  3-2  and  Figure  3-1).  At 
locus  XPL-2,  Aval,  and  HaeTR  revealed  polymorphisms.  Aval  generated  six  alleles  of 
which  one  was  European-specific  (F),  one  was  European-predominant  (£>),  and  two  were 
African-specific  (A  and  Q  (Table  3-2  and  Figure  3-2).  HaeTR  generated  twelve  alleles  of 
which  six  were  European-specific  (G,  H,  I,  J,  K  and  L)  and  three  were  African-specific 
(D,  E  and  F)  (Figure  3-2).  At  locus  XPL-2,  a  total  of  25  composite  alleles  were  found  of 
which  1 1  were  African-specific  and  14  were  European-specific  (Table  3-2).  At  locus 


XPL-5,  HaeJR,  Ddel  and  Spel  revealed  polymorphisms.  HaeTR  produced  10  alleles  of 
which  two  were  European-specific  (G  and  H)  and  six  were  African-specific  (A,  B,  C,  E,  F 
and  I)  (Table  3-2  and  Figure  3-3).  Ddel  generated  five  alleles  of  which  two  (B  and  E) 
were  European-specific  and  two  were  African-specific  (C  and  D)  (Table  3-2  and  Figure 
3-3).  Spel  produced  only  two  alleles  corresponding  to  the  presence  and  absence  of  one 
restriction  site  in  the  entire  region  (Figure  3-3).  Ddel  allele  B  and  Spel  allele  B  may  be 
west  European-specific.  From  each  of  5  colonies  of  east  (A.m.  ligustica)  and  west  (A.m. 
mellifera)  European  bees,  two  workers  were  selected  and  tested  (drones  from  Europe 
were  not  available  for  this  study).  Ddel  allele  B  and  Spel  allele  B  were  found  in  all  the 
west  European  bees  but  not  in  the  east  European  or  African  bees. 

At  locus  XPL-5,  a  total  of  18  composite  alleles  were  found  of  which  1 1  were 
African-specific  and  seven  were  European-specific  (Table  3-2). 

For  each  LEC,  the  diploid  genotypic  banding  pattern  was  constructed  for  every 
pairwise  combination  of  alleles  that  might  be  found  in  worker  bees,  and  the  pairs  of 
alleles  that  would  be  indistinguishable  from  others  were  determined.  Only  two  LECs 
(XPL-2  HaeHl  and  XPL-5  HaeHl)  produced  indistinguishable  genotypes,  probably  as  a 
result  of  the  complex  banding  pattern  and  the  large  number  of  low  molecular  weight 
fragments  generated  with  HaeHl.  For  LEC  XPL-2  HaeHl,  the  following  pairs  of 
genotypes  could  not  be  distinguished:  CC  and  CL,  IK  and  IL,  AJ  and  AK,  and ,  JK  and  JL. 
For  LEC  XPL-5  HaeHl,  genotypes  JJ  and  DJ  could  not  be  distinguished. 

Allele  Associations. 

To  determine  the  association  of  individual  alleles  generated  by  pairs  of  different 
enzymes  at  a  locus,  Fisher's  exact  tests  were  used  on  contingency  tables.  The  absence  of 


50 


independence  between  alleles  generated  by  pair  of  enzymes  at  a  particular  locus  was 
statistically  significant  for  all  five  comparisons  (P<0.01).  For  example,  at  locus  XPL-2, 
drones  with  allele  Aval-D  are  likely  to  have  allele  HaeJU-G  (Table  3-2).  High  frequency 
European  alleles  were  generally  associated  with  each  other.  For  example,  for  locus  XPL- 
2,  strong  associations  were  found  between  HaeJR-G,  and  Aval-D  (Table  3-2).  Other 
associations  between  European  alleles  were  found,  but  most  corresponded  to  low 
frequency  alleles.  African-specific  alleles  were,  in  general,  associated  with  common 
alleles  and  not  with  other  African-specific  alleles  (Table  3-2). 

Discussion 

As  an  alternative  to  Southern  blot  analysis  with  cloned  probes  and  as  an 
enhancement  of  standard  PCR-RFLP,  we  have  implemented  the  long  PCR  approach  in 
our  search  for  nuclear  DNA  markers  specific  to  African  and  European  bees.  As  would  be 
expected,  and  as  shown  by  our  results,  increasing  the  size  of  the  amplified  fragment 
increases  the  probability  of  detecting  RFLPs.  From  a  total  of  seven  LECs  that  generated 
informative  polymorphisms  between  African  and  European  samples,  six  came  from  long 
regions  compared  to  one  from  a  short  region.  The  sizes  of  the  long  PCR  products  in  this 
study  were  about  the  same  as  cloned  probes  used  in  earlier  studies,  and  alleles  were 
obtained  that  showed  similar  distributions  in  African  and  European  populations  (Hall, 
1992a;  McMichael  and  Hall  1996). 

The  amplification  of  long  DNA  fragments  is  more  demanding  than  standard  PCR. 
Higher  quality  and  larger  amounts  of  DNA  are  required.  Amplification  of  long  DNA 
fragments  is  more  sensitive  to  the  reaction  conditions  and  may  require  the  use  of  wax 
mediated  Hot  Start  PCR  to  increase  the  specificity  of  primer  annealing  and  the  yield  of 
PCR  product.  The  chances  of  false  primer  binding  increase  as  the  size  of  the  fragment  to 
be  amplified  increases,  and,  therefore,  longer  primers  are  required  to  increase  the 


specificity  of  annealing.  Autosegmented  programs  are  recommended  for  the 
amplification  of  long  DNA  fragments,  which  require,  after  an  initial  number  of  cycles, 
the  progressive  increase  of  the  extension  time  (usually  20  seconds  per  cycle).  A 
subroutine  system  can  be  used  instead,  in  which  the  extension  time  is  increased  100 
seconds  every  five  cycles.  This  modification  provides  comparable  results  to  the 
autosegmented  programs  and  is  an  acceptable  substitute  for  thermocyclers  that  do  not 
have  the  programming  capability  to  perform  autosegmented  profiles. 

For  the  regions  amplified  in  this  study,  a  number  of  different  fragment  patterns 
were  generated  by  single  restriction  enzymes  and  are  referred  to  as  individual  alleles. 
The  combinations  of  these  different  alleles  found  in  individual  drones  are  called 
composite  alleles.  A  composite  allele  is  African  or  European-specific  if  one  or  more  of 
its  individual  alleles  is  specific  and  its  other  individual  alleles  from  different  enzymes  are 
common.  Thus,  compared  to  individual  alleles,  a  greater  proportion  of  composite  alleles 
is  considered  specific,  but,  because  of  the  larger  number,  frequency  data  from  the  same 
number  of  individuals  are  not  as  accurate.  As  expected,  a  lack  of  independence  was 
found  between  pairs  of  individual  alleles  from  different  enzymes  at  the  same  locus.  In 
some  cases,  pairs  or  triplets  of  African  or  European-specific  individual  alleles  were 
strongly  associated,  which  reinforced  the  specificity  of  the  composite  alleles. 

As  employed  in  this  study,  haploid  drones  are  valuable  for  initially  identifying 
different  RFLP  patterns  and  deciphering  the  allelic  relationships  of  restriction  fragments. 
However  the  use  of  drones  for  determining  allele  frequencies  in  populations  has  some 
drawbacks.  Drones  are  usually  available  only  during  prosperous  times  and  not  from  all 
colonies.  Alleles  in  drones  from  the  same  colony  represent  those  of  only  one  individual, 
the  mother  queen.  Thus,  sample  sizes  are,  in  effect,  limited  to  the  number  of  colonies. 
Because  honey  bee  queens  mate  with  a  number  of  drones  from  different  colonies  (Adams 
et  ah,  1977;  Page,  1986),  worker  progeny  are  more  representative  of  the  local  population, 
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although  the  queen  genotype  would  still  be,  and  that  of  some  paternal  drones  may  be, 
over-represented. 

It  was  shown  that  screening  regions  of  DNA  amplified  by  long  PCR  is  an 
effective  approach  for  finding  useful  RFLP  markers.  As  discussed  later  (Chapter  5)  more 
narrow  sections  which  contain  the  informative  restriction  sites  were  isolated  from  the 
long  regions,  so  that  they  could  be  detected  using  standard  PCR.  The  resulting  simpler 
patterns  make  the  genotypes  of  workers  easier  to  identify,  and  the  less  demanding 
protocol  facilitates  the  analysis  of  large  numbers  of  samples.  More  samples  and 
populations  tested  in  the  future  should  result  in  more  robust  determinations  of  the 
specificity  and  frequencies  of  both  the  individual  and  composite  alleles. 
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Table  3-2.  Number  of  individual  and  composite  alleles  in  African  and  European 
populations.    Percentage  values  relative  to  each  population  are  indicated  in  parentheses. 


Locus  (Enzymes) 


Allele 


Population 


South  Africa 


United  States 


XPS-3  (Alul) 


B 

C* 

D* 

G 
H 


4(18) 
8(35) 
1(4) 
1(4) 
.3.(12). 
5(22) 
1(4) 


17(68) 
8(32) 


n 


23 


25 


XPL-104M)  B*  4(16) 

C*  3  (13) 

D*  3  (13) 

E*  1  (4) 

 G*  L(4)._ 


//t  2  (7) 

1  (4)                    21  (72) 
_A  1 1  (46)  6(21) 


n  24  29_ 

XPL-2  (Aval-HaeUT) 

A* -A  3  (12) 
A*-B  1  (4) 

A*-C  3  (12) 
A  -E  1  (4) 

fl-C  3  (12) 
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Table  3-2  continued. 

B-D 

* 

B-E 

B-F*  2  (8) 

C*-C  2  (8) 

D*-C  1  (4) 

__E-C_.   .114) 


1(3) 

5-/t 

2(6) 

B-tf 

2(6) 

D^-A 

1(3) 

dI-c 

1(3) 

15  (45) 

3(9) 

1(3) 

1(3) 

1(3) 

eJ 

1(3) 

1(3) 

5(20) 

1(3) 

1(4) 

3(9) 

n 

25 

34 

Locus  XPL-5  (i/aem-D^I-5^eI) 

* 

^  -A-A 

5(25) 

*  * 
^  -£>  -A 

2(10) 

* 

B  -A-A 

1(5) 

* 

C  -A-A 

1(5) 

D-A-A 

2(10) 

1(4) 
1(4) 


Table  3-2  continued. 


D-C  -A 

1(5) 

* 

D-D  -A 

1(5) 

* 

E  -A-A 

4(20) 

* 

F  -A-A 

1(5) 

* 

I  -A-A 

1(5) 

J-A-A 

 .L.C5). 

G^-A-A 

G^-B^-B^ 

G^-E^-A 

rf-A-A 

rf-E^-A 

J-B^-B^ 

 n  20 

(*):  Individual  African-specific  allele, 
(f ):  Individual  European-specific  allele. 
(%):  Individual  European-predominant  allele 
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Figure  3-1.  Alleles  found  in  region  XPL-1  with  Alul.  African-specific,  European- 
specific  and  European-predominant  alleles  are  indicated  by  a  solid  line,  a  dashed  line,  and 
a  dotted  line  respectively.  Bands  characteristic  of  African-specific  alleles  are  indicated  by 
closed  arrows.  Bands  characteristics  of  European-specific  and  European-predominant 
alleles  are  indicated  by  open  arrows.  European-specific  allele  H  and  European- 
predominant  allele  F  both  have  a  980  bp  fragment.  African  alleles  are  characterized  by 
fragments  of  1 120bp  (B),  1200bp  (Q,  1650bp  (D),  900bp  (E)  and  2200bp  (G).  Alleles 
without  an  African  or  European-specific  band  are  distinguished  by  their  unique  banding 
pattern.  Stl=Molecular  size  standard  (100  bp  DNA  ladder.  Lifetechnologies  Gibco 
BRL). 
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Figure  3-2.  Alleles  found  in  region  XPL-2  with  Aval  (I)  and  HaeW  (H).  Symbols  are  as 
indicated  in  Figure  3-1 .  (I)  European-specific  and  predominant  alleles  found  with  Aval 
are  characterized  by  a  600bp  fragment.  African  alleles  are  characterized  by  fragments  of 
1020bp  (A)  and  1500bp  (C).  (H)  The  European-specific  alleles  found  with  HaeHL  are 
characterized  by  the  presence  of  a  550bp  (G  and  H)  or  a  570bp  (/,  J  and  K)  fragment. 
African-specific  alleles  are  characterized  by  fragments  of  900bp  (£>)  and  300bp  (E). 
African-specific  allele  F  does  not  have  specific  fragments  and  is  distinguished  from  the 
others  by  its  unique  banding  pattern. 
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Figure  3-3.  Alleles  found  in  region  XPL-5  with  HaeUl  (I),  Ddel  (H)  and  Spel  (IE). 
Symbols  are  as  indicated  in  Figure  3-1.  (I)  European-specific  alleles  found  with  HaeTR 
are  characterized  by  a  1  lOObp  fragment.  African-specific  alleles  are  distinguished  by 
fragments  of  910bp  (5),  930bp  (1),  580bp  (Q  and  430bp  (F).  (H)  European-specific 
alleles  found  with  Ddel  are  distinguished  by  fragments  of  600  and  1090  bp  for  allele  B 
and  1975  and  880  bp  for  allele  E.  African-specific  alleles  are  characterized  by  a  1700  bp 
fragment  for  allele  C  and  a  580  bp  fragment  for  allele  D.  (HI)  Spel  produces  two  alleles 
of  which  one  (B)  is  probably  west  European-specific.  St2=Molecular  size  standard 
(lKbpDNA  ladder.  Promega). 


CHAPTER  4 
CHARACTERIZATION  OF  LOCUS  XPS-3 


Introduction 

This  chapter  describes  the  characterization  of  locus  XPS-3,  the  only  polymorphic 
locus  of  less  than  2kbp  found  when  comparing  long  PCR  and  standard  PCR  procedures  to 
search  for  nuclear  DNA  markers  (Chapter  3).  RFLP  analysis  of  locus  XPS-3  digested 
with  Alul  revealed  alleles  specific  to  or  predominant  in  African  and  west  European  bees. 

Materials  and  Methods 

The  samples  used  to  determine  allele  frequencies  for  locus  XPS-3  were  as  listed 
in  Chapter  2.  DNA  was  isolated  from  thoracic  muscles  of  adult  honey  bees  as  described 
in  Chapter  2.  Primer  sequences,  amplification  and  digestion  conditions,  restriction 
cleavage  mapping,  frequencies  and  estimates  of  observed  and  expected  heterozygosity, 
were  described  in  Chapters  3  and  Chapter  5  without  further  modification. 

Results 

Alul  digests  of  locus  XPS-3  revealed  seven  alleles  of  which  five  were  African- 
specific  {A,  B,  C,  D,  and  E),  one  was  west  European-predominant  (H)  and  one  was 
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common  (Table  4-1  and  Figure  4-1).  African-specific  alleles  were  collectively  found  at  a 
62.9%  frequency  in  the  South  African  samples.  Allele  Alul-Hwas  found  fixed  in  the 
west  European  subspecies  of  A.m.mellifera  and  A.m.iberica,  was  found  at  low  frequency 
in  the  east  European  subspecies  A.m.ligustica  and  A.m.carnica  (2.4%  and  7.4%, 
respectively),  and  was  found  absent  in  the  A.m.caucasica  samples.  Alul-H  was  also 
found  at  low  frequencies  in  the  South  African  samples  (2.3%). 

A  restriction  cleavage  map  of  locus  XPS-3  is  shown  in  Figure  4-2.  All  seven 
alleles  resulted  from  the  unique  combinations  of  four  polymorphic  Alu\  sites.  African- 
specific  alleles  were  characterized  by  fragments  of  1050  (A),  250  (B),  550  (Q,  940  (D) 
and,  900bp  (£).  The  west  European-predominant  allele  (H)  was  characterized  by  the 
presence  of  both  a  lOObp  and  a  150bp  fragment  (Figure  4-1).  As  with  other  polymorphic 
loci  with  African-specific  alleles,  genetic  diversity  values  were  higher  in  the  South 
African  population  compared  to  the  European  populations  (Table  4-1). 

Discussion 

Markers  specific  to  African  and  west  European  groups  of  subspecies  have  been 
more  difficult  to  find  than  to  the  east  European  groups.  Generally,  the  number  of 
African-specific  alleles  has  been  greater  than  that  of  east  European-specific  or 
predominant  alleles  resulting  in  a  high  genetic  diversity  in  the  African  populations. 
However,  the  collective  frequencies  of  the  African  alleles  tend  to  be  lower,  reported  here 
(Chapters  2,  3,  4  and  6)  to  average  40.5%  compared  to  80.0%  for  east  European  alleles. 
The  62.9%  collective  frequency  for  the  African  Alul  alleles  at  locus  XPS-3  is  among  the 
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highest  found  in  this  study  and  is  comparable  to  those  for  other  polymorphic  loci  found 
previously  (Hall,  1998:  McMichael  and  Hall,  1996).  The  ability  to  discriminate  African 
and  west  European  bees  from  east  European  bees,  the  simple  banding  patterns,  and  the 
ease  of  standard  amplification  conditions  make  this  polymorphic  locus  particularly 
useful. 
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Figure  4-1:  Restriction  cleavage  map  of  the  polymorphic  Alul  sites  at  locus  XPS-3. 
Positions  of  the  sites  are  indicated  to  the  left  (in  bp).  Letters  in  parentheses  to  the  right, 
are  the  alleles  that  have  the  restriction  sites.  Numbers  in  brackets  correspond  to  the 
primer  sequences  used  to  amplify  locus  XPS-3  (See  Chapter  3).  Numbers  in  parentheses 
indicate  the  sizes,  in  bp,  of  the  fragments  generated. 
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Figure  4-2:  RFLP  patterns  of  locus  XPS-3  digested  with  AM.  M100:  Molecular  size 
standard.  African  (AF)  and  European  (EU)  specific  alleles  are  marked  with  an  open  and 
closed  line  respectively.  Arrows  to  the  right  indicate  fragment  sizes,  in  bp,  characteristic 
of  each  allele. 


CHAPTER  5 

CONVERSION  OF  LONG  PCR  MARKERS  TO  STANDARD  PCR  FORMAT 


Introduction 


Nuclear  DNA-RFLPs  that  distinguish  groups  of  African  and  European  honey  bee 
subspecies  have  been  successfully  found  by  amplifying  and  digesting  anonymous  regions 
of  DNA  with  a  number  of  restriction  endonucleases  (Chapter  3).  Seeking  informative 
RFLPs  in  long  regions,  amplified  by  modified  procedures,  has  been  more  productive  than 
in  the  short  regions  amplified  by  standard  PCR.  However,  several  problems  remain 
concerning  the  amplification  of  large  DNA  fragments  and  the  analysis  of  RFLP  patterns 
of  the  amplifiable  PCR  products.  First,  high  quality  DNA  is  required.  DNA  from  poorly 
preserved  samples  is  difficult  and  usually  impossible  to  amplify.  Second,  as  the  length  of 
the  amplifiable  DNA  region  increases,  usually  the  number  of  alleles  increases,  and  their 
restriction  fragment  band  patterns  become  more  complex.  The  genotypes  of  diploid 
individuals,  that  is  workers  and  queens  in  bees,  become  more  difficult  to  determine. 
Therefore,  conversion  of  long  PCR  RFLP  markers  to  a  standard  PCR  format 
(amplification  of  fragments  of  up  to  2kbp)  facilitates  the  amplification  of  DNA  samples 
and  the  interpretation  of  band  patterns.  Conversion  of  RFLP  markers  found  with  probes 
on  long  PCR  to  a  standard  PCR  format  requires  subcloning  into  sections  and 
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amplification  of  the  sections  containing  the  informative  polymorphic  sites.  This  process 
is  easier  for  polymorphisms  found  in  amplified  regions  compared  to  those  found  with 
probes. 

In  this  chapter,  I  discuss  the  conversion  of  the  long  PCR  markers,  previously 
described  in  Chapter  3,  to  a  standard  PCR  format,  and  the  characterization  of  these 
markers  using  a  larger  sample  size  of  worker  honey  bees  from  different  populations. 

Materials  and  Methods 

Honey  Bee  DNA  Samples,  Restriction  Cleavage  Maps,  Subcloning,  and  DNA 
Amplification. 

Samples  of  adult  honey  bee  workers  came  from  South  Africa,  Europe,  and  North 
America  as  previously  described  in  Chapter  2.  Conversion  of  long  PCR  to  a  standard 
PCR  format  was  done  for  three  loci,  each  of  which  have  RFLP  alleles  generated  with  one 
to  three  enzymes:  XPL-1  with  AM,  XPL-2  with  HaelB.  and  Aval,  and  XPL-5  with 
HaeW.,  Ddel  and  Spel.  Restriction  cleavage  maps  were  obtained  for  the  informative 
polymorphic  sites  and  of  other  flanking  enzyme  sites,  which  enabled  the  regions 
containing  the  polymorphic  sites  to  be  subcloned  and  subsequently  amplified  using 
standard  PCR  procedures.  Maps  were  generated  using  the  procedure  reported  by  Her  and 
Weinshilboum  (1995).  This  approach  is  illustrated  in  Figure  5-1  for  two  hypothetical 
alleles.  For  this  mapping  procedure,  40ul  of  PCR  products  from  samples  corresponding 
to  each  allele  at  a  locus  were  partially  digested  with  two  units  of  the  respective  enzyme  in 
IX  restriction  digest  buffer  and  water  to  a  final  volume  of  50uT  Partial  digestions  were 
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obtained  by  aliquoting  13|a.l  of  the  digestion  mixture  in  a  tube  containing  12ul  of  0.5M 
EDTA  and  5ul  of  electrophoresis  loading  buffer  (20%  Ficoll,  0.5%  bromophenol  blue 
Tris-EDTA  pH  7.0)  at  2.5,  7.5, 15.0  and  30.0  minutes  after  initial  incubation  at  37°C. 
Partially  digested  products  were  electrophoresed  at  4V/cm  for  14-16  hours  in  a  2% 
agarose  gel  and  transferred  to  a  nylon  membrane  (Genescreen  plus)  by  capillary  blotting 
(Sambrook  et  al,  1989).  The  oligonucleotide  primers  used  to  amplify  each  region  served 
as  probes.  Primers  were  end-labeled  in  25ul  with  50ng  of  primer,  1  unit  of  T4 
polynucleotide  kinase,  IX  exchange  reaction  buffer  and  25uCi  of  y32P  (Lifetechnologies, 
Gibco  BRL,  Gaithersburg,  MD).  Labeling  reactions  were  performed  at  37°C  for  30min. 
Unincorporated  nucleotides  and  other  inpurities  were  removed  using  anionic  columns 
(Nucleotide  removal  kit,  Qiagen,  Palo  Alto,  CA).  Blots  were  dried,  then  prehybridized 
for  two  hours,  and  hybridized  for  three  hours  at  37°C  with  ExpressHyb  hybridization 
solution  (Clonetech,  Palo  Alto,  CA).  Following  hybridization,  the  blots  were  given  two 
washes  of  200ml  each  of  2XSSC  +  0.1%SDS  and  two  washes  of  200ml  each  of  0.1XSSC 
+  0.5%SDS  at  room  temperature  (Sambrook  et  al,  1989).  Each  wash  was  for  30  min. 
Membranes  were  sealed  in  plastic  bags  and  placed  against  X-ray  film,  held  between 
intensifying  screens  at  -70°C  overnight.  For  each  locus,  hybridization  was  done  first 
with  the  forward  primer.  Blots  were  stripped  and  then  hybridized  with  the  reverse 
primer.  By  using  both  forward  and  reverse  primers  as  probes,  the  position  of  the 
restriction  sites  relative  to  each  end  of  the  section  can  be  determined.  This  confirms  the 
position  of  the  site  and  allow  for  more  accurate  mapping  of  the  sites,  especially  at  the 
termini  of  the  sections.  Sites  for  enzymes  found  in  the  polycloning  section  of  the  plasmid 
vector  pGem3Z  (Promega  Corporation,  Madison,  WI)  were  also  mapped  within  the 
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sections.  For  subcloning  of  honey  bee  DNA,  each  region  was  first  digested  with  the 
enzymes  having  sites  both  within  the  polycloning  section  of  pGem3Z  and  that  flanked  the 
polymorphic  sites.  Digested  fragments  were  separated  in  a  1%  agarose  gel  and  purified 
(Agarose  gel  extraction  kit,  Qiagen,  Palo  Alto,  CA).  Purified  fragments  were  then 
subcloned  in  the  plasmid  vector  pGem3Z  by  conventional  cloning  procedures  (Sambrook 
et  al,  1989).  Sequences  for  the  terminal  400  bases  of  each  subclone  were  obtained  and 
used  to  design  primer  sequences  with  the  software  package  OLIGO  5.0  (National 
Bioscience,  Plymouth,  MN).  Amplifications  were  optimized  for  each  region  for 
magnesium  ion  concentration  and  cycling  parameters  (annealing  temperature  and  time). 
To  check  for  the  presence  of  polymorphic  sites  in  each  region,  5ul  of  PCR  products  were 
digested  with  the  respective  enzyme  as  explained  in  Chapter  2.  Digested  PCR  products 
were  electrophoresed  at  3V/cm  in  a  2%  agarose  gel  and  visualized  by  ethidium  bromide 
staining  over  UV  light. 

Allele  Determination  and  Data  Analysis 

Alleles  found  for  the  regions  amplified  by  long  PCR  were  named  using  the 
enzyme  name  followed  by  a  letter  (See  Chapter  3).  Amplifiable  sections  were  identified 
using  numbers  and  sometimes  an  abbreviation  based  on  their  position  relative  to  the 
entire  polymorphic  locus  (for  example,  XPL-2S1  for  section  1  of  locus  XPL-2,  and  XPL- 
2Slint  for  section  1  internal  of  locus  XPL-2).  RFLP  analysis  of  subclones  from  the 
original  polymorphic  loci  results  in  simpler  allelic  band  patterns,  but  because  usually 
only  a  subset  of  the  original  restriction  sites  are  present,  distinction  between  some  alleles, 
originally  found  with  long  PCR,  was  not  possible.  If  an  allele  could  be  distinguished 
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from  the  others,  it  retained  the  same  description  used  for  the  long  PCR  alleles. 
Otherwise,  groups  of  alleles  with  the  same  banding  patterns  in  a  subcloned  section  were 
named  with  an  UX"  followed  by  a  number.  For  locus  XPL-2  and  XPL-5,  one  to  three 
sections  were  amplified  to  determine  the  allele  identity:  two  sections  for  LECs  XPL2 
Aval  and  XPL-5  HaeUl  and  three  sections  for  LEC  XPL-5  Ddel.  When  two  or  three 
sections  were  used,  the  original  allele  identity  was  determined  by  the  combination  of 
band  patterns  found  for  the  different  sections.  For  each  locus  and  population,  the 
frequency  of  each  allele  and  genotype  was  determined  and  used  to  calculate  expected  and 
observed  heterozygosity  (Nei,  1978  and  Nei,  1987). 

Results 

For  each  polymorphic  locus  first  amplified  with  long  PCR  (XPL-1,  XPL-2  and 
XPL-5),  one  to  three  different  sections  were  reamplified.  These  sections  contained  the 
informative  polymorphic  sites  of  a  particular  restriction  enzyme  and  were  of  sizes 
amplifiable  by  standard  PCR  protocols:  0.48  kbp  to  3.19  kbp  (Table  5-1).  Amplification 
conditions  for  all  the  sections  within  each  locus  are  described  in  Table  5-1. 

Characterization  of  Polymorphic  Sections 

Locus  XPL-1.  As  detailed  in  Chapter  3,  RFLP  analysis  of  locus  XPL-1  revealed 
polymorphisms  with  the  restriction  endonuclease  Alul.  Eight  alleles  were  found  initially 
of  which  one  was  European-specific  and  one  was  European-predominant.  The  original 
European  alleles  were  all  characterized  by  the  presence  of  a  980  bp  fragment.  To 
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subclone  the  section  with  the  polymorphic  sites  responsible  for  this  European- 
predominant  fragment,  locus  XPL-lwas  divided  in  three  different  sections  (Figure  5-2). 
As  shown  in  Figure  5-2,  these  sections  are  defined  by  the  position  of  an  EcoSJ  site  at 
position  3 146  and  a  BamW  site  at  position  5852.  Positions  of  the  Alul  sites  that 
generated  the  different  alleles  were  also  mapped  (Figure  5-2).  The  Alul  site  at  position 
4990  is  present  in  the  European-specific  and  predominant  alleles  and  absent  in  the  other 
alleles  initially  found  with  long  PCR.  Therefore,  the  section  containing  this  polymorphic 
site  (XPL-1S2,  between  the  EcoRl  and  BamW  sites)  was  subcloned  and  amplified  for 
RFLP  analysis  using  standard  PCR  protocols.  To  verify  the  presence  of  the  Alul  site  in 
the  amplifiable  XPL-1S2,  DNA  from  drone  samples  containing  the  different  alleles, 
initially  found  using  long  PCR,  were  amplified  and  digested  with  Alul  (Figure  5-3). 
From  these  digests,  three  sub-alleles  were  found,  compared  to  eight  alleles  in  the  entire 
region;  one  African-specific  (G),  one  common  (XI  which  correspond  to  the  original 
alleles  Alul-A,  B,  C,  D  and  E  found  using  long  PCR)  and  one  European-predominant  (X2 
which  correspond  to  the  original  alleles  Alul-F  and  H).  Sub-allele  G  is  characterized  by  a 
2050bp  fragment.  The  European-predominant  sub-allele  is  characterized  by  fragments  of 
150  and  700  bp  (Figure  5-3). 

Frequency  data  of  the  three  different  sub-alleles  found  in  XPL-1S2  is  presented  in 
Table  5-2.  The  African-specific  sub-allele  G  was  not  found  in  any  of  the  worker  samples 
but  was  found  at  low  frequency  in  the  African  drone  samples  (Chapter  3).  Frequencies  of 
the  European  sub-allele,  X2,  were  high  in  the  east  European  samples  (97.4%,  and  92.5% 
in  A.m. ligustica,  and  A.m.carnica,  respectively)  and,  in  the  near  east  samples  (70.0%  in 
A.m.  caucasica).  Sub-allele  X2  was  also  found  in  high  frequencies  in  the  U.S.A  samples 
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(78.7%).  In  the  west  European  samples,  sub-allele  X2  was  absent  in  A.m.  iberica  but 
present  at  a  low  frequency  in  A.m.mellifera  samples  (4.5%)  and  in  the  South  African 
samples  (3.0%).  Therefore,  sub-allele  X2  was  predominantly  east  European. 

Locus  XPL-2.  Analysis  of  locus  XPL-2  revealed  polymorphisms  with  Aval  and 
HaeUl  (Chapter  3).  Two  sections  from  this  locus  were  amplified.  The  first  section, 
XPL-2Slint,  is  an  830bp  fragment  within  section  1  located  between  the  terminal  Pstl  site 
and  the  BamHl  site  at  position  1280.  A  second  section,  XPL-2S2,  is  a  3195bp  fragment 
located  between  the  EcoRl  site  at  position  2765bp  and  the  Sphl  site  at  position  6230 
(Figure  5-4).  Sub-alleles  for  Aval  were  identified  by  amplifying  both  XPL-2Slint  and 
XPL-2S2  (Figure  5-4  and  Table  5-3).  In  section  XPL-2Slint,  three  sub-alleles  could  be 
distinguished:  Aval-A  (African-specific),  characterized  by  an  830bp  fragment,  Aval-X2 
(European-predominant)  characterized  by  a  600bp  fragment,  zsxAAval-Xl  (common) 
(Figure  5-5  I).  XPL-2S2  has  three  different  sub-alleles:  Aval-C  (African-specific), 
characterized  by  a  1505bp  fragment,  and  Aval-Xl  and  X2  (common)  (Figure  5-5  II). 
Identification  of  the  six  alleles  originally  found  with  long  PCR  was  possible  by  looking  at 
the  unique  sub-allele  combinations  found  in  both  polymorphic  section:  XPL-2Slint  and 
XPL-2S2.  Combinations  A  JC1,  X1:C,  X2-X1,  X2-X2, Xl'Xl,  X1'X2  (first  letter  is  the  sub- 
allele  found  in  XPL-2Slint  and  the  second  letter  is  the  sub-allele  found  in  XPL-2S2) 
correspond  to  the  alleles  Aval-A,  C,  D,  F,  B  and  E,  respectively  revealed  using  long  PCR 
(Table  5-3).  For  HaeUl,  only  one  section  was  amplified  (XPL-2Slint).  Three  alleles  of 
the  original  12  alleles  found  with  long  PCR  were  clearly  distinguishable:  HaeTH-E,  F 
(African-specific)  and  HaelR-K  (European-specific).  The  remaining  alleles  were 


grouped  by  four  banding  patterns  of  which  two,  HaeHl-X2  and  X3,  were  European- 
specific  (Figure  5-6)  corresponding  to  four  of  the  original  six  European-specific  alleles 
seen  with  long  PCR  (HaelR-G,  H,  I  and  J).  The  European-specific  sub-alleles  HaelR- 
X2,  X3  and  K  were  characterized  by  a  520bp  fragment  (Figure  5-6).  Original  alleles 
HaeJR-D  (African),  L  (European)  and  C  (common)  were  seen  as  XI,  and  alleles  HaeTR-A 
and  B  as  X4. 

Sub-allele  frequencies  for  alleles  revealed  by  Aval  and  HaeHI  in  African  and 
European  populations  are  shown  in  Table  5-2.  Alleles  Aval-D  and  F  were  found 
collectively  at  frequencies  of  77.1%  and  97.6%  in  the  east  European  subspecies 
A.m.ligustica  and  A.m.carnica  respectively,  were  absent  in  west  European  samples  and 
present  at  very  low  frequencies  in  A.m.scutellata  (South  African  samples)  (1.1%)  and  in 
A.m.caucasica  (near-East)  (27.5%).  These  alleles  were  found  at  intermediate  frequencies 
in  samples  from  the  United  States  (65.6%).  African-specific  alleles  Aval-A  and  C  were 
found  at  a  51.1%  collective  frequency  in  the  A.m.scutellata  samples  and  were  absent 
from  the  east  and  west  European  samples.  However,  Aval-A  was  found  at  low 
frequencies  in  the  U.SA  samples  (3.7%).  Sub-alleles  HaeHi-X2,  X3  and  K  were  found  at 
collective  frequencies  of  76.2,  94.2  and  92.5%  m  A.m.ligustica,  A.m.carnica  and 
A.m.caucasica  samples,  respectively.  HaeW-K  was  found  only  in  A.m.caucasica 
(30.0%)  and,  therefore,  appeared  to  be  specific  to  this  subspecies.  These  sub-alleles  are 
east  European-specific,  found  to  be  absent  from  west  European  and  the  South  African 
samples.  African-specific  sub-alleles  HaeYR-E  and  F  were  found  at  a  low  frequency  in 
the  South  African  samples  (17.1%)  and  were  absent  in  all  the  European  and  the  U.SA 
samples. 
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Locus  XPL-5.  Analysis  of  locus  XPL-5  revealed  polymorphisms  with  the 
enzymes  HaeUi,  Ddel,  and  Spel  (Chapter  3).  Locus  XPL-5  was  divided  into  four 
sections.  The  first  section,  XPL-5S1  is  a  1499  bp  fragment  located  between  the  Hindi 
sites  at  positions  630  and  2129.  The  second  section  is  a  900  bp  fragment  located  between 
the  Hindi  sites  at  positions  3029  and  2129.  The  third  section  is  a  1501  bp  fragment 
located  between  the  Hindi  sites  at  position  4530  and  3029,  and  the  fourth  is  a  900  bp 
fragment  located  between  the  Xbal  site  at  position  5430  and  the  Hindi  site  at  position 
4530.  To  detect  the  polymorphic  HaeHi  site,  two  sections  were  amplified:  section  1 
extended  (XPL-5 Slxt),  a  1980bp  fragment  located  within  section  1  between  the  Pstl  site 
and  the  Hindi  site  at  position  2129,  and  section  1  terminal  (XPL-5Slter)  a  485bp 
fragment  that  overlaps  the  Hindi  site  near  section  2  and  the  polymorphic  HaeUl  site  in 
XPL-5Slxt  at  position  1880  that  is  absent  in  the  east  European  bees  (Figure  5-7).  Both 
sections  were  amplified  which  subdivided  the  complicated  banding  patterns  originally 
found  with  HaeUl.  From  the  original  ten  alleles  detected  using  long  PCR,  five  sub-allele 
groups  were  revealed  (Yl,  Y2,  Y3,  C  and  F)(resulting  from  the  combination  of  banding 
patterns  found  in  both  XPL-5  Slxt  and  XPL-5  S  Iter)  of  which  three  are  African-specific 
(72,  C  and  F)  and  one  is  European-predominant  (J3)  (Figure  5-8 1  and  II,  and  Table  5-2). 
The  "Y"  notation  has  been  used  in  this  case  to  avoid  confusion  with  the  conventional  "X" 
notation  used  previously  for  groups  of  banding  patterns.  In  this  particular  case,  groups  of 
banding  patterns  were  found  in  both  sections  and  were  named  with  an  "X"  followed  by  a 
number.  The  final  sub-allele  designation  is  designated  with  a  "Y".  The  African-specific 
sub-alleles  are  characterized  by  fragments  of  974,  600  and  450bp  respectively  (in  digests 
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of  XPL-5Slxt)  and  the  European-predominant  sub-allele  (Y3)  is  characterized  by  a  485bp 
fragment  detected  in  digests  of  XPL-5Slter  (Figure  5-8  IT). 

Sub-allele  Y3  is  East  European-predominant,  found  at  a  frequency  of  92.9  and 
98.1%  in  A.m.ligustica  and  A. m.carnica,  respectively,  at  100.0%  in  A.m.caucasica,  and  at 
intermediate  frequencies,  71.9%,  in  the  USA  samples.  Sub-allele  Y3  was  found  at  a  1.1% 
frequency  in  A.m.scutellata  (South  African). 

Identification  of  Ddel  alleles  (Ddel  A,B,C,D  and  E)  was  done  by  amplifying 
three  sections  of  locus  XPL-5:  section  1  extended,  section  2  (XPL-5S2),  a  715  bp 
fragment  located  between  the  Hindi  sites  at  positions  2129  and  3029,  and  section  4 
(XPL-5S4),  a  935  bp  fragment  located  between  the  Hindi,  site  at  position  4530  and  the 
terminal  Xbal  site  (Figure  5-9).  All  the  alleles  found  with  long  PCR  (Ddel  A,  B,  C,  D  and 
E)  could  be  identified  by  the  banding  patterns  detected  in  one  or  a  combination  of  two  or 
more  of  the  three  sections.  Polymorphisms  in  section  XPL-5  Slxt  distinguishes  allele  B 
from  the  others,  section  XPL-5  S2  distinguishes  alleles  C  and  E  and,  section  XPL-5  S4 
distinguishes  allele  D.  Allele  A  is  distinguished  from  the  others  by  its  banding  pattern 
combination  in  the  different  sections  (XPL-5 Slxt,  S2  and  S4)  (Figure  4-9).  Allele  B  is 
west  European-predominant  found  fixed  in  both  A.m.mellifera  and  A.m.iberica,  at  4.8  and 
1 1.5%  frequencies  in  the  east  European  samples  of  A.m.ligustica  and  A. m.carnica, 
respectively,  and,  at  a  frequency  of  5.0%  in  the  A.m.caucasica  samples.  Allele  E  is  east 
European-specific,  found  at  a  low  frequency  in  A. m.ligustica  (14.3%)  and  at  an 
intermediate  frequency  in  A. m.carnica  (48. 1%),  and  absent  in  the  African  and  west 
European  samples  (Table  5-2). 


Polymorphisms  revealed  by  Spel  were  detected  by  amplifying  one  section  (XPL- 
5S3)  of  locus  XPL-5.  Two  sub-alleles  were  detected  resulting  from  the  presence  (sub- 
allele  B)  or  absence  (sub-allele  A)  of  one  restriction  site  at  position  3290.  Sub-allele  B 
was  European-specific,  predominantly  found  in  the  west  European  subspecies:  in 
A.m.mellifera  at  a  73.1%  frequency  and  in  A.m. iberica  at  a  33.3%  frequency.  Sub-allele 
B  was  also  found  at  low  frequencies  in  the  east  European  subspecies  A.m.ligustica  (1.2%) 
and  A.m.carnica  (15.4%),  and  absent  in  A. m.caucasica  (Table  5-2). 

Discussion 

Making  long  PCR  markers  analyzable  by  standard  PCR  has  facilitated  their  use 
for  population  genetic  studies.  Amplification  of  small  regions  of  DNA  do  not  require 
high  quantity  and  quality  of  DNA  as  with  long  PCR.  Because  the  sizes  of  the  amplified 
regions  are  reduced,  the  resulting  band  pattern  of  the  alleles,  initially  identified  using 
long  PCR,  are  simplified.  Thus,  the  genotypes  of  workers  and  queen  (diploid)  honey 
bees  can  be  more  easily  determined.  The  quality  of  DNA  suitable  for  long  PCR  can  be 
consistently  obtained  when  extracted  from  entire  pupae  but  not  from  thoracic  muscles  of 
adult  bees.  The  difficulty  of  obtaining  such  high  quality  of  DNA  from  adults  restricts  the 
applicability  of  the  long  PCR  markers.  With  the  conversion  of  these  markers  to  a 
standard  PCR  format,  DNA  from  adult  honey  bee  thoracic  muscles  can  be  easily 
analyzed.  Moreover,  the  simplicity  of  the  banding  pattern  and  the  small  size  of  the 
fragments  amplified  make  it  easier  to  analyze  diploid  individuals.  Thus,  these 
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modifications  make  a  large  number  of  workers  from  different  populations  amenable  for 
analysis. 

To  determine  the  genotypes  of  worker  honey  bees,  one  to  three  sections  had  to  be 
amplified,  and  patterns  obtained  in  each  section  were  combined  to  obtain  a  final  allele 
identification.  In  testing  unknown  samples,  usually,  identification  of  all  the  original 
alleles  was  possible.  However,  sometimes  when  a  large  number  of  alleles  per  locus  was 
found  in  the  population,  the  identification  of  both  alleles  in  an  individual  worker  using 
the  standard  PCR  modification  was  not  possible,  because  the  phase  (cis  versus  trans)  of 
the  polymorphisms  relative  to  each  other  cannot  be  determined  (as  described  in  Chapter 
2).  For  example,  in  the  locus  enzyme  combination  XPL-2/Aval,  the  identification  of  each 
allele  is  done  by  amplifying  two  sections:  XPL-2Slint  and  XPL-2S2.  If  a  worker's 
genotype  for  XPL-2S 1  int  is  XI 1X2  and  for  XPL-2S2  is  XI 1X2,  two  possibilities  exist:  a 
cis  configuration  on  homologous  chromosomes  will  be  X1X1  on  one  chromosome  (Allele 
Aval-B)  and  X2X2  on  the  homologous  chromosome  (Allele  Aval-F)  (Table  5-3);  a  trans 
configuration  will  have  X1X2  (Allele  Aval-E)  on  one  chromosome  and  X2X1  (Allele 
Aval-D)  on  the  homologous  chromosome.  Thus,  the  individual  could  be  either  BIF  or 
EID.  In  this  example,  both  alleles  B  and  E  are  common  alleles  and,  D  and  F  are 
European-predominant,  therefore,  an  incorrect  identification  of  these  alleles  would  not 
affect  the  collective  frequencies  of  the  European  markers.  These  ambiguous  cases  were 
rarely  found  in  our  samples,  and,  typically,  most  of  the  alleles  could  be  correctly 
identified. 

As  in  any  PCR-based  method,  erroneous  readings  can  also  be  attributed  to 
preferential  amplification  of  certain  alleles  or  to  the  presence  or  artifactual  or 


heteroduplex  DNA  (Anglani  et  al.,  1990;  Arnheim  and  Erlich,  1992;  Hare  et  al,  1996). 
Partial  restriction  digestions  can  also  confuse  allele  identification. 

The  methodological  and  analytical  simplifications  of  the  markers  found  with  long 
PCR  procedures  greatly  increase  their  value  for  both  research  and  regulatory  purposes. 
These  genetic  markers  are  thus  an  important  addition  to  the  collection  of  markers 
becoming  available,  which  can  be  used  to  monitor  the  spread  of  the  African  bees  in  the 
United  States. 


80 


81 


O 

a 
o 


W 

H 

W 


C 

_o 

a 
o 


a 
o 


W 

to 


o 
< 


8 

CO  ~ 

Pi  « 
5 


i 

< 


f  I 


8 

•i  I 

CO  g 


I  s 


"O  c/j 

O  o 

a  -2 

P  C/5 


C3  5 

1 1 

o  S 


o 


O  O 

O  O 

O  O 

o  d 


o 


o  o 

r-  (N 

O  Os 

d  d 


m 


SO  -t 
O  On 

d  d 


o 
o 
o 
o 


m  in 

m  "si- 

OS  o 

d  d 


00 

so 


m 


-H  £j  fN 


O  O  O 
O  O  O 
O  O  (N 

r--  m 
odd 


o  o  t-» 
m  o  oo 

M  VI  Cl 
Os  <— <  i 

odd 


too 

t--  in  in 
Os  O  O 

odd 


in  os  oo 

m  o  so 

Tt  os  oo 

o  o  o 

odd 


00 

00 

00 

m 

m 

SO 

so 

00 

m 

00 

00 

o 

(N 

CO 

en 

d 

d 

d 

d 

d 

OS 

t— t 

(N 

OS 

o 

o 

O 

OO 

SO 

ro 

so 

I/O 

Os 

O 

o 

o 

o 

d 

d 

d 

d 

d 

o 


o  o  o  o 
in  o  o  <n 

M  iTl  b 
(N  O  SO  O 

dodo 


SO  SO  00 
00  t— 
(N  "<3-  (N 
Os  O  O 

odd 


<n 

m 


r-- 


so 
oo 

(N 
(N 

d 


O 

o 
o 
o 


o 
o 
o 
o 


o 

m 

r- 

SO 

SO 

r~ 

m 

in 

O 

ro 

o 

so 

O 

d 

d 

d 

d 

r- 

SO 


M  M  VI 
m  H  *  rt 

U  Tt  o  o 
odd 


o 
d  d 


S 


d 


5l 
o 


82 


O  O  00 
m  O  00 

r-  w> 
<N  w->  u-> 

odd 


CN  On  On 
NO  (N 
h  CI 

ON  —>  — ' 

odd 


Ti-  On  P- 

r~  V) 
m  m 

odd 


•o 

D 

.9 

I 
o 
o 

i 

u  U 

3 
n 
H 


h  r- 
vo  (N  m 
u->  oo  r- 

«  W 

odd 


m  ro  On 
— i  0O  (N 

O  NO  NO 

odd 


o 


o  o 

lO  o 

(N  O 

NO  m 

O  O 


o 

P* 

o 


NO 


On 

NO  I/O 

(N 
00 

o  d 


00  ON 

CO 

o  o 

d  o 


(N  on 
o 

o 
©  o 


CN  On 

NO  — < 

(N  — I 

(N  O 

o  o 


O 

o 
o 
o 


o 
o 
o 
o 


nO 

NO 


00  o 

o  o  o 

ro  M  ^ 

_!  q  — ; 
o  o 


o 

NO 

WO 

00 

r- 

<N 

o 

On 

00 

NO 

NO 

NO 

o 

o 

o 

(N 

NO 

NO 

NO 

o 

v—i 

NO 

NO 

NO 

o 

d 

d 

d 

d 

d 

d 

d 

s 


i 


o  o 
lO  o  m 
M  in  -h 
on 

odd 


c->  o 
<n  r-  h 

O  O 

On  rn  rn 
odd 


on  o  oo 

^  o  •-< 

no  O  *r> 
t-m  >n  -ft 

odd 


00 

cn 

o 

o 

NO 

wo 

m 

oo 

o 

CN 

ON 

m 

m 

rn 

NO 

d 

d 

d 

d 

d 

51  S| 


p. 


ON. 


NO 

04 


00  (N 

O  ON 

oo  ^ 

on  q 

d  d 


NO 

00  — < 

<n  r- 

On  O 

O  O 


o 
o 
o 
o 


NO 


Tf  NO 

in  O  on 

i-J  q  q 

d  d 


o 
o 
o 
o 


On  ~h 

O  On 

—  oo 

r-  (N 

d  d 


OO  H 

^ 

o  o  oo 
odd 


I 


2  gfct  *•  22  R 


ft 


83 


m  r-~ 

oo  r- 

o  o 

d  d 


d  d 


o 

00  — H 

d  d 


00  00 

f-  OS 

r-  <n 

d  d 


j 

I 

in 


o 


o 
o 
«n 
o 


o 
o 

in 

OS 


^  OO  00 

in  o  m 

h  OO  o 

^   Tj-  Tf 

d  d  d 


*0  On  m 

t  (N  On 

^  st  O 

O  ^  M 

odd 


in 


SO 
so 


Cs|  00  O 

r-H 

r~-  oo 
ro  — < 
odd 


ro  m 
o  o 

(N  O 

d  d 


so 
rs) 
OS 

r- 


-Si 


o  o  o 
o  o  m 
m  o  os 

O.-H  O 

odd 


oo       m  m 

o  m  oo  fN 

00  h  co  os 

i-h  «n  »n 

dodo 


On  NO  O 

csi  r-  cn 

^  ^  (N 

<— '  o  m  cn 
d  d  d  d 


o 
o 
o 
o 


o 
o 
o 
o 


00  <N  <N 

i— i  i— i  i— i  On 

CO  (N  M 

t-j  m  no  so 

d  d  d  d 


m 
o 

rsi 


r-  m 
m  oo  o 
o  CO 
o  m  m 

odd 


U  U  u  SB 


o 
o§ 


(N  00 
so  C") 

£  3 

M  00  -h 

d  d 


i-h  On 
n  00 
OS  O 

d  d 


f-  co 
no  ro 
«n  no  m 
^  no  m 

d  d 


<N  oo 

On  O 
cn  no  cn 
~  <N  r- 

d  d 


o  o 

^    Tf  NO 

^  OS  o 
NO  no  cn 

d  d 


o 

NO  O 

28 


# 

in  1 
hi  ' 


OO  \t 
CO  f-  o 

<n  o  no 
■h  m  ts 

odd 


os  oo  in 
h  n  n 

-h  <N  <N 
O  O  O 

odd 


cn  t>  -*t 
en  no  Ti- 
en no  tj- 

CO  NO 

odd 


oo  m 

o  oo 

cn  cn 

r-  in 

d  d 


o  o 

no  cn 

O  O 

CO  Tf 

d  d 


U 


3?  a? 


03 


£3 


*  8  S  & 


84 


Table  5-3.  Allele  determination  as  revealed  by  amplification  of  multiple  sections  using 
standard  PCR  procedures.  


Locus  bnzyme 

Allele  combinations 

XT              A  11  1 

New  Allele 

Long  PCR  a 

Combination 

round  in  each 

designation (1) 

section 

2\JrLi-Z.IAVai 

Mint/;>z 

A  1X1 

A* 

XI IC 

C* 

C* 

V")  /  V7 

A2  1  Al 

/  A2 

XI 1 XI 

B 

XI 1 X2 

E 

YPT  -S/DW/»T 

Q1  Yt/Q9/Q4 

J  1  Xli  oZ/  ot 

xiicixi 

C* 

C* 

XI 1X1  ID 

Z)¥ 

B 1X1 1X1 

5§§ 

XII  El  XI 

£« 

XI 1X1 1X1 

A 

XVL-5/HaeUl 

Slxt/Slter 

Yl 

X2IX1 

Yf 

CI  XI 

C* 

C* 

FIX1 

F* 

XI 1X2 

75§ 

(1):  Original  allele  found  at  locus  XPL-5  that  could  not  be  distinguished  by  the 
combination  of  banding  patterns  from  RFLP  analysis  of  section  1  extended  (Slxt)  and 
section  1  terminal  (Slter)  were  grouped.  Different  group  of  alleles  were  named  with  a  Y 
followed  by  a  number. 

(§):  East  European-specific  or  predominant  allele. 
(§§):  West  European-predominant  allele. 
(¥):  African-specific  or  predominant  allele. 
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Figure  5-2.  Restriction  cleavage  map  of  locus  XPL-1  with  the  restriction  enzyme  Alul. 
Numbers  indicate  the  position  of  monomorphic  (open  symbols)  and  polymorphic  (closed 
symbols)  Alul  sites.  Numbers  in  parentheses  indicate  sizes  of  fragments.  Letters  in 
parentheses  are  the  alleles  with  the  polymorphic  sites.  Numbers  enclosed  in  brackets  are 
of  the  primers  used  to  amplify  section  2  and  their  sequences  are  listed  in  Table  5-1 . 
Locus  XPL-1  was  divided  in  three  sections:  section  1  (XPL-1  SI)  from  the  Pstl  site  at 
position  0  to  the  EcoRl  site  at  position  3 146.  Section  2  (XPL-1  S2)  from  the  EcoRl  site  at 
position  3146  to  the  BamYU.  site  at  position  5852,  and  section  3  (XPL-1  S3)  from  the 
BamHl  site  at  position  5852  to  the  Xbal  site  at  position  8120.  I)  Relative  position  of 
Alul  sites.  H)  Position  of  Alul  sites  and  the  size  of  the  resulting  fragments  in  XPL1-S2. 
Indicated  with  an  asterisk  is  the  site  that  results  in  the  east  European  predominant  alleles 
FmdH. 
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Figure  5-3.  RFLP  patterns  revealed  byAlul  digests  of  section  2  of  locus  XPL-1  (XPL- 
1S2).  Three  patterns  are  seen:  the  first  pattern  {XI)  which  includes  most  of  the  African 
alleles  (B,  C,  D  and  E)  and  the  common  allele  (A).  The  second  pattern  corresponds  to  the 
African-specific  allele  G,  and,  the  third  pattern  (X2)  corresponds  to  the  European  alleles 
(F  and  H).  The  pattern  (X2)  is  characterized  by  the  presence  of  a  700  and  a  150bp 
fragment  resulting  from  one  restriction  site  located  in  the  850bp  fragment  found  in 
pattern  XI.  St:  Molecular  size  standard  (lOObp  ladder,  Lifetechnologies,  Gibco  BRL, 
Gaithersburg,  MD).  African  (AF)  and  European  (EU)  alleles  are  marked  with  an  open 
and  closed  line  respectively.  Letters  in  the  first  line  are  the  allele  designations  initially 
found  in  long  PCR  amplifications  and,  just  above  the  picture,  the  new  sub-allele 
designations  based  on  standard  PCR  amplification  of  this  one  region. 


Figure  5-4.  Restriction  cleavage  map  of  locus  XPL-2  with  the  restriction  enzymes  Aval 
and  HaeHl.  Numbers  and  symbols  to  the  left  and  right  of  each  map  correspond  to  the 
data  for  Aval  and  HaeHl,  respectively.  Circles  and  squares  correspond  to  Aval  and 
HaeHl  sites  respectively.  (Polymorphic  sites  are  indicated  by  a  closed  symbol  and 
monomorphic  sites  are  indicated  with  an  open  symbol).  Letters  enclosed  in  parentheses 
are  the  alleles  with  the  polymorphic  site.  Numbers  enclosed  in  brackets  are  of  the 
primers  used  to  amplify  each  region  and  their  sequences  are  listed  in  Table  5-1 .  Numbers 
in  parentheses  indicate  fragment  sizes  in  base  pairs.  Locus  XPL-2  is  divided  in  three 
sections:  Section  1  from  the  Pstl  site  to  the  EcoBI  site  (position  2765),  section  2  from  the 
EcoBl  site  to  the  Sphl  site  (position  6230)  and  section  3  from  the  Sphl  site  to  the  Xbal 
site  (7130).  The  section  between  the  Pstl  and  the  BamYH  site  is  named  section  1  internal 
(Slint)  as  it  includes  only  half  of  the  original  section  1.  I)  Relative  positions  of 
restriction  sites  for  Aval  (left),  HaeHl  (right)  and  the  enzymes  used  to  subclone  fragments 
of  sizes  amplifiable  by  standard  PCR.  II)  Positions  of  both  monomorphic  and 
polymorphic  Aval  and  HaeHl  sites.  Indicated  with  one  asterisk  are  those  polymorphic 
sites  that  generate  European-predominant  alleles,  and  with  two  asterisks  those  that 
generate  African-specific  alleles.  For  Aval,  the  polymorphic  sites  located  at  positions 
550  and  3835  generate  African-specific  alleles.  The  absence  of  the  Aval  site  at  position 
550  results  in  allele  Aval-A  and  the  absence  of  the  site  at  position  3835  results  in  allele 
Aval-C.  The  presence  of  the  site  located  at  position  1 150  results  in  the  generation  of  the 
east  European-predominant  alleles  Aval-D  and  F.  For  HaeHl,  two  sites  are  responsible 
for  the  African-specific  allele:  one  at  position  520  and  the  second  at  position  670.  Alleles 
HaeHl-E  and  F  lack  the  former  and  HaeHl-F  lacks  the  latter.  The  absence  of  the  HaeHl 
site  located  at  position  1 120  results  in  the  east  European-predominant  alleles  HaeHl-G 
and  H.  HI)  Regions  amplified,  using  standard  PCR  protocols,  that  contain  most  of  the 
polymorphic  sites  for  Aval  and  HaeHl. 
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Figure  5-5.  RFLP  patterns  revealed  by  Aval  digests  of  section  1  internal  (Slint)  and 
section  2  (S2)  at  locus  XPL-2.  I)  Aval  digests  of  section  1  internal  (XPL-2Slint).  Three 
patterns  are  revealed:  The  first  pattern  (A)  is  African-specific  and  is  characterized  by  an 
830  bp  fragment.  The  second  pattern  (XI)  is  common  and  the  third  pattern  (X2)  is 
European-predominant  and  is  characterized  by  fragments  of  600  and  90  bp.  H)  Aval 
digests  of  section  2  (XPL-2S2).  For  XPL-2R2,  three  patterns  are  seen:  XI  and  X2 
(common)  and  C  (African-specific),  characterized  by  a  1505  bp  fragment.  The  unique 
combination  of  patterns  revealed  in  XPL-2Slint  and  XPL-2S2  is  used  to  determine  the 
allele  designation  originally  found  in  RFLP  analysis  of  the  entire  locus.  African  (AF) 
and  European  (EU)  alleles  are  marked  with  an  open  and  closed  line,  respectively.  Letters 
in  first  line  are  the  allele  designations  initially  found  in  long  PCR  amplifications,  and  just 
above  each  picture,  the  new  sub-allele  designations  based  on  standard  PCR  amplification. 
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Figure  5-6.  RFLP  patterns  revealed  by  HaeTR  digests  of  section  1  internal  of  locus  XPL- 
2  (XPL-2Slint).  African  sub-alleles  E  and  F  and  the  European  sub-allele  K  are 
distinguishable  from  the  others  by  their  unique  patterns.  Four  additional  patterns  were 
seen:  XI  and  X4  (common)  andA2  and  JO  (European).  European-specific  patterns  are 
distinguished  by  the  presence  of  a  520bp  fragment.  African-specific  sub-allele  E  is 
characterized  by  a  260bp  fragment.  African  (AF)  and  European  (EU)  alleles  are  marked 
with  an  open  and  closed  line  respectively.  Letters  in  the  first  line  are  the  allele 
designations  initially  found  in  long  PCR  amplifications,  and  just  above  each  picture,  to 
the  new  sub-allele  designations  based  on  standard  PCR  amplification. 


Figure  5-7.  Restriction  cleavage  map  of  locus  XPL-5  with  the  restriction  enzymes 
HaeUl,  Ddel  and  Spel.  Numbers  and  symbols  to  the  left  of  each  map  correspond  to  Ddel 
and  Spel  sites.  Numbers  and  symbols  to  the  right  correspond  to  HaeUl  sites. 
Polymorphic  and  monomorphic  sites  are  indicated  by  a  closed  and  open  symbol 
respectively.  Circles  correspond  to  Ddel  sites,  triangles  to  Spel  sites  and  squares  to 
HaeUl  sites.  Letters  in  parentheses  are  the  alleles  that  contain  the  particular  restriction 
sites.  Numbers  in  brackets  are  the  primers  used  to  amplify  each  region.  Primer 
sequences  are  in  Table  5-1.  Numbers  in  parentheses  are  to  fragment  sizes.  I)  Relative 
position  of  the  Ddel,  Spel  and  HaeUl  and,  the  Hincll  sites  used  to  subclone  locus  XPL-5 
into  several  sections  amplifiable  by  standard  PCR  protocols.  H)  Relative  positions  of 
monomorphic  and  polymorphic  Ddel,  Spel  (left  side  of  the  map),  and  HaeUl  sites  (right 
side  of  the  map).  HI)  Sections  of  locus  XPL-5  amplified  with  standard  PCR  protocols. 
Locus  XPL-5  was  divided  in  four  sections:  section  1  between  the  two  Hindi  sites  at 
positions  630  and  2129;  section  2  between  the  Hindi  sites  at  positions  3029  and  2139; 
section  3  between  the  Hindi  sites  at  positions  4530  and  3029  and,  section  4  between  the 
Xbal  site  at  position  5430  and  the  Hindi  site  at  position  4530.  Sites  that  revealed  east 
European-predominant  alleles,  west  European-predominant  alleles,  and  African  alleles 
are  indicated  by  an  asterisk,  a  (+),  and  two  asterisks,  respectively.  For  Ddel,  the  site  at 
positions  2595  is  found  in  the  east  European-specific  allele  Ddel-E,  and  the  site  at 
position  1 1 00  in  the  west  European-predominant  allele  Ddel-B.  African-specific  alleles 
are  determined  by  the  presence  of  Ddel  sites  at  positions  5300  (Ddel-D)  and  2830  (Ddel- 
C).  For  Spel,  the  presence  of  one  site  at  position  3290  results  in  the  west  European- 
predominant  allele  Spel-B.  For  HaeUl,  three  sites  located  at  positions  200,  450  and  950 
are  responsible  for  African-specific  alleles.  The  absence  of  the  HaeUl  site  located  at 
position  1774  is  characteristic  of  the  east  European-predominant  alleles  HaeUl-G  and  H. 
Ddel  alleles  were  determined  by  amplifying  three  sections:  section  1  extended  (Slxt), 
section  2  (S2)  and  section  3  (S3).  Spel  alleles  are  determined  by  amplifying  one  section: 
section  3  (S3).  HaeUl  alleles  are  determined  by  amplifying  two  sections:  Slxt  and 
section  1  terminal  (Slter).  The  polymorphic  site  that  determins  the  banding  pattern  for 
the  European  predominant  alleles  HaeUl-G  and  H  is  included  in  S  lxt.  However,  the 
absence  of  the  site  at  position  1880  in  alleles  HaeUl-G  and  //results  in  a  200bp  fragment 
that  is  indistinguishable  from  the  200bp  fragment  generated  by  the  presence  of  and 
HaeUl  the  site  at  position  200.  Therefore,  a  small  region  located  at  the  terminal  end  of 
Slxt,  and  containing  the  HaeUl  polymorphic  site  (at  position  1880)  is  amplified  and  used 
to  determine  the  presence  or  absence  of  this  site  responsible  for  the  European- 
predominant  alleles  HaeUl-G  and  H. 
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Figure  5-8.  RFLP  patterns  revealed  by  HaeUl  digests  of  section  1  extended  and  section 
1  terminal  of  locus  XPL-5.  I)  Section  1  extended  (XPL-5Slxt),  and  H)  section  1  terminal 
(XPL-5Slter).  Four  patterns  are  seen  in  digests  of  XPL-5  SI  xt:  XI  (common),  and  X2,  C 
and  F  (African-specific)  which  are  characterized  by  fragments  of  974  (X2),  600  (Q,  and 
450  (F)  bp,  respectively.  In  digests  of  XPL-5Slter,  two  patterns  are  seen:  XI  (common), 
and^f2  (European).  X2  is  characterized  by  a  485bp  fragment.  African  (AF)  and 
European  (EU)  alleles  and  sub-alleles  are  marked  with  an  open  and  closed  line, 
respectively.  Letters  in  first  line  are  to  the  allele  designations  initially  found  in  long  PCR 
amplifications,  and  just  above  each  picture,  the  new  sub-allele  designations  based  on 
standard  PCR  amplification  relative  to  each  region.  Indicated  at  the  bottom  are  the  new 
sub-allele  designations  resulting  from  a  combination  of  RFLP  patterns  in  both  sections. 


Figure  5-9.  RF LP  patterns  revealed  by  Ddel  digests  of  sectionl  extended  (Slxt),  section 
2  (S2)  and  section  4  (S4)  of  locus  XPL-5.  I)  Patterns  revealed  in  section  1  extended 
(XPL-5Slxt),  II)  section  2  (XPL-5  S2)  and,  m)  section  4  (XPL-5  S4).  Four  of  the  initial 
five  alleles  (B,  C,  D  and  E)  can  be  distinguished  in  separate  regions.  Allele  Ddel-A  is 
distinguished  by  the  unique  banding  pattern  combination  revealed  in  all  three  sections. 
Allele  Ddel-B  is  distinguished  by  RFLP  patterns  of  XPL-5  Slxt  and  is  characterized  by  a 
1 100  and  a  630bp  fragment.  Alleles  Ddel-C  and  E  are  distinguished  in  RFLP  patterns  of 
XPL-5 S2  and  are  characterized  by  fragments  of  635  and  80  bp  (Ddel-C)  and  fragments  of 
400  and  315  bp  (Ddel-E).  Allele  Dde\-D  is  distinguished  by  RFLP  patterns  of  XPL-5 S4 
and  is  characterized  by  a  445bp  fragment.  The  common  allele,  {Ddel-A)  was 
characterized  by  the  combination  of  banding  patterns  in  the  three  sections.  African  (AF) 
and  European  (EU)  alleles  are  marked  with  an  open  and  closed  line  respectively.  Letters 
in  first  line  are  the  allele  designations  initially  found  in  long  PCR  amplifications,  and  just 
above  each  picture,  the  new  allele  designations,  relative  to  each  region,  based  on  standard 
PCR  amplification. 
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CHAPTER  6 
A  MODIFICATION  OF  THE  AFLP  PROTOCOL 
APPLIED  TO  HONEY  BEE  (APIS  MELLIFERA  L.)  DNA. 

Introduction 

Amplified  fragment  length  polymorphism  (AFLP)  (Vos  et  al.,  1995)  is  a  recently 
developed  PCR-based  technique  that  provides  genetic  markers  for  fingerprinting, 
mapping,  and  studying  genetic  relationships  among  organisms  (Becker,  et  al,  1995;  Hill, 
etal.,  1996;  Janssen,  etal,  1996;  Majer,  et  al,  1996;  Thome,  etal,  1996;  Vos  and 
Kuiper,  1997).  This  technique  is  popular  because  it  detects  a  high  amount  of 
polymorphism  and  is  reproducible. 

The  AFLP  protocol  developed  by  Vos  et  al,  1995  involves  digesting  genomic 
DNA  with  two  restriction  enzymes:  a  frequent  and  an  infrequent  cutter.  Following 
digestion,  specially  designed  adapters  for  each  restriction  enzyme  are  ligated  to  the  ends 
of  the  DNA  fragments.  Primers  are  made  complementary  to  the  adapter  sequences,  plus 
the  enzyme  sequences,  and  up  to  three  randomly  selected  nucleotides  at  the  3 '  end.  The 
primers  are  used  to  amplify  the  DNA  fragments  generated  by  the  restriction  enzymes  in 
two  reactions:  primers  with  one  extra  base  at  the  3'  end  of  the  enzyme  sequences  are  used 
in  a  preamplification,  and  primers  with  three  extra  bases  are  used  in  a  second 
amplification.  One  of  the  primers  in  the  second  reaction  is  radioactively  labeled.  The 
amplified  products  are  resolved  in  polyacrylamide  gels  and  exposed  to  film  for 
autoradiography. 
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Simplification  of  the  original  AFLP  protocol  would  facilitate  the  analysis  of  large 
number  of  samples  as,  for  example,  in  population  genetic  studies.  Alternative  procedures 
have  been  proposed  to  eliminate  the  use  of  radioactivity  and  polyacrylamide  gels,  but 
DNA  transfer  to  blots  followed  by  a  colorimetric  or  chemoluminescent  detection  are  still 
required  (Latorra  et  ah,  1994;  Muller  et  ai,  1996).  In  this  report,  we  present  an  AFLP 
protocol  simplified  and  optimized  for  honey  bee  (Apis  mellifera  L)  DNA.  This  procedure 
consists  of  three  steps:  1).  Digestion  of  genomic  DNA  with  a  single  restriction  enzyme 
and  ligation  to  one  type  of  adapter  are  accomplished  together  in  one  reaction.  2).  The 
digested-ligated  DNA  is  amplified  in  a  single  reaction.  3)  The  AFLP  products  are 
resolved  in  agarose-synergel  and  visualized  by  ethidium  bromide  staining. 

Amplification  conditions  were  optimized  for  concentrations  of  honey  bee 
template  DNA,  magnesium  ions,  primer,  and  Taq  DNA  polymerase.  This  procedure 
reduces  the  number  of  steps  and  expensive  reagents  required  in  the  original  AFLP 
protocol.  The  modified  method  detects  a  considerable  amount  of  polymorphism, 
sufficient  for  many  purposes,  but  not  as  much  polymorphism  as  the  original  procedure. 
With  this  modified  AFLP  procedure,  useful  honey  bee  DNA  markers  were  found. 

Materials  and  Methods 

DNA  Samples. 

DNA  was  isolated  from  drone  honey  bees  as  described  in  Hall,  1990.  Drones  are 
haploid  and  thus  display  simpler  banding  patterns  of  amplified  products  than  do  workers 
(diploid  females).  European  drone  samples  came  from  three  locations  in  the  United 
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States  (Florida,  Kansas  and  Tucson),  and  African  samples  came  from  three  locations  in 
South  Africa  (Louis  Trichardt,  Johannesburg  and  Pretoria). 

Digestion-Ligation  Reactions. 

EcoRI  Adapters  were  prepared  by  adding  equimolar  amounts  of  each  strand  (Vos 
et  al.,  1995),  and  the  mix  was  incubated  at  90°C  for  1  minute  followed  by  a  reduction  of 
1°C  every  20  seconds  until  room  temperature  was  reached.  Digestion  of  the  genomic 
DNA  and  ligation  with  the  adapters  was  done  in  one  reaction.  Two  ug  of  genomic  DNA 
were  added  to  20  ul  of  lOmM  Tris.HAc  pH  7.5,  lOmM  MgAc,  50mM  KAc,  5mM 
dithiotreitol  (DTT),  50ng/ul  bovine  serum  albumin  (BSA)  and  lmM  ATP  as  suggested 
by  Vos  et  al.,  1995.  Twenty  units  of  EcoBI,  one  unit  of  T4  DNA  ligase  and  50  pmol  of 
the  EcoRl  adapter  were  added  and  the  mixture  was  incubated  at  37°C  for  five  hours.  To 
remove  unincorporated  adapters  and  residual  enzymes,  80ul  of  water,  lOOu  of  7.5M 
NtLt.Ac  and  600(4.1  of  95%  ethanol  were  added.  The  annealed  adapter  was  precipitated, 
the  pellet  washed  with  500ul  of  75%  ethanol,  air  dried  for  10  minutes,  resuspended  in 
100  ul  of  lOmM  Tris.HCl,  0.1  mM  EDTA  (TE)  and  stored  at  -20°C.  Enzymes  were 
obtained  from  Boehringer  Mannheim  (Indianapolis,  IN).  Primers  were  synthesized  by 
Lifetechnologies,  Gibco  BRL(Gaithersburg,  MD). 
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DNA  Amplification  and  Electrophoresis. 

As  described  in  the  Results  section,  amplification  conditions  were  optimized  for 
concentrations  of  four  reagents:  template  honey  bee  DNA,  magnesium  ions,  primer,  and 
Taq  DNA  polymerase.  PCR  reactions  were  performed  in  50ul  total  volume  in  thin 
walled  0.5ml  microcentrifuge  tubes  in  an  MJ  Research  PTC- 100  thermocycler 
(Watertown,  MA).  The  PCR  profile  consisted  of  one  step  at  95°C  for  1  min  followed  by 
35  cycles  of  94°C  for  lmin,  58°C  for  1  min  and  72°C  for  2  min  with  a  final  extension  at 
72°C  for  5min.  After  optimization  of  the  PCR  conditions,  African  and  European 
honeybee  drone  DNA  samples  were  amplified  with  three  different  primers  (EcoRI  + 
AGT,  EcoKl  +  ATA  and  EcoBl  +  AGC),  and  PCR  products  were  resolved  in  a  2.5% 
agarose-synergel  (0.7%  agarose  and  0.9%  synergel:  Diversified  Biotech,  Boston,  MA)  in 
0.5X  TPE  buffer  or  IX  TBE  buffer  (buffer  formulas  in  Sambrook  et  al,  1989.  Gels  were 
run  for  at  4V/cm  for  14  hours,  stained  with  ethidium  bromide,  and  viewed  over  UV  light 
to  visualize  the  PCR  products. 

Results 

Optimization  of  PCR  Reagents. 

DNA  concentration  was  tested  from  lng  to  lOOng  in  a  50ul  reaction  volume  (in 
triplicate).  Above  lOng,  high  reproducibility  with  no  variability  in  the  banding  pattern 
was  seen,  but  intensification  of  some  bands  was  apparent  as  the  DNA  concentration  was 
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increased.  Magnesium  ion  concentrations  were  tested  from  0.5mM  to  3.5mM  (in 
duplicate).  At  0.5mM,  no  bands  or  a  very  small  number  of  bands  were  seen.  At  1.5mM, 
the  number  of  bands  increased,  and  differences  in  band  intensity  were  seen.  Results  were 
reproducible  between  2.0  and  2.5mM.  As  the  concentration  increased  above  that  level, 
no  new  bands  were  formed,  and  bands  began  to  disappear  until  at  3.5mM  no  bands  were 
seen.  Primer  concentrations  from  100  to  500nM  in  lOOnM  increments  were  tested.  Low 
molecular  weight  bands  were  not  amplified  at  lOOnM.  As  primer  concentration 
increased,  more  bands  were  visible,  and  low  molecular  weight  bands  increased  in 
intensity.  No  differences  in  banding  patterns  were  seen  when  primers  at  400  to  500nM 
were  used.  Taq  DNA  polymerase  concentration  was  tested  from  0.5  to  3.5  units  in  a  50ul 
reaction  volume.  No  amplification  was  seen  at  0.5units.  Amplification  results  were 
consistent  at  enzyme  concentrations  between  1.0  and  3.0  units.  At  3.5  units,  a  small 
amount  of  smear  was  formed.  Nucleotide  concentrations  from  200  to  400 uM  did  not 
result  in  differences  in  the  amplification  when  2.5mM  magnesium  ions  was  used.  Thus, 
the  optimized  conditions  for  amplification  using  honey  bee  DNA  were  set  at:  2.0mM 
MgCl2,  2.0  units  of  Taq  DNA  polymerase,  50ng  of  DNA,  400nM  of  primer,  and  250uM 
of  each  dNTPs  in  lOmM  Tris-HCl  pH  8.3,  50mM  KC1  in  a  50ul  reaction  volume. 

Positive  and  negative  controls  were  used  for  all  amplifications.  The  positive 
controls  included  undigested,  unligated  DNA  with  all  the  other  reagents,  which  checked 
for  any  possible  bands  that  could  result  from  primer  annealing  to  a  genomic  DNA  region 
and  not  to  the  adapters.  The  negative  control  was  a  reaction  with  only  the  reagents 
without  template  DNA  and  was  intended  to  detect  possible  artifacts  generated  from  the 
primers.  No  visible  bands  were  found  in  either  positive  or  negative  controls,  therefore, 
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all  the  bands  were  produced  from  the  annealing  of  primers  to  the  adapters.  Artifactitious 
bands  were  seen  in  the  negative  control  when  more  than  35  cycles  were  used  in  the  PCR 
profile  but  were  not  seen  in  any  of  the  reactions  that  included  template  DNA. 

AFLP  Using  Honey  Bee  DNA. 

Three  primers  were  used  to  amplify  DNA  of  samples  from  African  and  European 
honey  bees.  The  DNA  of  two  drones  from  each  of  six  colonies,  three  African  and  three 
European,  was  amplified  using  the  optimized  conditions.  The  amplified  products, 
separated  by  electrophoresis  and  stained  with  ethidium  bromide,  are  seen  in  Figure  6-1. 
Substantial  amount  of  polymorphism  was  seen  for  all  the  primers  used.  A  few 
monomorphic  bands  (one  or  two  per  primer)  were  produced.  Each  sample  had  a  unique 
banding  pattern,  thereby  showing  the  effectiveness  of  the  method  for  individual 
identification  (i.e.  fingerprinting).  For  example,  with  one  primer,  individuals  from 
different  colonies  could  be  distinguished  by  the  presence  of  unique  bands  (Figure  6-1 -A). 
African  or  European-specific  bands  were  also  found  (Figure  6-1 -A;  Figure  6-1-B). 

Discussion 



Population  genetic  studies  have  been  enhanced  with  genetic  markers  detected 
through  methods  based  on  the  polymerase  chain  reaction  (PCR),  such  as  microsatellites 
(Tautz,  1989),  restriction  fragment  length  polymorphisms  in  PCR  products  (PCR-RFLP) 
(Karl  et  al,  1992;  Karl  and  Avise,  1993),  random  amplified  polymorphic  DNA  (RAPD) 
(Williams  et  al.,  1990),  and  more  recently,  AFLP  (Vos  etal,  1995).  AFLP  uses  arbitrary 
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sequences  to  generate  polymorphic  amplified  fragments  from  restriction  enzyme  digested 
DNA.  Thus,  this  method  combines  the  simplicity  of  RAPDs  with  the  reliability  of 
RFLPs.  As  with  RAPD  analysis,  AFLP  does  not  require  prior  knowledge  of  the  DNA 
sequence.  The  longer  primers  used  for  AFLP  require  higher  annealing  temperature  and 
binding  is  more  specific  and  stable. 

The  simplified  AFLP  procedure  described  here  generates  a  substantial  amount  of 
polymorphism,  although,  not  at  levels  as  high  as  those  generated  with  the  original 
procedure.  The  amount  of  polymorphism  can  be  increased  by  making  primers  less 
specific,  (using  two  extra  bases  instead  of  three)  by  changing  the  enzyme,  (using  a  four 
base  cutter  instead  of  six)  and,  by  changing  the  source  of  DNA  (e.g.  diploid  honey  bee 
workers  rather  than  haploid  drones).  The  digestion  of  DNA  and  ligation  of  adapters  is 
completed  in  five  to  eight  hours.  After  the  digestion-ligation  step,  the  same  sample  can 
be  used  with  different  primers,  in  which  case  the  remaining  procedure  can  be  completed 
in  less  than  1 8  hours. 

In  optimizing  the  reaction,  we  found  a  low  sensitivity  to  the  concentrations  of 
DNA  and  other  PCR  reagents.  High  DNA  concentrations  in  AFLP  do  not  affect  the 
reproducibility  of  the  reaction,  but  the  reaction  may  need  to  be  optimized  for  DNA  from 
different  organisms.  In  negative  controls  (all  the  reagents  except  DNA),  artifactitious 
bands  were  detected  when  more  than  35  cycles  were  used.  When  genomic  DNA  was 
amplified,  these  bands  did  not  appear  even  after  35  cycles.  Possible  problems  with  AFLP 
include  the  coamplification  of  DNA  from  endosymbionts  such  as  bacteria  or  other 
microorganisms  associated  with  the  organism  of  interest.  Partial  digestion  of  the  DNA  in 
the  digestion-ligation  reaction  can  lead  to  erroneous  results  and,  therefore,  an  excess  of 


106 


restriction  enzyme  to  digest  the  DNA  to  completion  should  be  used.  Also,  comigration  of 
amplified  products  can  make  patterns  difficult  to  interpret. 

Markers  found  with  AFLP  are,  in  general,  dominant.  Codominant  markers  are 
more  valuable  for  population  genetic  studies,  because  homozygous  and  heterozygous 
individuals  can  be  distinguished.  Codominant  markers  can  be  found  with  AFLP  if  the 
polymorphic  restriction  sites  fall  within  fragments  sufficiently  short  to  be  amplified. 
Thus,  the  fragments  lacking  the  sites  and  the  fragments  resulting  from  cuts  at  the  sites  can 
be  seen.  By  modifying  this  AFLP  procedure  to  amplify  longer  DNA  fragments  (Barnes, 
1994;  Cheng  et  al.,  1994a  and  b),  the  possibility  of  finding  codominant  markers  would 
increase. 

The  AFLP  procedure  typically  displays  multiple  polymorphisms  which  can  serve 
as  genetic  fingerprints  of  individuals.  However,  the  system  can  also  be  used  as  a  source 
of  single  genetic  markers.  DNA  fragments  containing  polymorphic  sites  can  be  isolated 
from  the  agarose-synergel,  cloned,  sequenced  and  used  to  design  primers  that  amplify 
only  one  polymorphic  region. 


Figure  6-1.  Amplification  of  African  and  European  honey  bee  drone  samples  with 
primers  EcoW  +  AGT  (A),  EcoRI  +  ATA  (B)  and  EcoBl  +  AGC  (C).  M=  molecular  size 
standard  (lOObp  ladder,  Lifetechnologies,  Gibco  BRL).  Two  drones  per  colony  (1-6)  and 
three  colonies  of  each  bee  type  (C1-C3  for  the  African  bees  and  C4-C6  for  the  European 
bees)  were  used.  Note  the  few  monomorphic  bands  and  the  substantial  amount  of 
polymorphism  found  in  the  amplification  products  for  all  the  primers.  Bands 
characteristic  of  some  colonies  were  found,  for  example,  band  "a"  for  colony  1  and  band 
"b"  for  colony  2  and  4  using  primer  EcoRl  +  AGT.  African  or  European-specific  bands 
were  also  found,  for  example,  band  "c"  using  primer  EcoRI  +  AGT  is  European-specific, 
and  bands  "a"  and  "b"  using  primer  EcoBJ  +ATA  are  African-specific. 
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CHAPTER  7 

GENE  INTROGRESSION  AND  HYBRIDIZATION  STUDIES  OF  NEW  WORLD 
HONEY  BEE  POPULATIONS  USING  NUCLEAR  DNA  RFLP  MARKERS 


Introduction 

In  the  first  part  of  this  dissertation,  nuclear  DNA  RFLP  markers  specific  to  groups 
of  African,  East  and  West  European  subspecies  were  sought.  Different  methods  were 
employed,  including  Southern  blot  analysis  with  randomly  cloned,  anonymous,  honey 
bee  DNA  used  as  probes,  and  analysis  of  anonymous  regions  amplified  by  Long  PCR. 
Markers  found  were  made  analyzable  by  standard  PCR  protocols  to  facilitate  their  use. 
These  markers,  along  with  nuclear  DNA  RFLP  markers  previously  described,  were 
employed  to  study  honey  bee  populations.  The  main  objectives  of  the  work  reported  in 
this  chapter  were  to  genetically  characterize  New  World  African  and  European  honey  bee 
populations,  and  to  determine  levels  of  hybridization  and  patterns  of  gene  flow. 

Materials  and  Methods 

Honey  bee  DNA  Samples. 

Honey  bee  DNA  samples  were  workers,  adults  or  larvae  from  Old  and  New 
World  African  and  European  populations.  The  Old  World  African  samples  were  bees  of 
the  subspecies  A.m.scutellata  from  South  Africa.  The  Old  World  European  samples 
included  subspecies  from  three  major  groups:  west  European  (A.m.mellifera  and 
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A.m.iberica),  east  European  {A.m.ligustica  and  A. m.carnica),  and  near-eastern 
(A.m.caucasica).  Samples  from  Old  World  African  and  European  populations  have  been 
previously  described  (See  Chapter  2)  and  were  used  as  references  to  which  New  World 
populations  were  compared.  New  World  African  samples  included  bees  from  the  Rio 
Grande  Valley  in  Texas,  Mexico,  Honduras,  Costa  Rica,  and  Venezuela.  Samples  from 
Texas  were  collected  from  July  to  September  1995  and  from  March  to  May  1997,  in  what 
was  the  northernmost  limit  of  the  African  bee  front  in  the  United  States  and  were 
provided  by  O.  R.  Taylor  (University  of  Kansas,  Lawrence,  Kansas).  African  bees  were 
first  reported  in  this  area  in  1990.  Mexican  samples  were  collected  in  1989  by  H.  Glenn 
Hall  (University  of  Florida,  Gainesville,  FL)  from  swarm  bait  hives  in  Tapachula, 
Mexico.  Honduran  samples  were  collected  in  1989  by  ASC  from  four  different  locations: 
El  Paraiso,  Choluteca,  Olancho,  and  Francisco  Morazan.  Samples  from  Costa  Rica  were 
collected  in  1989  by  H.  Glenn  Hall  and  H.  Arce  (Universidad  Nacional  de  Costa  Rica) 
near  Cerro  de  la  Muerte  and  San  Isidro  del  General.  Samples  from  Venezuela  were 
collected  in  1986  and  1988  and  were  provided  by  R.  Hellmich  Jr.,  J.  Villa,  A.  Collins, 
and  T.  Rinderer  (USDA-ARS,  Baton  Rouge,  Louisiana).  Details  of  number  of 
individuals  and  colonies  used  in  each  of  the  New  World  populations  are  provided  in 
Table  7-1 .  The  New  World  European  samples  included  bees  from  a  collection  of  bee 
stocks  across  the  United  States  (Page  et  al,  1982;  Severson  et  al,  1986)  and  from  Kansas 
(provided  by  O.R.  Taylor,  University  of  Kansas,  Lawrence,  KS).  DNA  extraction  was 
performed  as  described  in  the  Material  and  Methods  section  of  Chapter  2. 
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Mitochondrial  and  Nuclear  DNA  Markers. 

Mitochondrial  DNA  markers  used  in  this  study  were  RFLPs  in  three  regions.  One 
region  was  within  the  cytochrome  c  oxydase  subunit  I  (CO-I)  gene,  which  contains  a 
Hindi  site  in  west  European  bees  but  not  in  African  and  east  European  bees.  A  second 
region  includes  portions  of  the  cytochrome  c  oxydase  subunit  I  and  II  (CO-II)  gene  which 
has  length  polymorphisms  in  African  and  west  European  samples.  A  third  region  is 
within  the  cytochrome  b  region  which  contains  a  Bgffl  site  present  in  east  and  west 
European  samples  but  absent  in  African  samples.  Primers  and  amplification  conditions 
used  were  as  described  in  Hall  and  Smith  (1991)  and  in  Lee  and  Hall  (1996)  without 
further  modification. 

Nuclear  DNA  markers  used  were  those  reported  earlier  in  this  dissertation 
(Chapter  2,4  and  5)  and  summarized  in  Table  7-2.  Amplification  conditions,  restriction 
enzyme  digestions,  and  electrophoresis  were  performed  as  described  in  Chapters  2, 4,  and 
5.  In  addition  to  the  markers  reported  here,  three  loci  were  used:  locus  123 1  region  E1N 
digested  with  Hhal  (Hall,  1998),  a  section  of  locus  130  digested  with  ,4/wI,  and  a  section 
of  locus  178  digested  with  Mspl ,  which  revealed  three  African-specific  alleles,  one  east 
European-specific  allele,  and  one  west  European-predominant  allele,  respectively. 
Primer  sequences  and  amplification  conditions  for  locus  1231  were  as  described  in  Hall 
(1998)  without  further  modification.  Loci  130  and  178  were  amplified  in  50mM  Tris- 
HC1  pH  9.2, 16mM  (NH4)2S04  and  1.5mM  MgCl2  with  200uM  of  each  of  the  four 
deoxynucleotides  (dNTPs),  1  unit  of  Taq  DNA  polymerase,  and  250nM  of  each  primers 
in  a  total  volume  of  25  ul.  Primers  used  for  130  were  [27]  5'- 


TCCGGTAC AAAGATAACTAACG-3 '  and  [28]  5 ' -AACGGAGGTTGTTTGAC AC-3 ' . 
Primers  used  for  178  were  [29]  5 '-TGCGTGATAACGTGGAAATAGCGGA-3 '  and  [30] 
5 '-ACGCCGAATCTGCCTTTTATTTGC A-3 ' .  Locus  130  was  amplified  with  the 
following  PCR  profile:  95°C  for  1:00  min,  followed  by  35  cycles  of  94°C  for  45sec, 
58°C  for  45sec,  and  72°C  for  1 :00  min  and  a  final  extension  at  72°C  for  5:00  min.  The 
amplification  profile  for  locus  178  was  as  for  locus  130  but  with  an  annealing 
temperature  of  60°C  and  an  extension  of  72°C  for  2:00  min.  Enzymes  and  reagents  were 
obtained  from  Gibco  BRL  (Lifetechnologies,  Gaithersburg,  MD). 

Allele  Counts  and  Data  Analysis. 

Genotypes  of  all  the  workers  for  each  locus  in  every  population  were  tabulated. 
Allele  frequencies  were  calculated  from  this  table  using  the  computer  software  package 
GENETDC  and  GENEPOP  version  1.2  (Raymond  and  Rousset,  1995).  Frequency  data 
for  loci  1231,  130,  and  178  have  been  already  reported  in  both  Old  and  New  World 
honey  bee  populations  (Hall,  1990;  Hall,  1992b;  Hall  and  McMichael,  1996)  but  were 
obtained  again  in  this  study  with  the  African  and  European  samples  from  the  Old  and 
New  World  previously  described  in  Chapters  2  and  5. 

Fifty-three  alleles  from  eight  loci  (12  locus  enzyme  combinations,  LECs)  were 
scored,  of  which  seventeen  were  African-specific,  ten  were  east  European-specific  and 
five  were  west  European-specific.  Collective  frequencies  of  African,  east  and  west 
European  alleles  were  calculated  for  each  locus  and  averaged  across  loci  to  obtain  mean 
values  for  each  population  (Table  7-2).  Markers  were  tested  in  13  Old  and  New  World 


populations  for  a  total  of  617  individuals  from  322  different  colonies.  Allele  frequency 
data  were  also  used  to  calculate  observed  and  expected  heterozygosity  as  a  measure  of 
genetic  diversity  as  described  (Nei  1978;  Nei,  1987)  using  the  computer  software 
package  GENETTX.  Haplotypic  diversity  values  or  gene  diversity  values  (HdiV),  which 
correspond  to  the  proportion  of  heterozygotes,  were  calculated  as  a  measure  or  genetic 
variation  in  each  population.  These  values  provide  the  probability  of  randomly  selecting 
two  different  alleles  in  a  population  and  were  calculated  as  described  in  Nei,  (1987) 
considering  all  the  alleles  at  all  loci  examined.  Values  of  Haw  will  vary  from  "0"  (no 
diversity,  population  fixed  for  one  allele  at  each  locus,  thus  only  with  individuals 
homozygous  at  every  locus)  to  "1"  (high  diversity,  populations  only  with  individuals 
heterozygous  at  every  locus).  Exact  tests  of  Hardy- Weinberg  equilibrium  were 
performed  using  the  computer  software  package  Arlequin  VI .  1 .  Departures  from  Hardy- 
Weinberg  equilibrium  were  tested  using  a  modified  Markov-chain  algorithm  (Guo  and 
Thompson,  1992).  Table  wide  significance  values  were  tested  and  adjusted  using  a 
sequential  Bonferoni  procedure  (Holm,  1979;  Rice,  1989;  Sokal  and  Rohlf,  1989). 

For  each  bee  in  every  population,  the  genetic  composition  was  calculated  as  the 
proportion  of  African  and  European  specific  alleles  across  loci.  Frequencies  of 
individuals  with  the  same  genetic  composition  were  used  to  determine  the  degree  of 
hybridization  in  each  population.  Hybrid  index  distributions  (Mallet,  1995)  were  also 
used  to  determine  levels  of  hybridization  at  the  population  level.  African-specific  alleles 
were  given  a  value  of  "-1"  and  European  alleles  "+1".  Alleles  common  to  African  and 
European  populations  were  scored  with  a  "0".  A  score  for  each  allele  at  each  locus  was 
assigned  to  every  worker.  Hybrid  index  values  were  obtained  for  each  individual  by 
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summing  the  scores  across  loci.  Frequency  distributions  of  these  scores  were  generated 
for  each  population.  Mean  value  comparisons  across  populations  for  the  different 
coefficients  studied  were  tested  for  significance  using  Duncan  multiple  range  tests  (SAS 
institute).  Genetic  variation  was  determined  using  an  analysis  of  molecular  variance 
(AMOVA),  which  partitions  the  total  genetic  variation  into  a  hierarchical  structure 
(Excoffier  et  al,  1992)  (between  groups,  among  populations  within  groups,  and  within 
populations).  AMOVA  analysis  produced  O  statistics,  analogous  to  Wright's  F's 
statistics,  which  characterized  the  variation  between  groups  (Oct),  between  populations 
within  groups  (Osc)  and  within  populations  (Ost).  Six  groups  were  defined  for  the 
AMOVA  (populations  in  parentheses):  African  (A.m.scutellata),  west  European 
{A.m.mellifera  and  A. m.iberica),  east  European  (A.m.ligustica  and  A. m.carnica),  New 
World  European  (United  States),  New  World  African  (Venezuela,  Costa  Rica,  Honduras, 
and  Mexico)  and  Texas  (Texas  1995  and  Texas  1997).  Population  pairwise  Ost  values  (or 
Fst  values)  were  calculated  as  a  measure  of  population  differentiation.  Significance  of  O 
statistics  and  population  pairwise  Fst  values  were  established  with  permutation  tests  using 
1000  replicates.  Gene  frequency  data  were  also  used  for  cluster  analysis  of  Old  and  New 
World  populations.  From  these  frequencies,  a  genetic  distance  matrix  was  calculated 
using  Cavalli-Sforza  genetic  distance  parameter  in  NTSYS-pc  v  1 .2  (Exceter  software) 
(Cavalli-Sforza  and  Edwards,  1967).  Cluster  analysis  of  the  different  populations  was 
performed  using  the  unweighed  pair  group  method  with  arithmetic  average  (UPGMA) 
and  by  the  neighbour-joining  method  (NT)  (SAHN  and  NJOIN  functions  in  NTSYS-pc). 
To  evaluate  the  goodness  of  fit  of  the  trees  produced  by  the  two  methods,  a  cophenetic 
matrix  was  calculated  from  both  trees  produced  by  the  UPGMA  and  the  NJ  method  using 
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the  function  COPH  in  NTSYS-pc.  These  matrices  were  compared  to  the  genetic  distance 
matrix  used  to  generate  the  trees.  Pearson's  correlation  coefficients  were  calculated  for 
each  matrix  comparison  and  used  as  a  measure  of  goodness  of  fit  (Smouse  et  al,  1986). 
High  values  of  this  correlation  coefficient  correspond  to  better  fits  of  the  data  with  the 
tree  produced.  Graphical  representation  of  the  genetic  relationship  among  populations 
was  also  done  using  a  principal  coordinate  analysis  (PCO)  with  the  functions  DCENTER 
and  EINGENinNTSYS. 

Results 

Mitochondrial  DNA  Allele  Frequencies. 

Frequencies  for  African,  east  European,  and  west  European  mitochondrial  DNA 
in  New  World  African  samples  are  listed  in  Table  7-3.  In  the  U.S.A.  samples,  east  and 
west  European  mitotypes  were  found  in  79  and  21%  of  the  colonies  respectively.  West 
European  mitotype  frequencies  in  New  World  European  populations  were  very  similar  to 
those  found  with  nuclear  DNA  markers,  reflecting  the  contribution  of  the  west  European 
subspecies  in  the  North  American  population  (see  below).  West  European  mitotypes 
were  found  in  all  the  feral  samples  from  the  U.S.A.  In  the  Texas  1995  samples,  93  and 
7%  of  the  colonies  had  African  and  west  European  mitotypes,  respectively.  East  and 
west  European  mitotypes  were  absent  in  all  the  Texas  97  and  the  New  World  African 
populations. 
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Population  Genetic  Parameters:  Allele  Frequencies,  Haplotypic  Diversity,  Hardy- 
Weinberg  Equilibrium 

Frequencies  for  LEC's  1231E1N/ZM,  P130A4M,  and  ?\7S/Mspl  in  the  African 
subspecies  A.m.scutellata,  the  European  subspecies  A.m.mellifera,  A.m.iberica, 
A.m.ligustica,  A.m.carnica  and  A.m.caucasica,  and  the  U.S.A.  are  given  in  Table  7-4. 
Allele  frequencies  for  all  LECs  across  New  World  honey  bee  populations  are  shown  in 
Table  7-5. 

Collective  mean  frequencies  of  African  alleles  increase  in  the  populations  moving 
toward  the  south,  coinciding  with  the  length  of  time  since  African  bee  invasion. 
However,  no  significant  differences  were  found  among  the  collective  mean  frequencies  in 
the  Old  and  New  World  African  samples  including  Texas  1995  and  Texas  1997  (P>0.05) 
(Table  7-6).  African  allele  frequencies  from  Texas  1995  were  also  not  significantly 
different  from  the  U.S. A  population.  Compared  to  the  Old  World  African  samples,  the 
New  World  European  samples  (United  States)  showed  significant  differences  in  the 
collective  mean  frequency  of  African  alleles  (P<0.05). 

Collective  mean  frequencies  of  east  European  alleles  decreased  from  57.9%  in  the 
U.S.A  populations  to  4.2%  in  the  Venezuelan  population,  which  was  the  most 
Africanized  New  World  population  tested  here.  For  east  European  alleles,  mean 
frequency  values  from  the  U.S.A  populations  were  significantly  different  from  the  Old 
World  African  populations,  and  the  two  most  Africanized  populations  in  the  New  World: 
Venezuela  and  Costa  Rica  (P<0.05).  Frequencies  in  the  Old  World  African  population 
were  comparable  to  those  in  Venezuela  and  Costa  Rica.  East  European  allele  frequencies 


in  the  Texas  population  were  comparable  to  those  in  Mexico  and  Honduras  and  were 
significantly  different  from  those  in  Venezuela  and  Costa  Rica. 

Collective  mean  frequencies  of  west  European  alleles  between  New  World 
African  and  New  World  European  populations  were  not  significantly  different  (P>0.05) 
but  were  significantly  different  from  the  Old  World  African  population  in  which  west 
European  alleles  were  found  at  a  mean  frequency  of  1 .5%.  Mean  frequencies  of  west 
European  alleles  were  maintained  relatively  constant  in  New  World  African  and 
European  populations. 

Haplotypic  diversity  values  for  Old  and  New  World  African  and  European 
populations  are  shown  in  Table  7-7.  New  World  African  and  European  populations  had 
significantly  higher  haplotypic  diversity  values  than  the  Old  World  European  populations 
represented  by  east  and  west  European  honey  bee  subspecies  (P<0.05)  but  were  not 
significantly  different  from  the  Old  World  African  population.  No  significant  differences 
were  found  between  the  U.S.A.  populations  and  the  New  World  African  populations. 

No  significant  departures  from  Hardy- Weinberg  equilibrium  were  found  in  any  of 
the  Old  World  European  populations.  For  all  the  remaining  populations,  including  Old 
World  African,  significant  departures  from  Hardy- Weinberg  equilibrium  were  found  in 
36  tests  (23.1%)  from  a  total  of  156  possible  tests  (12  LEC's  x  13  populations). 

Population  Genetic  Structure. 

Results  from  partitioning  the  total  genetic  variation  revealed  by  AMOVA  is 
presented  in  Table  7-8.  Of  the  total  genetic  variation,  69.4%  was  found  within 
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populations,  23.7%  among  groups  and  6.9%  among  populations  within  groups.  This  high 
within  population  variation  results  in  a  high  level  of  population  differentiation  as  revealed 
from  the  Fst  value  (Fst=0.3056  PO.001).  Significant  variation  among  groups  (Fct=0.2365 
PO.001)  and  among  populations  within  groups  (Fsc=0.09  PO.001)  was  also  found. 
Population  differentiation  as  revealed  by  Fst  values  from  pairwise  comparisons  between 
populations  shows  similarities  in  genetic  variation  between  Texas  1 995  and  Mexico,  and 
between  the  populations  from  Honduras,  Costa  Rica,  and  Venezuela  (Table  7-9). 

Cluster  and  Principal  Coordinate  Analyses. 

Dendrograms  produced  from  the  genetic  distance  matrix,  for  all  Old  and  New 
World  honey  bee  populations,  with  the  UPGMA  and  NJ  methods,  are  shown  in  Figure  7- 
2.  Three  major  groups  are  clearly  distinguishable  by  both  methods:  I)  An  African  group 
formed  by  all  the  New  World  African  samples  (T95,  T97,  MEX,  HON,  CRI  and  VEN) 
and  the  Old  World  African  samples  of  A.m.scutellata  (SAF).  In  both  methods,  the 
samples  for  Texas  1995  are  grouped  closer  with  the  samples  from  Mexico.  U)  A 
European  group  formed  by  the  east  European  samples  of  A.m.ligustica  (LIG)  and 
A.m.carnica  (CAR),  by  the  New  World  European  population  (USA),  and  by  the  near-east 
samples  of  A.m.caucasica  (CAU).  HI)  A  west  European  group  formed  by  the  samples  of 
A.m.mellifera  (MEL)  and  A.m.iberica  (IBE).  The  west  European  group  is  separated  and 
clearly  distinguishable  from  the  other  European  group  and  the  African  group  in  the 
UPGMA  dendrogram  but  clusters  more  closely  to  the  African  group  in  the  NJ 
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dendrogram.  Despite  the  similarities  found  in  both  UPGMA  and  NJ  methods,  the 
UPGMA  clustering  provided  a  better  goodness  of  fit  (r=0.88)  than  the  NJ  clustering 

(r=0.71). 

Principal  coordinate  analysis  of  the  genetic  similarities  of  Old  and  New  World 
populations  is  presented  in  Figure  7-3.  The  first  three  principal  coordinate  axes 
accounted  for  72.1%,  33.2%  and  8.6%  of  the  total  variation.  Clusters  defined  were 
similar  to  those  found  with  UPGMA  and  NJ  dendrograms  and  corresponded  to  the  main 
groups:  west  European,  east  European,  the  near-east  lineage,  New  World  European,  New 
World  African,  and  the  Old  World  African  group.  Texan  populations  were  grouped  with 
the  populations  from  Honduras  and  Mexico  as  in  the  UPGMA  and  NJ  clustering 
methods. 

Hybridization  in  New  World  Honey  Bee  Populations. 

Absolute  hybrid  index  values  in  all  populations  studied  are  listed  in  Table  7-10, 
and  distributions  of  the  corresponding  relative  values  are  shown  in  Figure  7-4.  The  Old 
World  African  population  has  a  range  of  hybrid  index  values  from  -10  to  0  compared  to 
the  east  European  populations  which  has  a  range  from  5  to  17.  The  distribution  of  hybrid 
index  values  in  New  World  European  samples  is  similar  to  that  of  the  east  European 
samples,  whereas  the  distribution  of  hybrid  index  values  in  New  World  African  samples 
is  closer  to  the  Old  World  African  samples.  Samples  from  the  two  most  Africanized 
populations  in  the  New  World  (Venezuela  and  Costa  Rica)  have  distributions  similar  to 
those  of  the  Old  World  African  samples,  whereas  samples  from  Texas  have  distributions 
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that  are  intermediate  between  the  Old  World  African  and  the  New  World  European 
samples,  suggesting  strong  hybridization.  Moving  toward  the  south,  a  shift  of  the  hybrid 
index  distributions  toward  the  African  distribution  is  seen  in  the  populations,  and  most  of 
the  intermediate  hybrid  index  values  disappear.  This  trend  is  also  evident  in  the 
genotypic  composition  of  workers  in  the  different  populations  (Figure  7-5)  and  the 
number  of  hybrid  workers  (Figure  7-6).  A  dramatic  reduction  in  the  frequency  and  type 
of  hybrids  is  evident  in  the  populations  that  had  been  Africanized  for  longer  periods  of 
time.  Frequencies  of  bees  with  predominantly  African  genotypes  and/or  African  and 
west  European  genotypes  are  mainly  found  in  highly  Africanized  populations  in  the  New 
World.  Although  hybrids  with  east  European  genotypes  are  also  found,  their  frequencies 
were  low. 

Discussion 

Frequency  and  Distribution  of  Mitotypes  in  New  World  Populations. 

The  high  levels  of  African  mitochondrial  DNA  found  in  neotropical  populations 
is  consistent  with  previous  studies  (Hall  and  Muralidharan,  1989;  Smith  et  al,  1989). 
Populations  in  the  subtropical-temperate  region  in  Texas  also  showed  high  levels  of 
African  mitochondrial  DNA,  pointing  to  the  continuous  expansion  of  the  African  bee  as 
unbroken  maternal  lineages,  that  is  as  swarms  with  African  queens  (Hall  and 
Muralidharan,  1989). 
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Nuclear  DNA  RFLP  Allele  Frequency  Distributions  in  the  New  World 

Asymmetric  gene  flow  from  feral  African  to  managed  European  neotropical 
honey  bee  colonies  was  reported  from  previous  studies  using  markers  specific  to  African 
and  east  European  groups  of  subspecies  (Hall,  1990).  European  markers  were  present 
only  at  low  levels  in  feral  African  swarms  near  the  expanding  front  and  were  virtually 
absent  behind  the  front  in  highly  Africanized  populations.  Managed  European  colonies 
behind  the  front  had  a  greatly  reduced  level  of  European  alleles,  apparently  due  to 
considerable  African  paternal  gene  flow.  In  the  current  study,  only  feral  swarms  were 
tested,  and  levels  of  introgression  into  European  colonies  were  not  measured.  Therefore, 
possible  asymmetric  gene  flow  was  not  investigated.  However,  a  more  accurate  measure 
of  the  composition  of  the  feral  population  was  obtained.  More  samples  and  populations 
were  tested  with  a  larger  collection  of  markers  with  broadened  specificity. 

Earlier  studies  used  east  European  nDNA  markers  and  did  not  have  nDNA 
markers  specific  to  west  European  bees  (Hall,  1990).  Thus,  any  west  European  paternal 
gene  flow  could  not  be  detected.  In  a  more  recent  study,  a  west  European  contribution  to 
the  feral  African  population  was  suggested  (McMichael  and  Hall,  1996).  This  current 
study  included  west  European-specific  or  predominant  markers.  In  populations  from 
tropical  and  temperate  regions  in  the  New  World,  west  European  alleles  were  found  at 
constant  frequencies  even  in  highly  Africanized  populations.  The  results  indicate  that 
west  European  and  African  bees  do  not  exist  as  separate  populations.  Virtually  all  west 
European  alleles  in  the  neotropical  populations  were  found  associated  with  African 
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alleles  in  individual  bees,  pointing  to  admixture  between  these  two  subspecies.  The 
original  honey  bee  population  in  Brazil  was  largely  of  the  subspecies  A.m.mellifera 
(Kent,  1988;  Lobo  et  ai,  1989;  Lobo  and  Krieger,  1992).  The  west  European  alleles  in 
feral  neotropical  African  populations  may  be  due  to  the  predominance  of  this  subspecies 
in  the  area  of  Brazil  at  the  time  African  bees  were  introduced.  The  nearly  constant  levels 
of  west  European  markers  in  feral  neotropical  African  populations  suggest  that  these 
markers  had  introgressed  early  when  the  African  population  became  established  and  have 
persisted  as  the  population  expanded. 

In  European  populations  of  North  America  (U.S.A.),  east  European  alleles  were 
found  at  a  70%  frequency.  West  European  alleles  were  found  at  a  20%  frequency.  One 
A.m.caucasica  specific  allele  and  some  African  alleles  were  also  found  but  at  very  low 
frequencies.  This  genetic  contribution  from  different  subspecies  in  the  North  American 
population  is  probably  a  reflection  of  previous  importations  of  honey  bee  subspecies, 
possibly  including  non  -scutellata  African  subspecies  such  as  A.m.lamarckii  (Oertel, 
1976;  Pellet,  1938;  Schiff  and  Sheppard,  1993). 

Just  south  of  the  expanding  African  bee  front  in  Texas,  all  the  feral  honey  bee 
colonies  had  African  mitochondrial  DNA.  In  these  colonies,  African  nuclear  DNA  allele 
frequencies  were  high,  and  east  European  allele  frequencies  were  low.  East  European 
allele  frequencies  increased  from  south  to  north.  The  frequencies  of  east  European  alleles 
and  levels  of  hybridization  were  highest  in  the  feral  African  swarms  caught  in  Texas. 

These  results  are  consistent  with  what  had  been  found  previously,  although 
generally  a  higher  level  of  east  European  alleles  were  found  in  neotropical  bees  than  that 
suggested  earlier.  The  increase  in  frequency  of  the  east  European  markers  from  south  to 


north  may  reflect  ongoing  European  paternal  introgression  (east  European  drones  from 
apiaries  mating  with  feral  African  queens)  and  accumulation  of  European  alleles  into  the 
feral  African  population  as  it  expanded  and  encountered  additional  managed  European 
colonies.  This  introgression  could  have  been  enhanced  by  continuous  importation  of 
European  queens  into  commercial  apiaries  and  migratory  beekeeping  practices,  giving 
the  impression  that  African-east  European  hybrids  persist  in  the  neotropics  at  higher 
levels  than  they  actually  would  under  natural  conditions.  Alternatively,  the  lower  levels 
of  east  European  alleles  in  the  southern  populations  may  reflect  their  loss  over  time  as  the 
African  bees  became  established.  An  apparent  loss  could  be  due  to  continued  migration 
of  African  bees  from  more  Africanized  areas  to  the  south.  A  loss  could  be  due  to 
selection  against  hybrids  between  African  and  European  bees  because  of  poor  adaptation 
to  tropical  environments  or  other  deficiencies.  A  loss  of  European  alleles  would  coincide 
with  changes  seen  in  morphology  (Boreham  and  Roubik,  1987),  in  allozyme  frequencies 
(Taylor  et  al,  1991),  and  in  DNA  allele  frequencies  between  Texas  1995  and  Texas  1997 
reported  in  this  chapter.  As  African  bees  move  into  an  area,  frequencies  of  African 
alleles  should  increase  and  those  of  European  alleles  decrease,  reaching  an  equilibrium 
point  (assuming  random  mating)  which  would  be  affected  by  continuous  introduction  of 
European  bees,  continuous  African  bee  migration  from  the  south,  the  fitness  of  hybrids. 

As  the  expanding  African  bee  front  reaches  its  northernmost  limit,  a  hybrid  zone 
may  be  formed  and  maintained,  as  in  northern  Argentina  where  African  bees  have 
reached  their  southernmost  limit  in  the  Americas  (Taylor,  1977;  Taylor  and  Spivak, 
1984).  Conditions  will  eventually  favor  European  genotypes  as  the  African  bees  reach 
temperate  regions,  and  hybrids  may  also  be  better  adapted  to  the  environment.  The  high 
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degree  of  hybridization  found  in  Southern  Texas  and  the  prevalence  of  hybrids  in  a 
region  that  has  been  Africanized  for  five  years  is  probably  an  indication  that  the  African 
migratory  front  is  reaching  its  limits  in  North  America.  However,  this  has  to  be 
determined  from  further  monitoring  of  the  region  over  a  period  of  several  years. 

Several  mechanisms  have  been  suggested  to  explain  the  low  levels  of  fitness  of 
hybrids.  These  include  differences  in  metabolic  capacity,  perhaps  due  to 
incompatibilities  between  mitochondrial  and  nuclear  DNA  genes  (Harrison  and  Hall, 
1993).  Selection  against  European  genotypes  and  hybrids  might  be  the  result  of 
reproductive  isolating  mechanisms  that  include  both  behavioral  and  physiological 
differences  and  a  better  fitness  of  African  bees  to  tropical  environments.  Such 
mechanisms  include  assortative  mating  (preferential  mating  of  drones  of  one  type  with 
queens  of  the  same  type),  queen  developmental  time  (daughter  queens  with  African 
paternity  emerge  sooner  than  those  with  European  paternity),  and  kin  selection  (workers 
raise  queens  of  their  own  genotype)  (DiGrandi  Hoffman  et  al,  1996;  Kerr  and  Bueno, 
1970;  Page  et  al,  1989;  Winston  et  al,  1983).  Queen  developmental  time  and/or  kin 
selection  might  select  against  European  genotypes  and  would  be  an  effective  mechanism 
to  eliminate  European  alleles  from  the  reproductive  population  that  may  be  expressed  in 
the  workers.  Hybrids  have  been  found  in  neotropical  conditions,  suggesting  a  lack  of 
selection  against  them  (Rinderer  et  al,  1991;  Sheppard  et  al,  1991).  However,  and  as 
explained  previously,  these  hybrids  could  be  the  result  of  constant  European  gene  flow 
from  commercial  beekeeping  operations  into  the  feral  African  swarms.  In  this  study,  a 
high  level  of  hybridization  has  also  been  found  in  the  subtropical-temperate  region  in 
Texas,  which  would  be  expected  if  African  bees  enter  an  area  densely  populated  with 
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commercial  European  colonies.  This  however,  does  not  necessarily  reflect  the 
adaptability  and  survival  of  these  hybrids,  although  higher  hybrid  fitness  is  expected  in 
more  temperate  regions.  In  highly  Africanized  populations  in  the  neotropics,  hybrids  are 
found  in  very  low  frequencies,  and  most  correspond  to  the  African-west  European  type. 

Further  DNA  studies  should  help  resolve  whether  or  not  various  processes  affect 
hybrid  formation  and  survival.  For  example,  the  availability  of  multiple  markers,  more 
easily  analyzable  with  PCR,  enables  parental  analysis.  Genotypes  from  workers  can  be 
used  to  determine  that  of  their  queen  mother  and  the  haplotypes  of  their  drone  fathers. 
Hybrid  maternity  or  paternity  would  suggest  hybridization  events  that  happened  in 
previous  generations.  Hybrid  paternal  drones  would  come  from  hybrid  queens,  pointing 
to  hybridization  events  over  at  least  the  previous  two  generations.  The  absence  of 
parental  hybrids  and  the  presence  of  hybrid  workers  would  suggest  recent  hybridization. 
Paternity  or  maternity  not  found  in  expected  proportions  might  point  to  preferential 
selection  against  hybrid  drones  or  queens,  thereby  pointing  to  the  processes  that  might  be 
responsible. 

The  genetic  interactions  between  African  and  European  bees  can  be  measured 
with  gametic  phase  disequilibria  (random  association  of  gametes).  Gametic  phase 
disequilibria  values  approach  "0"  if  populations  are  randomly  mating.  The  absence  of 
gametic  phase  disequilibria  therefore  would  suggest  admixture  between  populations, 
whereas  the  presence  would  indicate  strong  interactions  of  parental  types  and  possibly 
selection  against  hybrids.  Gametic  phase  disequilibria  may  be  an  effective  measure  to 
detect  selection  due  to  factors  other  than  the  environment,  such  as  hybrid  dysfunction  or 
other  isolating  mechanisms. 
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Table  7-1.  Number  of  colonies  and  individuals  used  in  Old  and  New  World  honey  bee 


populations. 


Population 

No.  of 
individuals 

Number  of 
colonies 

Number  of 
individuals  per 
colony 

Date  collected 

Old  World  African 

A.m.scutellata 

136 

67 

2-3 

Jan- 1990 

New  World  European 

United  States 

69 

44 

1-3 

New  World  African 

Texas  1995 

42 

14 

3 

Jul-Sep  1995 

Texas  1997 

48 

16 

3 

Anr-Mav  1 997 

Mexico 

38 

19 

2-3 

Jan  1988 

Honduras 

96 

31 

2-4 

Nov  1989 

Costa  Rica 

23 

14 

2-3 

May  1989 

Venezuela 

48 

19 

2-3 

Jan-Dec  1988 

West  European 

A.m.mellifera 

13 

9 

1-2 

A.m.iberica 

15 

15 

1 

East  European 

A.m.ligustica 

42 

41 

1-2 

A.m.carnica 

26 

26 

1 

Other  European 

A.m.caucasica 

21 

7 

2-3 

Totals  617  322  1-4 
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Table  7-2.  Nuclear  DNA  markers  used  to  study  hybridization  and  genetic  diversity  in 
New  World  honey  bee  populations.  Number  of  alleles  per  locus  and  their  corresponding 
frequencies  (in  parentheses)  are  shown.  


Locus/Enzyme 

Total 
number 
of  alleles 

African(a) 

East  European^ 

West  European(c) 

XPS-3A4M 

7 

4  (0.6288) 

1  (1.0000) 

XPL-l/AluI 

2 

1  (0.9497) 

XPL-2/Aval 

6 

2(0.5114) 

2  (0.8738) 

XPL-2/Haem 

7 

2(0.1708) 

2(0.8521) 

XPL-5/Haem 

5 

3(0.1148) 

1  (0.9547) 

XPL-5/Ddel 

4 

1  (0.2037) 

1  (0.3119) 

1  (1.0000) 

XPL-5/Spel 

2 

1  (0.5321) 

P227-S3xt/y4/«I 

7 

2  (0.5280) 

1  (0.8960) 

P227-S3xt/#/#I 

5 

1  (0.6140) 

1  (1.0000) 

miElWHhdl 

4 

3  (0.6767) 

?\30/AluI 

2 

1  (0.9378) 

Pm/Mspl 

2 

1  (0.9121) 

Totals0 

53 

17(0.4049) 

10(0.7988) 

5(0.8888) 

(a)  :  Collective  frequency  of  African-specific  alleles. 

(b)  :  Average  collective  frequencies  in  the  east  European  populations  A.m.ligustica  and 
A.m.carnica.  Values  were  calculated  by  averaging  the  collective  frequencies  of  east 
European-specific  and  east  European  predominant  alleles  in  both  populations. 

(c)  :  Average  collective  frequencies  in  the  west  European  populations  A.m.mellifera  and 
A.m.iberica.  Values  were  calculated  as  in  the  east  European  populations. 

(d)  :  Values  in  parentheses  correspond  to  the  averages  across  loci. 


Table  7-3.  Frequencies  of  mitochondrial  DNA  types  in  New  World  honey  bee 


populations 

Mitotype 

Population 

Number  of 
colonies 

African 

East  European 

West  European 

United  States 

44 

0 

0.79 

0.21 

Texas  1995 

14 

0.93 

0.07 

Texas  1997 

16 

1.00 

Mexico 

19 

1.00 

Honduras 

31 

1.00 

Costa  Rica 

14 

1.00 

Venezuela 

19 

1.00 
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Table  7-6.  Average  collective  frequencies  of  African,  east  and  west  European  alleles  in 
New  World  and  South  African  honey  bee  populations.  Means  with  the  same  letters  are 

not  significantly  different.  

Allele  specificity^ 


Population  African  (n=7)  East  (n=8)  West  (n=5) 


United  States 

0.006 

± 

0.000 

b 

0.579 

± 

0.034  a 

0.276  ± 

0.006 

a 

Texas  1995 

0.255 

± 

0.045 

ab 

0.194 

± 

0.011  a 

0.210  ± 

0.005 

a 

Texas  1997 

0.329 

± 

0.047 

a 

0.208 

± 

0.013  a 

0.179  ± 

0.010 

a 

Mexico 

0.298 

± 

0.054 

a 

0.203 

± 

0.007  a 

0.179  ± 

0.006 

a 

Honduras 

0.339 

± 

0.068 

a 

0.123 

± 

0.005  ac 

0.151  ± 

0.005 

a 

Costa  Rica 

0.361 

± 

0.090 

a 

0.069 

± 

0.002  cd 

0.161  ± 

0.003 

a 

Venezuela 

0.345 

± 

0.079 

a 

0.042 

± 

0.002  cd 

0.170  ± 

0.005 

a 

South  Africa 

0.403 

± 

0.048 

a 

0.012 

± 

0.000  d 

0.016  ± 

0.000 

b 

Average  0.292  ±  0.054  0.179  ±0.009  0.168  ±0.004 


(1):  n=number  of  locus  enzyme  combinations  (LEC's) 
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Table  7-7.  Haplotypic  diversity  values  in  Old  and  New  World  honey  bee  populations. 
Means  with  the  same  letter  are  not  significantly  different.  


Population 

n(1) 

Haplotypic  diversity  (HdiV)  ±  s.e 

Costa  Rica 

12 

0.5296 

± 

0.0313  a 

Texas  1995 

12 

0.5279 

± 

0.0376  a 

Texas  1997 

12 

0.5136 

± 

0.0330  a 

Mexico 

12 

0.5097 

± 

0.0396  a 

Honduras 

12 

0.4722 

± 

0.0427  a 

United  States 

12 

0.4396 

± 

0.0359  a 

Venezuela 

12 

0.4101 

± 

0.0706  a 

A.m.scutellata 

12 

0.3484 

± 

0.0988  ad 

A.m.carnica 

12 

0.2296 

± 

0.0340  cd 

A.m.ligustica 

12 

0.1735 

± 

0.0337  cd 

A.m.caucasica 

12 

0.1590 

± 

0.0450  c 

A.m.mellifera 

12 

0.1020 

± 

0.0293  c 

A.m.iberica 

12 

0.1006 

± 

0.0315  c 

(1):  n=number  of  locus  enzyme  combinations  (LEC's). 
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Figure  7-1.  Collective  mean  frequencies  of  African  (A)  East  (•)  and  west  (■)  European 
alleles  in  New  World  African  and  European  populations  and  in  Old  World  African 
populations. 


Figure  7-2.  Cluster  analyses  of  Old  and  New  World  honey  bee  populations.  A)  Analysis 
using  the  unweighted  paired  group  mean  arithmetic  (UPGMA)  and  B)  analysis  using  the 
neighbour  joining  method  (NJ).  Trees  were  constructed  using  the  Cavalli-Sforza  genetic 
distance  coefficient.  OWA=01d  World  African,  NWA=New  World  African,  NWE=New 
World  European,  EE=east  European,  WE=west  European  and  0=Other  (near-east 
lineage:  A.m.caucasica).  Groups  represented  are:  African  (I),  east  and  near-east 
Europeans  (IT)  and  west  European  (HI). 
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Figure  7-3.  Three  dimensional  representation  of  principal  coordinate  analysis  of  New 
and  Old  World  honey  bee  populations  as  inferred  from  the  genetic  distance  matrix  using 
the  Cavalli-Sforza  distance  parameter.  Population  abbreviations  are  as  in  Table  7-9. 


Figure  7-4.  Hybrid  index  distributions.  A)  A.m.scutellata  (dark  bars)  and  United  States 
(light  bars),  B)  Texas  1995,  C)  Texas  1997,  D)  Mexico  (light  bars)  and  Honduras  (dark 
bars),  E)  Costa  Rica  (light  bars)  and  Venezuela  (dark  bars),  and  F)  Old  World  European: 
A.m.mellifera  (hatched  bars),  A.m.iberica  (open  bars),  A.m.ligustica  (crossed  bars), 
A.m.carnica  (dark  bars)  and  A.m.caucasica  (light  bars). 


Hybrid  index 


Figure  7-5.  Bubble  plots  representing  the  genotypic  composition  of  Old  and  New  World 
honey  bee  populations.  Frequencies  of  workers  with  a  specific  genotypic  composition 
are  in  proportion  to  the  size  of  the  bubble.  New  World  African  populations  (Texas  95  to 
Venezuela)  are  arranged  according  to  a  north  to  south  geographical  location  coinciding 
with  increasing  times  since  the  African  bees  invasion. 
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Figure  7-6.  Frequencies  of  African-European  hybrid  workers  in  New  World  honey  bee 
populations.  African  (A),  African-east  European  (O),  African- west  European  (■)  and, 
African-east-west  European  (+)  hybrid  distributions  correspond  to  the  proportion  of 
workers  with  a  corresponding  hybrid  type  relative  to  the  nuclear  DNA  genotypes. 


CHAPTER  8 
CONCLUDING  REMARKS 


The  arrival  of  the  African  honey  bee  to  the  United  States  and  its  dispersal  to  most 
of  the  southern  states  will  have  a  negative  impact  in  the  beekeeping  industry  in  this 
country.  Behavioral  characteristics  including  high  swarming,  absconding,  and  stinging 
make  African  bees  difficult  to  manage  (Michener  et  ah,  1972;  Michener,  1975). 
Commercial  beekeeping  operations  with  African  bees  will  require  different  methods  of 
handling  and  managing  colonies  and  will  be  subject  to  restrictions  on  locations  and 
movement  of  beehives.  Alternatively,  beekeepers  will  have  to  rely  on  certified,  desirable 
breeding  stocks.  Although  income  from  honey  production  amounts  to  only  150  million 
dollars  a  year,  the  most  important  value  of  beekeeping  in  the  United  States  is  in 
pollination,  which  directly  and  indirectly  represent  benefits  estimated  at  10  to  20  billion 
dollars  a  year  (McDowell,  1984;  Levin,  1983). 

Accurate  identification  of  African  and  European  bees  is  necessary  for  research, 
regulation  and  commercial  beekeeping.  Prior  identification  has  been  limited  to 
morphometries,  cuticular  hydrocarbons,  allozymes,  and  mitochondrial  DNA.  These 
characteristics,  however,  have  their  limitations  as  was  explained  in  previous  chapters. 
With  previous  genetic  markers,  hybridization  studies  have  been  limited  due  to  the  low 
levels  of  polymorphisms  detected  (allozymes)  and  the  uniparental  mode  of  inheritance 
(mitochondrial  DNA).  To  understand  better  the  processes  involved  in  the  spread  of  the 
African  bees  and  the  genetic  interactions  between  African  and  European  subspecies, 
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nuclear  genetic  markers  with  high  levels  of  polymorphism,  that  are  biparentally  inherited, 
are  needed.  Nuclear  DNA  RFLP  markers  specific  to  different  subspecies  of  honey  bees 
including  east  European  (A.m.ligustica  and  A.m.carnica),  west  European  {A.m.iberica 
and  A.m.mellifera),  near-eastern  (A.m.caucasica)  and  African  (A.m.scutellata)  have  been 
sought  using  different  methodologies.  RFLP  markers  were  initially  found  using 
laborious  procedures,  that  is  with  Southern  blots  and  randomly  cloned  honey  bee  DNA 
probes.  RFLP  markers  first  found  with  probes  were  made  analyzable  by  standard  PCR 
procedures  (Hall,  1998),  which  do  not  require  high  quality  and  quantities  of  DNA.  Thus, 
large  number  of  samples  could  be  tested  at  a  relatively  low  cost. 

As  reported  in  later  chapters,  markers  were  found  in  DNA  first  amplified  with 
PCR.  To  enhance  this  search,  markers  were  found  in  DNA  amplified  by  long  PCR  and 
were  subsequently  made  analyzable  by  standard  PCR.  The  nuclear  DNA  PCR-RFLP 
markers  specific  to  groups  of  honey  bee  subspecies  described  in  this  dissertation,  make 
for  an  important  addition  to  those  that  already  exist  or  that  are  still  being  converted  to  a 
PCR  analyzable  format.  Despite  the  collection  of  markers  now  available,  more  markers 
specific  to  African  and  west  European  bees  will  be  needed  in  the  future.  Most  of  the 
markers  now  available  are  specific  for  the  subspecies  groups,  but  no  single  allele  is 
diagnostic  (Ayala  and  Powell,  1972)  that  is,  specific  and  nearly  fixed  in  a  population. 
African  and  west  European-specific  markers  are  difficult  to  find  due  to  the  high  genetic 
diversity  in  the  African  population,  and  to  the  effects  of  genetic  pollution  that  result  from 
migratory  beekeeping  practices  and  importation  of  queens  from  other  environments. 

The  markers  specific  to  different  honey  bee  subspecies  were  successfully 
employed  in  this  research  to  characterize  the  subspecies  composition  and  the  extent  of 
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hybridization  in  New  World  populations  including  a  region  in  southern  Texas  close  to  the 
predicted  hybrid  zone.  Evidence  of  a  high  degree  of  Africanization  in  New  World 
populations  was  obtained,  however,  a  high  degree  of  hybridization  was  also  apparent 
during  the  first  years  after  African  bee  invasion.  An  apparent  lost  of  east  European 
alleles  was  seen,  greater  in  populations  farther  south  that  had  been  Africanized  for  longer 
periods  of  time,  suggesting  selection  against  bees  with  east  European  genes.  West 
European  allele  frequencies  however,  remained  relatively  constant,  suggesting  the 
possibility  of  an  admixed  African- west  European  population. 

These  markers  will  continue  to  be  valuable  for  ongoing  studies  in  the  neotropics, 
for  studies  of  the  hybrid  zone  expected  to  form  in  the  southern  portion  of  the  United 
States,  and  for  studies  of  African  gene  introgression  north  of  the  hybrid  zone  where 
European  bees  are  predominantly  found.  Future  research  results  may  point  to 
mechanisms  that  might  isolate  the  African  and  European  subspecies  and  may  reveal  if  the 
apparent  loss  of  hybrids  in  the  neotropics  reflects  poor  adaptation  to  the  environment  or  if 
other  hybrid  deficiencies  are  involved.  Whatever  the  cause,  hybrid  survival  is  an 
important  factor  to  consider  in  hybridization  programs  proposed  to  minimize  the  impact 
of  Africanization  in  the  United  States,  for  example,  genetic  dilution  of  feral  African  bees 
or  the  development  of  African-European  hybrids  for  commercial  beekeeping.  Also,  such 
studies  of  African-European  interactions  should  provide  a  better  understanding  of  the 
evolutionary  processes  involved  in  subspecies  differentiation  resulting  from  their 
geographical  distribution. 
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